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PREFACE 



This paper constitutes a deliverable to Task T-R2-597.01, "SDI Battle Management/ 
C3 Studies'' in accordance with Section 5.0, "Schedule," of the task order dated 1 October 
1988, It presents an overview of tracking methods and issues writen for the nonexpert, 
and an overview of tracking algorithm architectures. We summarize the algorithms 
surveyed to facilitate the understanding of survey responses and underscore the algorithms* 
general features, information flow, and calculational techniques. 

The paper endeavors to serve both the nonexpert and the expert. For the nonexpert, 
the survey of methods explains tracking problems and solutions to provide both a sunmiary 
of the state of the art and the technical background for the algorithm survey. For the 
expert, the survey of algorithms provides a catalog of approaches and results in tracking 
and a community of designers with whom to interact. 
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EXECUTIVE SUMMARY 



A. ALGORITHM SURVEY 

Table 1 provides one summary of the 19 algorithms surveyed in this report The 
general information about an algorithm that we felt should be immediately available 
includes: What phase of a ballistic missile's trajectory does the algorithm track?; Does the 
algorithm use single or multiple sensors?; Does it track clusters of objects as well as 
individual objects?; Does the algorithm rely on track information being handed over by 
another source? and What is the status of the algorithm? Other information of a more 
technical nature is provided in Table 2: In what manner does the algorithm share 
information among multiple sensors?; Does the algorithm perform measurement-to-track 
association by assignment or multiple hypotheses? For a detailed description of these temais 
and summary of algorithms please see Chapters 2-4. ^ 

The tables illustrate that tracking algorithm activities, at least those we surveyed, are 
concentrated in boost and midcourse phases, mosdy multiple sensor tracking of individual 
objects. In addition, a few organizations are developing algorithms for tracking closely 
spaced objects or clusters. The manner in which the tracking information is processed 
varies, as can be seen from the different types of algorithm architectures implemented. 
Last, and most important, both major approaches to the association problem are being 
addressed. 

B. CRITICAL ISSUES! 

SDI tracking algorithms will face demonstration and evaluation milestones in the 
near future. Solutions to the most difficult tracking problems considered both in isolation 
and as part of a surveillance system remain to be successfully demonstrated. In particular, 
the future activities in SDI tracking algorithms must focus on five critical issues: 



! This section draws heavily on discussions by the SDI Panels on Tracking and from Dnimmond, OB., 
•'Multiple Target Tracking Lecture Notes," 18 March 1988, Technology Training Corporation, 
Torrance, CA. 
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• Quster tracking 

• Scan-to-scan association 

• Perfonnance evaluation and prediction methods 

• Testbed 

• Signal processing. 

!• Cluster Tracking 

Qust^ tracking methods are currendy under development at a few organizations, as 
demonstrated by Table 3. 



Table 3. Cluster Tracking Activities 



Organization 


Description of Alaorithm 


Advanced Systems Architectures 


Tracks spawned from a common source are combined to 
form a cluster track. Cluster tracks are also formed from 
objects with similar tracks. No mention of unresolved 
objects. 


Hughes Aircraft Company 


State-of-the-art multiple sensor tracking of unresolved or 
resolved clusters of objects. Group-to-object transifion 
included. 


Lockheed Missiles & Space Company 


S;tate-of-the-art multiple sensor tracking of unresolved or 
resolved clusters of objects. Group-to-object transition 
included. 


MIT Lincoln Laboratory 


Individual object tracks for members of a group of closely 
spaced objects (CSO) are initiated from a cluster track that 
is generated by the edges of the cluster. 


MindGate Technologies 


Pattern matchirig of clusters of co-moving objects. 


Space Computer Corporation 


Track-before-detect approach for determining clusters in 
velocity sj3ace. Potentially very useful as a method for 
performirig computationally affordable individual object 
track initiation. 


TRW 


Cluster tracking to implement a pattem-matching track 
initiation alaorithm. 



The many conceptual difficulties in approaches to cluster tracking are just beginning 
to be addressed Much more work in this area is needed. However promising they may 
be, none of the algorithms listed in Table 3 have been demonstrated to be a solution to the 
cluster tracking problem. The SDI Panels on Tracking are actively pursuing this issue. 
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Individual object tracking will be impossible, too difficult, or too expensive to be 
practical, or not necessary during portions of the ballistic missile defense engagement. 
Shortly after deployment from post-boost vehicles, reentry vehicles and decoys may be so 
closely spaced as to be unresolvable to tracking sensors presenting as point or extended 
objects. Even when resolved, RVs and decoys could be so closely spaced that the 
computational resources required for individual object tracking would be prohibitively 
expensive because of scan-to-scan association difficulties. A sensor's ability to resolve 
objects depends on more than its optical qualities and signal processing. Resolution is also 
a function of the viewing geometry and range. For this reason, as the sensors move along 
their orbits, the closely spaced objects may unresolve as well as resolve. Therefore, cluster 
tracking will play a very important role in any tracking algorithm architecture. 

There are critical operational requirements for maintaining tracks on individual 
targets, including discrimination of RVs from decoys, threat assessment, and attack 
execution. As the threat resolves, cluster tracks spawn individual object tracks, that is, 
individual object tracks are initialized from the cluster track. Therefore, cluster tracking 
should be evaluated, in part, on the quality of the initial estimates for the spawned 
individual object tracks. 

2. Scan-to-Scan Association 

There have been no full-scale demonstrations and evaluations of the two principal 
competing conceptual approaches to this problem. Much more work is needed in this area. 

The great challenge of SDI tracking results from the high density of target and 
clutter observations reporting out of the sensor's signal processor. A high density means 
that the association of tracks to measurements cannot be made without significant 
uncertainty or error. At large but achieveable computational cost, tracks can be assigned to 
one scan's worth of observations in an optimal fashion. Incorrect assignments can lead the 
tracking system into estimatmg the^uality of track predictions as better than they actually 
are. Since track prediction affects the assignment of weapons to targets and the ability of 
weapons to autonomously locate their targets, poor performance in this regard would have 
critical implications for the management of tiie ballistic missile defense engagement. 
Misassociations may also result in the loss of track as the filter follows an incorrect 
sequence of observations. 

Misassignments that may occur in a high-density environment may result in tracks 
based on measurements from more than one target Impure tracks over the course of n:iany 
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scans randomly mix measurements from several targets. In this case, any phenomenology 
used to discriminate between classes of objects, such as RVs and decoys, that depends on 
repeated measurements will be of limited usefulness. In other words, poor track purity 
performance because of misassignments liniits multiple scan discrimination techniques. 

Assignment algorithms are computationally affordable but may not provide the 
necessary performance. On the other hand, multiple hypothesis algorithms should provide 
superior performance but their computational requirements may not be affordable. 
Rigorous, full-scale testing is critical 

3. Performance Evaluation and Prediction Methods 

Much more work is needed in developing methodologies for fairly scoring tracking 
algorithms that use different conceptual and mathematical approaches. Until recentiy, there 
was no such scoring methodology. One scoring methodology has been agreement 
developed by the SDI Panel on Tracking Parameters. Development of appropriate scoring 
criteria must keep pace with new approaches to tracking. 

The complexity of complete SDI tracking algorithms is such that analytic track 
performance predictions do not exist. More work must be done in this area. The only 
alternative is to run computer simulations, which are costiy and sometimes difficult to 
interpret. 

Much more work also must be done on determining the required computational 
resources of tracking algorithms before expensive simulations are run. Computer 
throughput and memory demands will play a critical role in selecting tracking algorithms. 

4 . Testbed 

There is a critical need for a portable testbed that can be used in the development of 
tracking algorithms, not just for evaluation purposes. Contractors are naturally reluctant to 
bring their algorithms to a central testbed facility during development to avoid revealing 
proprietary details and embarrassing algoritiim performance, which is to be expected during 
development 

The portable testbed should consist of a complete set of library modules for those 
functions that support a tracking algorithm. Library modules would consist of accepted, 
standardized models for such things as possible threat trajectories and signatures, 
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background signatures, and signal processing.^ The tracking algorithm designer should be 
able to select individual modules to plug into his sensor system in support of developing a 
tracking algorithm, 

A modest step in this direction is threat data developed by members of the SDI 
Panel on Track Parameters that has been distributed to 13 organizations by IDA, A second 
generation of threat data is nearly ready for distribution. An approved signal processing 
model could be used to supply measurements for the tracking algorithm. In particular, 
realistic, accepted models of background clutter in which tracking algoridims are expected 
to operate need to be made available as soon as possible, 

Currentiy, ekch algorithm design team must provide their own testbeds. This is 
done at great cost The government pays for the redundant effort of contractors developing 
their own testbeds with their own supporting models. 

The efforts of some contractors could be gready handicapped by poor fidelity 
supporting models. This could, periiaps, prevent good ideas from getting adequate testing. 
In this case, the government would suffer the opportunity cost of losing a good idea. 

The government pays an additional cost in that low fidelity support function models 
used by tracking algorithm designers may lead to misleading performance results that might 
not be discovered until late in the development cycle and at the great expense of full-scale 
computer simulations. 

A portable testbed is critical for the ability to run full-scale simulations during 
development to provide tracking algorithm designers an opportunity to study the problenis 
they face and to adequately exercise different approaches. 

S. Signal Processing 

There may be large potential for tremendous growth in tracking performance by 
better or new, innovative approaches to signal processing. For instance, if a signal 
processor could greatly decrease the high density of observations passed to the SDI 
tracking system by removing stars, persistent background, and decoys, the processing load 
for a space surveillance and tracking system would be substantially lessened. Although 
such an approach has been hypothesized, we have seen no evidence of such a bulk filter. 



2 In addition, modules with alternative levels of fidelity will expand the usefutaess of the testbed to 
various purposes and stages in the development pnx:ess. 
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1. INTRODUCTION 



The Institute for Defense Analyses (IDA) was tasked by the Strategic Defense 
Initiative Organization (SDIO) to monitor, evaluate, and facilitate the development of 
tracking algorithms. Among other things, IDA was asked to survey tracking algorithms 
under development for or applicable to SDI in order to ascertain the status of activities in 
this critical area. This paper reports the results of ID A*s survey. 

As part of IDA's overall task, three working panels were established to provide a 
forum for addressing problems in tracking. These panels,^ which meet bimonthly for three 
consecutive days, are staffed by tracking algorithm designers from Federal Contract 
Research Centers, many companies, and each of SDFs sensor program elements. The SDI 
Panel on Critical Issues in Tracking developed a common survey format for describing 
tracking algorithms to ensure that key questions were answered in a succinct manner and to 
simplify the process of understanding the details of the activities. One part of that format 
was a decomposition of all tracking algorithms into architectures consisting of four-track 
initiation and four track maintenance generic processing chains. It was expected that all 
responses would conform to the survey format and algorithm architectures. 

To a large extent, the accomplishment of the survey was dependent on the tracking 
algorithm community responding to ID A*s request for information. Given that replying to 
such a survey is not a contractual obligation, the number of answers received is gratifying 
and IDA expresses its appreciation to those who took the time to carefully and thoughtfully 
respond. This survey report does not contain algorithms for some SDI development efforts 
because of classification or proprietary restrictions. 

The survey benefited greatly from the collective expertise of the SDI Panels on 
Tracking. A list of the contributors from these panels is included in Appendix C. The 
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The SDI Panel on Tracking Parameters, the SDI Panel on Critical Issues in Tracking, and the SDI 
Panel on Advanced Concepts. For a descriptipn of their activities see the Proceedings of the SDI 
Panels on Tracking. 
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survey also benefited greatly from three useful expositions on multiple-target tracking.'* 
The reader desiring more depth on this subject is directed there. 

The paper consists of an executive summary, four chapters, and three appendices. 
Chapter 2 presents an overview of tracking methods written for the non-expert. Chapter 3 
presents an overview of tracking algorithm architectures developed by the SDI Panels on 
Tracking. Chapter 4 contains sunmiaries of the surveyed algorithms that are intended to 
facilitate the understanding of the unprocessed responses and underscore the algorithms* 
general features, information flow, and calculational techniques. The survey form and the 
* responses can be found in the appendices. Our conclusions and recommendations are 
contained in the Executive Summary. 

The paper endeavors to serve both the nonexpert and the expert. For the nonexpert, 
the survey of methods explains tracking problems and solutions to provide both a summary 
of the state of the art and the technical background for the algorithm survey. For the 
expert, the survey of algorithms provides a catalog of approaches and results in tracking 
and a community of designers with whom to interact. 



Samuel S. Blackman, Multiple-Target Tracking with Radar Applications, Artech House Inc., 1986; 
OE, Drummond, "Multiple Target Tracking Lecture Notes." UCLA October 1985, Revised 
18 March 1988, Tfcchndlogy Training Corporation. Torrance, CA.; and Yaakov Bar-Shalom and 
Thomas E. Fortmann, Tracking and Data Association, Academic Press, Inc. 



1-2 



11. OVERVIEW OF TRACKING METHODS 



Before summarizing the results of the algorithm survey, it is worthwhile to provide 
the nonexpert reader with an overview of tracking methods. We begin by explaining the 
tracking process in a general, introductory fashion to acquaint the reader with teraiinology 
and fundamental concepts. The difficulties of SDI tracking are described next. Finally, the 
principal tracking approaches are described in depth but at a level appropriate for the 
nonexpert. 

A. INTRODUCTION 

1 . What is Tracking? 

A track is an estimate, based on sensor measurements, of the kinematic components 
(position, velocity, and acceleration) that describe the motion of an object. These 
components are collectively referred to as the state of the object; the state's evolution in 
time describes the object's motion. Therefore, a track is an estimate of the object's state 
derived from sensor measurements of it.^ 

The tracking process as generally practiced today consists of the interrelated 
functions of association and estimation. A (radar, optical, acoustic) sensor system 
uses its observations of an object's reflections or emissions to derive measurements of 
the object's state. The portion of the state that is measured depends on the sensor type. 
For instance, a passive optical observation cannot provide range measurements because it 
relies entirely on emissions. In general, the measurements are some (not necessarily linear) 
function of the state. 

Association is the decision process of Unking observations or tracks of a 
common origin. Links can be made observation to observation, observation to track, or 
track to track. Observations taken at (nearly) the same time by multiple platforms or from 
. one platform's different sensor systems can be linked together as assumed to have a 
common origin for the purpose of sensor fusion. 
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The state could also include quaniiiies other than kinematic components, such as temperature. 
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A sensor's observations can be linked across frames^ to form a time sequence of 
measurements. This sequence without further processing is a type of track, one that lacks 
the ability to predict the future and is limited to the measurements rather than the full state. 
A time sequence of linked observations, however, can be processed, that is, statistically 
filtered, to transform a measurement sequence into an estimate of the full state's evolution 
as a function of time. Statistical filtering is the estimation part of tracking. 

In the next section we will discuss in detail the Kalman filter, which currently is the 
most general and powerful track estimation method commonly used. Here it suffices to 
note that the Kahnan filter refers to the algorithm that produces the statistical estimate of the 
state and covariance. 

A major task in tracking algorithms is the linking of observations to tracks. A gate 
is a region in the sensor's field of view, determined in part by the prediction from the 
Kalman filter, where the subsequent track measurement is likely to fall. As a rule, only 
observations in the gate are considered for association with that track, thus greatiy reducing 
the number of computations. Often in SDI applications there will be more than one 
observation within a gate and, therefore, several possible observation-track pairings. Since 
the output from a Kalman filter participates in determining the size of the gate, the 
association and estiniation functions are interrelated. We will discuss this further below. 

Last, as with observations, tracks from multiple sensors can be linked as assumed 
to have a common origin, also for the purpose of sensor fusion. 

2. What are the Difficulties in SDI Tracking? 

The essential difficulties in SDI tracking are tiie large number of objects to be 
tracked, tiie high density of observations, and die inability of sensors to resolve individual 
objects from closely spaced neighbors. 

The massive number of objects that have to be tracked in SDI scenarios requires 
huge computing resources. To understand this, consider the computational burden from 
processing one Kalman filter for each object being tracked. Since all the tracks update on 
each frame, every few seconds the information processor must perform die necessary, 
highly nontrivial, update calculations of die Kalman filter. We see tiiat tius computer 



^ A frame is defined as one data collection survey of the surveillance region. In this form, the definition 
is independent of whether the sensor surveUs by mechanically sweeping the field of view with detectors 
or surveils elecUDnically with staring detectors. 
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burden scales linearly with the number of tracks. Faster, more efficient Kalman filters can 
mitigate this huge computer load Unfortunately, this is the least of the difficulties. 

In common SDI scenarios, the density of observations originating from targets and 
clutter is likely to be high. Recall from our discussion above, the observation-to-track 
association function establishes a gate in the sensor's field of view around where it expects 
to find the observation for the track. A high density means that there is likely to be more 
than one observation in a gate. When this occurs, the tracking system cannot know with 
certainty which observation if any originates with the target. Handling this association 
problem is the most computationally intensive aspect of tracking. Incorrect 
observation-to-track association can lead to poor track performance, loss of 
track, and tracking errors far worse in reality than those predicted by the 
Kalman filter. This is the most critical SDI tracking difficulty. 

One approach to managing the high density threat is to temporarily forego tracking 
individual targets and rely instead on tracking the group or cluster in which they are 
traveling. There are substantial computational advantages to this approach. A track can be 
established for some carefully chosen parameters of the group, such as the groiip centroid 
and extent. This saves on the computer resources required when many targets are so close 
that it is not practical to process many individual tracks and circumvents the problem of 
large numbers of misassociations likely to occur if the individual tracks were maintained. 

Unresolvable closely spaced objects cause another major difficulty. The sensor's 
ability to resolve neighboring objects, of course, depends on the sensor, the range, and the 
viewing geometry. If unresolved, a group of closely spaced objects may appear as a 
relatively large (compared to the signal from individual objects) extended object on the 
sensor's detectors. A track, however, can be established on the extended object 

Another major tracking problem caused by CSOs is that the resolution can be 
unstable from frame to frame, e.g., two targets may be resolved on one frame and not 
resolved on the next and possibly resolved on the following. The instability of the 
measurements stresses the association and estimation processes. 



B. ESTIMATION, ASSOCIATION, AND DECISION 



1. Estimation: The Kalman Filter^ 

A statistical estimation filter consists of three parts: models for the dynamics and 
measurement processes, statistical assumptions, and an optimality criterion. The result is 
an algorithm for transforming measurements of a state into an estimate of it. The model for 
the objea's dynamics describes its time evolution, which may contain elements, referred to 
as process noise, that are unknown or unpredictable, except for their statistics. One 
example is the random changes in acceleration typical of rocket boosters. A measurement 
model must both specify the relation between the state and the measurements and account 
for the generally random inaccuracies, referred to as measurement noise, always present in 
measurements. For instance, we may wish to determine the object's position, velocity, and 
acceleration from noisy position measurements alone. The most simplifying assumption is 
that the system noises are statistically independent, white processes. 

Optimality criteria establish a measure for the "goodness" of an estimate. Once 
selected, an optimality criterion, such as maxindzing the likelihood function orminimizing 
the mean square error between truth and estimate, leads to a procedure for transforming a 
set of measurements into an estimate for the state. When the random processes in the filter 
are assumed gaussian, or we require the estimate to be linear in the data, all optimality 
criteria lead to the same estimator. 

The term filter can be thought of as a generic term for the process of recovering 
information from noisy measurements. Statistical filtering develops information from 
noisy measurements by assuming that the desired signal and unwanted noise can be 
distinguished by their statistical properties. 

The term filter also refers to a type of infonnation processing that is distinct from 
two related types, smoothing and prediction. Filtering means the recovery at some 
particular time, tk, of information about the system using measurements up to and including 
that time. Smoothing differs in that tiie information about the system need not become 
available at tk, and measurements derived later than tk can be used in obtaining information 



This secuon is based on material in the foUowing pubbcauons; Andrew P. Sage and J™fs L Melsa. 
Estimation Theory With Applications to Communications and Control, McGraw-HUl Book Company. 
1971; Brian D.O. Anderson and John B. Moore, Optimal Filtering, Prenuce-Hall Inc., IV /y. 
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about the system at tk. Prediction is the forecasting type of information processing in 
which the aim is to obtain information at tk about the state of the system at a later time- 

The Kalman filter is the mathematically optimal estimator for deriving at some 
particular time; tk, an estimate of the state and its covariance from measurements of the 
states. 

The filtered estimate of the state is processed in two stages: a time update and a 
measurement update. The first step, the time update, is the filtered estimate of the state at 
the preceding time, tk-i, predicted ahead one step. The measurement update involves the 
difference between the associated measurement and the predicted measurement, referred to 
as the innovation, multiplied by the Kalman gain. Calculated by the filter, the gain 
determines the weight given to the new measurement information. The prediction of the 
state at some future time is computed from the present filtered estimate, without employing 
the innovation process. 

The filter also provides for the time and measurement update of the covariance of 
the state estimate. In addition, the algorithm calculates the covariance of the innovation 
process, which is the measurement prediction uncertainty, indicating the quality of the 
prediction. 

2 . Association: The Concept and Role of a Gate 

The innovation chi-square is derived from the innovation and its covariance 
matrix. It specifies an elliptical volume in measurement space known as a gate that is an 
indication of die track prediction uncertainty. The gate establishes an acceptance or 
validation region into which observations considered for association widi the track must 
fall. Landing within the gate is a necessary but not sufficient condition for an observation 
to be considered as having originated from tiie track because incorrect measurements may 
also fall in die gate or the target may not have been detected in die gate. The sole purpose 
of the gate is to decrease the processing load by decreasing the number of possible 
observation-to-track association pairs by limiting the number of candidates. Observations 
within the gate are often called validated 

The size of the gate is determined by fixing the probability that the correct 
measurement will fall witfun its volume. Since this probability is set to less than one there 
is a nonzero probability tiiat tiie correct observation will not be a candidate for association: 
A larger ellipse enjoys a higher probability of capturing die measurement that originates 
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with the target but at the expense of possibly increasing the number of association 
candidates and thereby increasing the processing burden. 

If a gate is empty, there is no measurement to update that particular track. In this 
case, the Kalman filter, which predicted the center of this gate, predicts ahead an additional 
frame but now without the benefit of measurement update information. The measurement 
prediction uncatainty is increased; thus, the gate is larger for the next frame. 

There is some measurement origin uncertainty even if there is one measurement in 
the gate because the target may not be detected in the gate and the observation could be 
from a different target or a false signal. When updated, the measurement prediction 
uncertainty of the Kalman filter typically decreases regardless of the true source of the 
updating measurement because it fails to account for its origin uncertainty. In other words, 
the measurement prediction uncertainty typically decreases after the update regardless of the 
origin of the measurement. 

For this reason, misassociations unless compensated for cause the filter to estimate 
the quality of its prediction as better than actual. Misassociations may also cause poor track 
performance, that is, a large gap between estimated and true tracks, and loss of track as the 
filter follows an incorrect sequence of observations. Discrimination and battie management 
is also adversely affected. 

3* Decision: Strategies for Managing Measurement Origin Uncertainty 

The fundamental multiple target tracking dilemma is deciding which measurement to 
use, if any, in updating a track. An observation has three possible sources: an individual 
target, clutter, or a set of unresolved targets. Following Bar-Shalom and Fortmann,^ we 
define clutter as observations from background stationary objects, interference, 
environmental anomalies, false alarms, etc., that are generally random in number, location, 
and intensity. 

The difficulties caused by not knowing the source of a measurement are shown by 
considering tiie tracking situations in Figure 1. In the single-target case, there are two 
measurements witfun the gate, each, or none, possibly originating from the target The 
correct decision is not obvious. In the multiple-target case, the gates overlap, with 
measurement 3 falling in both. It is usually assumed tiiat one measurement cannot be 



8 Yaakov Bar-Shalom and Thomas E. Fortmann, Tracking and Data Association, Academic Press. Inc. 
p. 153. 
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simultaneously assigned to more than one track, that is, each observation is uniquely 
assigned. If we decide that measurement 3 is associated with track 2, then measurement 1 
is the only candidate for association with track 1. Thus, associations over multiple 
targets are interdependent. 

The manner in which the problem of uncertain measurement origin is managed can 
be used to categorize approaches to multiple target tracking. Assignment methods make a 
definitive decision, typically at each frame, on the origin of the measurement. One 
measurement (or none) from those within the gate is selected as having originated from the 
target. 

Instead of selecting a single observation, the probabilistic data association (PDA) 
approaches avoid selecting by averaging over all decisions. The track is updated, using all 
measurements in the gate, weighted by the probability that they are correct. This is referred 
to by Blackman as an all-neighbors approach. 

Multiple hypothesis tracking (MHT) algorithms defer a decision on the origin of the 
measurement. Multiple alternatives are retained as distinct tracks until later information 
improves the probability of the correct measurement-to-track association. 

C. MULTIPLE TARGET TRACKING APPROACHES: ASSIGNMENT, 
PDA, MHT 

1 • Assignment 

The simplest association scheme, known as the nearest-neighbor algorithm, assigns 
a track to the measurement "nearest" to the predicted measurement, such as calculated by 
the innovation chi-square. All other observations in the gate are disregarded. 

One approach to multiple target tracking is to run a nearest-neighbor algorithm for 
each track, independent of all other tracks. This is referred to as uncoordinated nearest 
neighbor. In the situation shown in Figure lb, the processor would decide for track 1 
whether to update with measurement 1 or 3, or nothing, independent of track 2 
assignments. 

We have commented that this is unsatisfactory for high-density environments for 
the reason that associations over multiple targets are interdependent when we require 
unique assignments. This requirement is fundamental to simplifying the calculation 
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because its effect is to make the different states statistically independent Therefore, a 
nearest-neighbor algorithm is generally executed in a coordinated manner as follows. 

A cost matrix is defined by all possible track assignments, including that the 
measurement is from a hew target or false alarm, and all possible candidate measurements, 
including the case that the correct track measurement is not detected in the gate. The cost 
matrix entries are proportional to the probabilities of the assignments. These can involve 
the innovation chi-square for the measurement-track association pair, the probability of 
detection, the probability of the gate, the probability that the observation is fit>m a new 
« source, and the probability of choosing no observation for association with the track. 

An algorithm, such as the Munkres algorithm, is run that assigns measurements to 
tracks in a coordinated fashion by maximizing the sum of matrix entries subject to the 
constraints that no track is updated by more than one measurement and one measurement is 
not assigned to more than one track. The results of such an algorithm is a unique pairing of 
tracks to observations. 

Assignment algorithms are used not only to associate one list of observations with 
one list of tracks as just described. They can also be used to associate two lists of 
observations or two lists of tracks. Furthemiore, there are assignment algorithms tiiat can 
be used to associate data among more than two data lists, for instance linking several 
frames worth of observations to tracks. In this manner, assignment algorithms can 
generate multiple hypotheses, in the sense that more tfian one viable alternative per track is 
retained over a number of frames. The general distinguishing characteristic, however, of 
assignment algorithms is that of a definitive decision. For examples and discussion of 
multiple frame assignment algorithms see the Alphatech survey. 

2. Probabilistic Data Association^ 

The Probabilistic Data Association Filter (PDAF) applies to a single target and is 
strictiy a method for handling tiie problem of multiple observations witiiin the gate of an 
established track. The fundamental idea is to exploit die association probabilities of 
the complete set of observations within the gate for the target 

In the PDAF, each observation in turn is considered as originating from the target 
Also, die case tiiat the observation originating from the target is not detected is considered. 



9 This section bonows heavily from Bar-Shalom and Fortmann, ibid. 
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An association hypothesis consists of assigning one observation (or none) to the track and 
considering all others as statistically independent clutter. 

The PDAF procedure first multiplies the probability of each association hypothesis 
with the updated state estimate that assumes that hypothesis is true; thereby weighting that 
state estimate. Then tiie final PDAF state estimate is formed as the weighted average, that 
is, the sum over the association hypotheses of the weighted average of the state estimates 
for each hypodiesis. The estimate can be shown to be equal to the predicted state plus the 
weighted average of the individual innovation chi-squares multiplied by the standard 
« Kahnan gain. 

The covariancc of the final PDAF state estimate follows immediately as the average 
over the covariances for each hypothesis. This can be seen to be equal to a sum of three 
terms. One term is the prediction covariance multiplied by the probability that no 
observation originated with the target A second term consists of the covariance of the state 
updated with the correct measurement multiplied by the probability that the target-originated 
observation is available. The last term increases the covariance of the total updated state to 
account for the uncertainty in the origin of the observation. 

The Joint Probabilistic Data Association Filter (JPDAF) extends the fundamental 
idea of the PDAF to multiple targets by computing the association probabilities joindy 
across all targets rather than for each individual track. The final JPDA state estimate is 
calculated as before as an average over the association hypotheses. 

To summarize, the PDA state estimate is an average over observation-to-track 
association hypotiieses. Each hypothesis consists of an unique assignment of die track to 
an observation. The total state estimate is an average over the many feasible assignments 
for the one track. 

The five principal distinguishing characteristics of the PDA approaches are the 
assignment of one track to many observations, one per hypothesis, the exploitation of 
association probabilities, the calculation of state estimates as averages over association 
hypotheses, an adjustment to the covariance for source uncertainty, and a lack of organic 
track initiation logic. The association probabilities are calculated v^dth Bayes Theorem from 
Probability Theory. For this reason, PDA approaches are one memb« of a class of 
tracking approaches referred to as Bayesian Tracking Algorithms. We will discuss 
other members of this class shortiy. 
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3. Multiple Hypotheses Tracking 

An intuitive approach to managing multiple observations in a gate is to split the 
original track into many tracks, one for each validated observation. This process is known 
as track splitting. Each track is updated with the associated observation and carried 
forward to the next frame in the standard fashion. The fate of these many tracks depends 
on whether the associated observations arise from clutter or targets. 

If the source of the associated observation is clutter, then subsequent observations 
for this track will be randomly detected and located. The quality of the track, therefore, is 
expected to decrease markedly. For this reason, a pruning mechanism is usually contained 
in track splitting algorithms to drop low-quality tracks. 

An observation originating from targets may have four sources. First, the 
observation could be from the original target. In this manner, this algorithm ensures 
maintaining the original track. "Extra" observations within the gate could be from objects 
just released by a comimon carrier vehicle, such as reentry vehicles and decoys released 
from a post-boost vehicle. Also, the new observations could be newly resolved, closely 
spaced objects from what had been a single unresolved CSO. In these two cases, tracks 
are split, also referred to as spawned, by initializing the new tracks with the original 
track's state estimate on the previous frame. 

Last, the extra observations could be from new targets just moving into detection 
range or field of view that happen to fall within the gate. Track splitting may provide poor 
estimates in this case because the original state estimate may have very litde to do with the 
newly detected targets, except for their location on the sensor focal plane. 

A track splitting algorithm has two limitations. First, the algorithm disregards all 
observations that fall outside the gate. For this reason, a separate track initiation algorithm 
must be included in the battie management system that uses this approach. 

The second and major limitation of the track splitting algoritiim is that association 
over multiple tracks is performed in an uncoordinated manner. There is no conflict 
resolution logic that manages the problem of observations within multiple gates. 

Track splitting algorithms are distinguished by the assignment of multiple 
observations in the gate to one track, tiie deferring of difficult assignment decisions until 
more information is generated, and the absence of association probabilities based on global 
hypotheses. 
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Deferring difficult assignment decisions is prototypical of multiple hypothesis 
tracking algorithms but other MHT approaches use association probabilities based on 
global hypotheses. Recall that PDA algoritiims use the probabilities of the multiple 
association hypotheses, each one consisting of unique assignments of the latest set of 
observations to the tracks from the prior fi^e. 

Both the PDA and track splitting algorithms assign many observations to one track, 
but very differently. Track splitting generates one candidate track for each 
validated observation assignment. PDA yields a final state estimate that is 
averaged over all candidate hypotheses. 

Two avenues for generalization are suggested. First, use association probabilities 
based on global hypotheses but maintain individual state estimates rather than average state 
estimates. Second, extend the association hypotheses over many frames rather than just the 
most recent. An Optimal Bayesian Tracking algorithm would generate hypotheses across 
all frames from the first through to the current 

Bar-Shalom and Fortmann^^ relate the PDA single frame approaches to an optimal 
Bayesian algorithm as follows. Consider a time isequence of observations, one observation 
per frame, from the initial to the present time. Such a sequence forms one possible target 
history, tiiat is, one possible track. Consider all possible such sequences. The set of all 
possible assignments at the current frame can be decomposed into tracks at the previous 
frame associated with some observation from the current frame. A few moments thought 
reveals this is to be an abstract description of track splitting. 

An association probability for each observation sequence, tiiat is, a probability for 
each track, can be calculated, conditioned on the entire set of observations. As in PDA, the 
conditional probability for each hypothesis multiplied by the state estimate tiiat assumes that 
hypotiiesis is true is summed over all possible hypotiieses. Thus, the updated state 
estimate for a track is an average over the diffa-ent possible association hypotheses. 

Optimal PDA associates over all frames, not just the most recent Its computational 
expense may be prohibitive. A suboptimal approach looks back N frames, referred to as 
N-backscan, rather than all the way to the initial frame. The original PDA is the zero- 
backscan suboptimal version. 



10 Ibid. 
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The optimal PDA shares the four principal distinguishing characteristics of the zero- 
backscan PDA approaches. It generates track-oriented hypotheses, that is, every 
observation is considered for association with each track from the previous frame. This is 
the reason for the absence of organic track initation logic: no observation is considered for ~ 
association with a track that did not exist on the previous frame, that is, a new target 

Reid's^^ multiple hypothesis tracking algorithm remedies the absence of organic 
track initiation logic by generating observation-oriented hypotheses. Each 
observation is associated with a false alarm, as a feasible continuation of a previous track, 
or as a new target, in the following manner. 

Start with the hypotheses generated on the previous frame. Consider the first new 
observation. Generate a new hypothesis/or each possible assignment of the observation: 
as a false alarm, as a feasible continuation of a previous track, or as a new target Take this 
new set of hypotheses and repeat this procedure with the second observation, except that 
more than one observation cannot be assigned to one track. Continue in this way until 
every cuirent observation has been assigned. 

For instance, Reid*s algorithm applied to Figure 1 would generate eight hypotheses 
in the single target case and 30 hypotheses in the multiple target case. Reid refers to these 
as cluster hypotheses. See Table 4. 



Table 4. Results of Reid*s Algorithm Applied to Figure la. 



Cluster Hypothesis 


Measurement 1 


Measurement 2 


1 


FA 


FA 


2 


T1 


FA 


3 


T2 


FA 


4 


FA 


T1 


5 


T2 


T1 


6 


FA 


T3 


7 


T1 


T3 


8 


T2 


T3 



Notes: FA « false aiarm; 

T1 m the original track; 

T2 ■ a possible new track originating with obsen^ation 1 ; 
T3 » a possible new track originating with observation 2. 



1 1 Donald B, Reid, "An Algorithm for Tracking Multiple Targets," IEEE Trans. Auto. Control, Vol. 
AC-24, No. 6, December 1989, pp. 843-854. 
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While the total number of cluster hypotheses generated can be quite large, the 
number of track assignment hypotheses is relatively few. This is important because the 
number of computations in this approach can be greatly reduced, as we now demonstrate. 

Observe from Table 4 that a cluster hypothesis consists of one possible set of track 
assignment hypotheses. The original track, Tl, is either associated with measurement 1 or 
2 or with no measurement. New track 2, T2, may or may not be generated by 
measurement 1 and similarly for new track 3 and measurement 2. This yields seven track- 
oriented hypotheses. A similar calculation shows that there are ten track-oriented 
hypotheses in the multiple target case in Figure lb. 

A particular track assignment hypothesis can appear in many different cluster 
hypotheses. Each track assignment hypothesis is followed by a Kalman filter update 
computation. If the track update computations were performed for each cluster hypothesis, 
then the same filter update computation would be repeated many times. Instead, 
association probabilities are calculated over alternative target assignment hypotheses and 
then mapped onto the larger set of cluster hypotheses. It must be reemphasized that each 
association hypothesis assumes unique track assignments so that the -association 
probabilities are calculated over statistically independent states. 

The optimal implementation of Reid's algorithm would require ever-increasing 
computer memory as more hypotheses are generated on each frame. A practical version 
must limit the number of hypotheses. One method is to divide the set of tracks and 
observations into independent groups, which Reid calls clusters, requiring conflict 
resolution. Hypotheses are also limited by pmning and merging. Hypotheses considered 
unlikely, say those below some threshold, are dropped while those that are "similar" 
according to some criteria are combined. These operations are suggestive of track spUtting 
but in that case there were no association probabiUties and there were multiple assignments 
of tracks to observations. The Reid algorithm generates individual state estimates that are 
scored by association probabiUties. It is foraiulated in optimal and suboptimal versions, 
which can be weU implemented. Clustering, pruning, and merging can be adjusted to fit 
hypothesis growth to track density and throughput and memory computer restrictions. 

The major limitation of Reid's algorithm is that it does not include multiple 
assignments of tracks to observations, such as may occur in merged measurements, or 
observations to tracks, such as may occur in track spawning. The fundamental reason for 
this is the manner in which the association probabUities are calculated. One association 
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hypothesis consists of a set of unique assignments. The probability of the assignnient 
hypothesis, calculated by Bayes rule, when the states are statistically independent, 
decomposes into the products of probabilities for the individual tracks. Unique 
assignments under each hypothesis ensure the statistical independence of the states. 

Kovacich of Lockheed Missiles and Space Company recently described a Bayesian 
multiple hypothesis tracking algorithm that remedies the defects in Reid's approach to 
produce the most advanced MHT algorithm to date.^^ The key idea is to use a Bayesian 
network architecture (also known as influence diagrams) to provide a calculus to represent 
and manipulate joint probability distributions such as those that occur in multiple target 
tracking. Rather than decompose the association-to- track problem into assignment 
hypotheses, the fundamental unit in Lockheed's approach is the scene which is defined as 
the joint set of observation-oriented hypotheses, track-oriented hypotheses, and track 
spawning outcomes for different clusters. The probability for each individual possible 
outcome is calculated by the Bayesian network. 

D. GROUP/CLUSTER TRACKING 

Thus far the discussion has been limited to tracking of individual objects. In 
Section 2.1 of this chapter we pointed out that one approach to managing the high density 
SDI threat is to forego tracking individual objects and instead track groups. In this section 
we will describe the issues and methods of such an approach. 

1. Definitions 

To begin, we need to define what we mean by group and cluster. In the previous 
section the term clustering referred to collecting interacting observations and tracks, that is, 
track clusters. In the discussion of group tracking, the term cluster is defined differentiy 
and is used in reference to nearby objects, that is, target clusters. Before stating the 
definition, a comment on what is to be included. Consider two simple examples that maik 
the extremes: (1) a long line of equally spaced objects; (2) a sphere of objects. In the 
second example, each object is within some radius of the center of the sphere. In contrast, 
each object in the first example is within some metric distance of at least one other object 



Michael Kovacich, -Application of Bayesian Networks to Midcourse Multi-Target Tracking," 
presentation to the SDI Panel on Advanced Concepts in Tracking, Proceedings of the SDI Panels on 
Tracking, No. 4, 1989. 
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Another issue that often generates confusion is whether the objects in the collection 
are based on truth or measurement data. 

We will follow the lead of the SDI Panels on Tracking in our definitions.^^ 

Cluster: A maximal collection of objects each of which is within some metric 
distance of at least one other object in the collection. 

Group: To be used at the author's discretion but to include a cluster. 

Group track: Track established to represent a cluster of objects. 

Qump: A single observation arising from two or more objects. 

2. Types of Group Tracking 

Following Drummond^"* we identify individual and group tracking as the endpoints 
of a spectrum: 

• Group tracking without individual target tracks 

• Group tracking with simple individual target tracks 

• Individual target tracks supplemented with simple group information 

• Individual target tracking without group tracks. 

This ordering is suggestive of a logical sequence of operations that might occur in 
midcourse tracking. During deployment of reentry vehicles (RVs) and decoys from post- 
boost vehicles (PBVs), the threat initially consists of closely spaced objects. The objects 
may resolve with increasing time from deployment as the threat cloud disperses. As the 
sensors move along their orbits, however, the resolution of objects is a function of the 
sensor resolution and the viewing geometry and range. For this reason, the threat could be 
resolving or unresolving during the course of the sensor's observations. 

These considerations lead to the conclusion demonstrated below that the 
information processing in group tracking algorithms must begin by deciding, based on 
data, computational, communication, and operational considerations what combination of 
group and individual target tracking to execute. 



Proceedings of the SDI Panels on Tracking. 

Oliver E. Drummond, Hughes Aircraft Company, presentation to the SDIO Panels on Tracking, 
Proceedings of the SDI Panels on Tracking. 
^ ^ This is not to suggest that the threat density will not or cannot be increased later in the flight 
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It is not possible to establish individual object tracks on members of an unresolved 
CSO, Extended object tracking must be executed. Only when the objects in a CSO are 
resolved is individual object tracking possible. The density of observations and the high 
rate that new objects are resolved, however, may compel group tracking as the only 
practical alternative because of great computational expenses in track initiation and 
misassociation. If objects unresolve, the tracking architecture must extrapolate individual 
object tracks or establish group or extended object tracks. 

There are critical operational requirements for maintaining tracks on individual 
targets, including discrimination of RVs from decoys, threat assessment, and attack 
execution. Group tracking is performed when individual object tracking is impossible or 
too expensive. As the threat resolves, group tracks spawn individual object tracks, that is, 
individual object tracks are initialized by the group track. Therefore, group tracking should 
be evaluated based on its relatively inexpensive computation and communication 
requirements and the quality and the processing load required of the initial estimates for the 
spawned individual object tracks. 

3. Single and Multiple Sensor Group Tracking 

Blackman^^ describes a single sensor group tracking algorithm that tracks the group 
centroid position and velocity. A gating logic that is a generalization of the gate for an 
individual object track determines which observations will be considered for updating the 
group track. A conflict resolution logic is required for all observations that satisfy multiple 
group track gates. All observations assigned to a group track are used to form a group 
observation consisting of a measurement centroid and dispersion ellipse. The measurement 
centroid updates the group centroid state in the standard manner of Kalman filtering. 
Tracks for objects splitting off the group are initialized by the group centroid state. 

Drummond, Blackman, and Hell,!'^ have extended these ideas to multiple sensor 
group tracking. The principal difficulty in multiple sensor group tracking is that the size, 
shape, and composition of the group varies from sensor to sensor. For this reason, 
multiple sensor group tracking must have more information than just the location of the 



Samuel S. Blackman, Multiple-Target Tracking with Radar Applications, Artech House, Inc., 1986, 
Chapter 11. 

OE. Drummond, S.S. Blackman, K.C. Hell, "Multiple Sensor Tracking of Clusters and Extended 
Objects," Technical Proceedings 1988 Tri-Service Data Fusion Symposium, Laurel, Maryland, May 
1988. 
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group. Drummond ct aL*s approach is to model the group as an ellipsoid in three 
dimensions. Separate filters are established for the group centroid and the ellipsoid extent 
parameters. 

The group centroid state estimate initializes tracks for objects that split away from 
the group, as before. The ellipsoid extent state estimate permits sensors in different 
locations to associate groups and facilitates handing over group data to other sensor 
systems. 

£• CONCLUSIONS 

We have reviewed the major difficulties and methods in SDI tracking. We have 
chosen to organize the discussion around the critical problem of uncertain measurement 
origin that arises during the association of observations to tracks. This is sometimes 
referred to as the frame-to-frame association problem. We have not discussed other types 
of association, such as observation-to-observation and track-to-track, that are typically 
performed in sensor fusion, even though real-world SDI tracking algorithms involve 
multiple sensors. We have said littie about filtering beyond defining it and developing its 
relationship with observation-to-track association. 

There are many components in a complete tracking system. Before tracks can be 
maintained by association of validated observations and updating they must be initiated. 
Track initiation refers to the formation of the first or initial estimate of the state of an object 
We have found that, for the nonexpert, following the information flow among these many 
components is often one of the most significant impediments to understanding particular 
tracking algorithms. For this reason, the next chapter presents an overview of tracking 
algorithm architectures developed by the SDI Panels on Tracking. 
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III. OVERVIEW OF TRACKING ALGORITHM 
ARCHITECTURES 



The SDI Panels on Tracking have standardized on a small set of algorithm 
architectures as high-level descriptions of the logical flow of information in tracking 
algorithms. It is expected that most algorithms can be decomposed in terms of these 
algorithm architectures, first identified by DmmmondJ^. 

A natural taxonomy of multiple sensor, multiple target tracking algorithms is 
organized by dividing the problem into track initiation and track maintenance algorithm 
architectures for individual objects and clusters. First consider the algorithm architectures 
for individual object track maintenance contained in Figures 2-5, 

Figure 2 establishes the basic functions of single sensor track maintenance 
algorithms. Figure 3 represents an architecture in which individual sensor tracks are. fused 
together. Instead of combining tracks, a frame's worth of measurements from multiple 
sensors could be combined before being filtered, as in Figure 4. Figure 5 differs from 
Figure 3 in that the individual sensors no longer maintain individual tracks. The system 
uses only central track files. 

Figures 6-9 depict individual object track initiation architectures. Figure 6 
establishes the basic functions of single sensor track initiation algorithms. Figures 7, 8, 
and 9 are similar to Figures 3, 4, and 5, demonstrating that tracks or one frame's worth of 
measurements coiild be combined and that the system could use individual sensor tracks or 
centralized tracks. 

Cluster tracking algorithm architectures involve more than the basic individual 
object functions. The data emerging from the signal processor and requirements from the 
battle manager determine the type of tracking to be performed. This is depicted in 
Figure 10. Figures 1 1-18 are essentially equivalent to Figures 2-9. 



^ ^ O.E. Drummond, "Multiple Target Tracking Lecture Notes." UCLA October 1985, Revised 18 March 
1988, Technology Training Coiporaiion, Torrance. CA. 
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Figures 2-18 establish the functions that make up tracking algorithms. The 
descriptions in the figures use terminology that is indicative of but not specific to 
techniques in particular architectures. This was done in order to remain at a sufficiently 
high level of generality so that these architectures were applicable to most algorithms... 
Hence, the terms association, filter, and track were used in their most general sense. 

There are two major areas of detail lacking in these figures. First, specific 
definitions and descriptions of techniques for each function of the process. Second, the 
distribution and location of computer resources for carrying them out. 

The SDI Panels on Tracking have moved to describe some of the specific 
techniques, which are the fundamental tracking algorithms, most of which are described in 
Chapter 2. For instance, Figure 19 summarizes track maintenance functions and Figure 20 
the elements essential to the description of measurement-to-track association. Finally, 
Figure 21 indicates specific algorithms for this purpose. Figures 22 and 23 repeat this for 
track initiation. 

The overview of methods and tracking architectures presented in the last two 
chapters should prepare the nonexpert reader for the algorithm survey in tiie remainder of 
this report. 
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Figure 2. Single Sensor Track Maintenance (Type I) 
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Figure 3- Single Sensor Then Fuse Track Maintenance (Type 11) 
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Figure 4. Sensor-to-Sensor Then Scan-to-Scan Track Maintenance (Type III) 
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Figure 5. Centralized Sensor-by-Sensor Track Maintenance (Type IV) 
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Figure 6. (Cold Start) Single Sensor Track Initiation (Type 1) 
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Figure 7. (Cold Start) Scan-to-Scan then Sensor-to-Sensor Track Initiation 

(Type II) 
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Figure 8- (Cold Start) Sensor-to-Sensor Then Scan-to-Scan Track Initiation 

(Type III) 
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Figure 9. (Cold Start) Centralized Track Initiation (Type IV) 
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Figure 10. Group Tracking Hierarchy 
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Figure 11. Single Sensor Group Track Maintenance (Type I) 
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Figure 12. Single Sensor Then Fuse Group Track Maintenance (Type II) 
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Flaure 13. Sensor-to-Sensor Then Scan-to-Scan Group Track Maintenance 

(Type III) 
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Figure 14. Centralized Sensor-by-Sensor Group Track IMaintenance (Type IV) 
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Figure 15. (Cold Start) Single Sensor Group Track Initiation (Type I) 
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Figure 16. (Cold Start) Scan-to-Scan Then Sensor-to-Sensor Group 

Track Initiation (Type II) 
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Figure 17. (Cold Start) Sensor-to-Sensor Then Scan-to-Scan Group 

Track Initiation (Type III) 
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Figure 18. (Cold Start) Centralized Group Track Initiation (Type IV) 



• MEASUREMENT-TO-TRACK ASSOCIATION 

• FILTERING 

• TRACK-TO-TRACK ASSOCIATION AND FILTERING 

• MEASUREMENT-TO-MEASUREMENT ASSOCIATION 

• TRACK MANAGEMENT 

Figure 19. Tracic (Maintenance Functions 

m-ii 



IS THE ASSOCIATION PERFORMED IN AN INDEPENDENT MANNER OR IN 
A COORDINATED MANNER? 

IS THE ASSOCIATION OVER ONE SCAN OR MULTIPLE SCANS? 

WHAT IS BEING TRACKED? 
OBJECTS 
CLUSTERS 

TO WHAT EXTENT IS THE SIGNAL PROCESSING AND TRACKING 
SEPARATED VERSUS INTEGRATED? 



Figure 20. Measurement-to-Track Association Descriptors 
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Figure 22. Functions and Descriptors of (Cold Start) Track Initiation 
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IV. SURVEY OF TRACKING ALGORITHMS 



This chapter contains a summary of the algorithm survey responses to facilitate the 
understanding of the responses and to succincdy point out the algorithms* general features, 
information flow and calculational techniques- 

The SDI Panel on Critical Issues in Tracking developed a common survey format 
for describing tracking algorithms to ensure that key questions were answered in a succinct 
manner and to simplify the process of understanding the details of the activities, One part 
of that format was a decomposition of all tracking algorithms into a taxonomy consisting of 
four track initiation and four track maintenance generic processing chains. It was expected 
that most responses would conform to the survey format and processing chains. The 
survey forms and responses are contained in Appendix A and B, respectively. 

SUMMARY OF ALGORITHMS 

1 • O rganization : Advanced System Architectures Ltd. 

Algorithm: An Object-Oriented Architecture for Sensor Data Fusion/ 
Tracking in Dense Threat Environments. 

Submitter: Edward J. G. Goodchild 

Description: The algorithm is designed to perform multiple sensor, birth-to- 
death, three-dimensional tracking of individual objects and clusters. Single sensor scan-to- 
scan association and track-forming functions are excluded. Unassociated two-dimensional 
observations are fused across multiple sensors. Tracks are initiated directly into three 
dimensions from three two-dimensional unassociated observations. New tracks are formed 
by track spawning of boost phase tracks and are assigned to members of a cluster. New 
tracks formed by spawning during die post-boost phase are assigned to members of the 
cluster. Cluster tracks are also formed later by association of a number of non-cluster 
member tracks, all having near identical trajectories. Cluster tracks determine rectangular 
gates that are used to partition new observation data. A six-state Kalman filter is used for 
individual target track state estimation and a second six-state Kalman filter is used to 
maintain die cluster track trajectory and extent. Track pruning tests implemented and 
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planned include: consistent and sensible track behavior, behavior of the error covariance, 
and continued updating with new sensor data. The algorithm is in a concept-proving stage. 
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2, Organization:- Alphatech _ 

Algorithm: Multiple Information Set Tracking Correlator (MISTC) 
Submitter: Robert B. Washburn, Jr. 

Description: This work investigated eight different data association 
algorithms focused on the midcourse problem of multiple sensor individual object track 
initiation and track maintenance. Tracks are initiated as two-dimensional tracks until 
multiple sensor associations are made; three-dimensional tracks are then initialized. A 
single scan of observations from multiple sensors (two or three), containing an assumed 
100,000 objects, was partitioned into spatially separated groups assumed to consist of 
about 100 objects. After partitioning, each data group was processed by the 
tracker/association algorithm. Of the eight association algorithms, tiiree were zero-scan, 
pairwise approaches and five were N-scan, multiple hypotheses approaches. The same 
association algorithm performed scan-to-scan and sensor-to-sensor association. Ten 
hypotheses were permitted per target, with each hypothesis scored by likelihood ratios. 
Extended Kalman filters for state estimation assume the targets travel along Keplerian 
trajectories. The algorithms did not handle track spawning. Algorithms have been 
implemented for sequential, off-line processing in FORTRAN. Plans include incorporating 
CSO tracking and resolution into the algorithms and implementing them on different 
parallel processors. 

3 . Organization: Ball Systems Engineering (VERAC, Incorporated) and 

Daniel H. Wagner, Associates 

Algorithm: SDI Midcourse Tracker/Correlator Algoritiun 
Submitter: LanyRlippelli 

Description: This single sensor, midcourse individual object tracking 
algorithm does not perform cold start track initiation. Instead, warm start track initiation of 
a six-state Kalman filter is accomplished by handover of boost phase information. The 
algorithm relies on the assumption tiiat objects tend to cluster into spatially inseparable 
groups to reduce the combinatorial explosion tiiat results from high target density. Gutter 
is assumed removed by the signal processor and not passed to the information processor. 
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Multiple observation-to-track association hypotheses are formed by track splitting and 
scored using a Bayesian approach that takes into account probability of detection, 
probability of false alarm, false and new target densities, and observation-to-track 
association scores. Hypotheses are deleted if their score is below a percentage of the score 
of the best hypothesis. There is also a maximum number of hypotheses that can be saved. 
A second-generation algorithm has been completed and is in the testing phase. The 
algorithm has been installed at NRL and LANL. 

4 . Organization: CALTECH Jet Propulsion Laboratory 

Algorithm: CALTRAX: The Tracking Program for Simulation-88 
Submitter: James Ortolf (Applied Research Associates) 
Developer: Thomas D. Gottschalk 

Description: This algorithm tracks individual objects in the boost and post- 
boost phases fix)m multiple sensors. Single-sensor, two-dimensional tracks are initiated by 
a three-scan batch processor. Two-dimensional tracks are maintained by a four-state 
Kalman filter and track splitting. Tracks are deleted on the occurrence of a single empty 
gate (probability of detection is assumed equal to one) and merged if of conmion history. 
Mature tracks are propagated to a common time and exchanged with a stereo partner sensor 
to determine the three-dimensional state vector in earth-centered inertial coordinates. The 
two sets of tracks are associated by a modified nearest neighbor algorithm to initiate three- 
dimensional tracks, which are used to solve for launch parameters according to a powered 
flight model. Once initialized, launch parameters are updated on subsequent scans by 
means of extended Kalman filters. Individual sensor observations are associated with 
three-dimensional tracks by a global modified nearest-neighbor algorithm. Any 
observations unassociated with a three-dimensional track on one sensor, if part of a mature 
two-dimensional track, are associated with those on the stereo partner sensor to initiate new 
three-dimensional tracks. The design of the three-dimensional tracker is based entirely on 
parameterized trajectories, with all updates of existing tracks done using extended Kalman 
filters for assumed trajectory models. The next generation tracker will perform parameter 
estimations on arbitrary flight models. 
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5 . Organization: ESL 

Algorithm: Tracking Algorithm for Projea Swat 
Submitter: Jack Liu 

Description: The algorithm, produced to evaluate the applicability of 
DARPA's MOSAIC architecture developed by ESL, performs single platform tracking of 
individual objects during midcourse. The survey response provides no specific 
information on track initiation procedures. Track maintenance is accomplished by a single 
target, multiple hypotheses scheme: nearest neighbor observation-to-track association in a 
rectangular gate updates the track; remaining observations in the gate are used to split new 
tracks; and detections within overlapping gates are utilized by all affected tracks. Tracks 
are scored by the log-likelihood function of the detection relative to tiie prediction plus 
maintenance bias. The scores are cumulative. When a set track redundancy is reached, the 
low score tracks are deleted. Tracks that are very close in terms of estimated object state 
and uncertainty are merged in a probabilistic fashion based on their track scores. A six- 
state extended Kalman filter is the track estimator witii an earth gravity noodel selectable up 
to J6. Multiple platform tracking will be addressed in the near fuftire. 

6. Organization: Hughes Aircraft Company 
Algorithm: Multiple Sensor Ouster Tracking 
Submitter: Oliver E. Drummond 

Developers: Oliver E. Drummond and Samuel S. Blackman 

Description: This algorithm, for use in the early midcourse phase, 
accomplishes multiple sensor tracking of multiple clusters. Ouster tracks can be used to 
initiate individual object tracks as the closely spaced objects become resolved, by track 
spawning based on the estimated PB V track. This permits a smooth transition from PBV 
cluster deployment to individual target tracking. The filtering segment of the algoritiim 
estimates the cluster centroid position and velocity in inertial space and the cluster extent, 
the second central moment in inertial space, of the objects in the group. The extent, which 
establishes the cluster size and shape in inertial space, not just relative to a particular 
sensor, is used to determine which observation belongs to which group. Based on the 
projection of the predicted extent on to the field of view of a sensor, a gate is computed for 
a cluster. The filtering is composed of two filters, one for the state of the centroid and 
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another for the extent The filter for the 6-D centroid state is a simplified extended Kalman 
filter, the extent state has 6 elements and the filter is a pseudo-linear filter. Feasibility tests 
have been conducted successfully for cluster tracking under realistic conditions and further 
testing is under way. 

7 . Organization: Lockheed Missiles and Space Company 

Algorithm: SSTS Tracking and Association Algorithm 
Submitter: Michael Kovacich 

Description: This algorithm accomplishes multiple sensor, birth-to-death, 
group-to-object, midcourse tracking of groups, clumps, and objects in the presence of 
clutter, including stars, nuclear redout, and a structured background. Initial coarse three- 
dimensional track estimates are produced by each sensor after four-to-six updates using an 
iterated maximum likelihood passive ranging algoritiim. Precision ECI tracks are initiated 
by multiple sensor triangulation. The track initiation process completes with ±t formation 
of precision tracks. Data association and track maintenance are accomplished by a multiple 
hypotiieses approach known as a Bayesian network architecture. Pruning, merging, and 
clustering are used to control the combinatoric explosion. Hypothesis scoring accounts for 
clutter density, new track density, missed detections, and cumulative chi-squares. 
Thresholds based on a fraction of the best track score are used for track promotion and 
deletion. A variant of tiie A* search algoritiun is used to find likely hypotheses. Tracks are 
merged and the covariance matrix adjusted accordingly, in a manner similar to PDAF. The 
Bayesian network approach includes tiie multiple assignment of observations to tracks and 
the multiple assignment of tracks to observations. An extended Kalman filter performs 
track estimation. The algorithm is currentiy being implemented in ADA in preparation for a 
1989 demonstration. 

8. Organization: M.I.T. Lincoln Laboratory 

Algorithm: A Mid-Course Track Initiation and Maintenance Algorithm 
Submitter: MingJ.Tsai 

Developers: M.J. Tsai, K.P. Dunn, L.C. Youens, and C.B. Chang 

Description: This midcourse, individual object and cluster track initiation and 
track nuiintenance algorithm accomplishes single sensor track initiation of clusters by 
forming track files for edges of clusters. Tracks on cluster members are then initiated by 
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assuming that targets within the same cluster travel in parallel. Tracks initiated by two 
sensors are merged and track maintenance performed in a sensor-by-sensor centralized 
fashion. Tracks are niaintained by an extended Kalman filter and nearest neighbor 
observation-to-track association. The algorithm has been implemented, tested, and ran in a 
number of simulated threat/sensor scenarios. Currendy, it is being integrated with 
discrimination algorithms and radar tracking functions. 

9 . Organization: McDonnell Douglas Space Systems Company 
Algorithm: Integrated Correlation, Track 
Submitter and Developer: Thomas R. Blackburn 

Description: This algorithm was developed for cold start track initiation and 
track maintenance of individual objects in the late midcourse phase of their trajectory, when 
objects are resolvable but not undergoing the influence of the atmosphere. Two sensor, 
sensor-to-sensor, observation association is performed before track initiation in order to 
initiate three-dimensional tracks. Sensor-to-sensor association is performed by matching 
the pseudo-elevation angle generated from interpolated line-of-sight measurements taken 
from two frames of data. The track file is initiated with a square-root information filter. A 
nearest-neighbor observation-to-track association is used in track maintenance to provide 
measurement updates to a six-state Kalman filter. The algoritiun is in the late conceptual 
development stage and has been tested and debugged ranning against threats consisting of 
about 600 objects. 

1 0 • Organization: MindGate Technologies, Inc. 
Algorithm: Ouster Map Tracking 
Submitter and Developer: Lawrence M. Beyl 

Description: The cmx of diis midcourse tracking algoritiim is the supposition 
that there are patterns within a threat that are naturally formed by the objects dispersed from 
the same PB V and tiiat are heading toward the same target. A collection of co-moving 
objects is termed a cluster. The pattern of angle measurements of objects in a cluster can be 
traced from one scan to anotiier by using tiie previous scan's two-dimensional map of the 
cluster as a pattern for die next scan's data associations. Thus, tiiis is a pattern-matching 
algorithm. The interlocking of the cluster to its source data permits each map to be used as 
a filter to remove the cluster's new measurements from die field-of-view for the sensor 
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with the extracted data set then used to replace the existing set as the new cluster map. 
Each time a cluster map is updated, the associated cluster centroid state is updated via a six- 
state extended Kalman filtering, where the data used in the update is a calculated pseudo- 
measurement created from the collection of angle measurements that define the cluster map. 
The transition to individual object tracking is accomplished by fonning an initial state vector 
and covariance for each object within the cluster through the centroid state, the angular 
separation of the objects within the cluster from the centroid, the cluster spread and spread 
rate, and the angular measurement accuracy. Splitting, merging, and other phenomena are 
handled within the individual object tracking environment, but are always restricted to the 
domain of the cluster. A cluster*s centroid state can be estimated initially based on the PBV 
state at the time of deployment fi-om boost phase information. 

11. O r ganization: MITRE Corporation, Bedford 

Algorithm: The MTTRE Experimental Version Prototype (EVP) Ballistic 
Tracker 

Submitter and Developer: J. A. Krajewski 

Description: This algorithm performs individual object tracking during post- 
boost and midcourse using sensor-by-sensor centralized track maintenance. There is no 
cold start track initiation capability. Instead, the algorithm relies on warm start track 
initiation from handover of boost-phase tracks and spawning of RV/decoy tracks (assumed 
distinguishable from PBV tracks). The data are partitioned based on a user defined angular 
distance threshold and a distance matrix calculated for each group. Hungarian and Greedy- 
type algorithms are used for observation-to-track association. For each association pair, 
the distance is tested against a threshold and, when greater, the association is suppressed. 
A six-state extended Kalman filter is used to estimate the Keplerian three-dimensional 
ballistic tracks. The state is updated using two-dimensional observational data from one 
SSTS satellite at a time. All algorithms have been coded in FORTRAN and are being 
written in Ada. Test cases are currendy being run and analyzed 

1 2 • Organization: MTTRE Corporation, Bedford 

Algorithm: The MTTRE Experimental Version Prototype (EVP) Boost- 
Phase Tracker, 

Submitter and Developer: J. H. Latimer 
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Description: This is a boost phase, individual object, two sensor, track 
initiation and track maintenance algorithm that develops three-dimensional position 
estimates by associating observations from the asynchronous sensors before initiating and 
updating tracks. Cold start track initiation is performed from a single tiiree-dimensional 
position estimate and the a priori assumptions of a three-degree-of-freedom, reference 
rocket trajectory model A Newton's method itCTative procedure is used to fit the reference 
trajectory to tiie three-dimensional target position. The value of die reference trajectory 
velocity and acceleration at the point of the fit is used to initialize the track state estimate. 
The data from two sensors are paired by mapping die data from one sensor onto die focal 
plane of the odier. The position estimate-to-track association is performed by a Greedy- 
type algoritimi. Tracks are estimated by a nine-state Kalman filter for die position, 
velocity, and acceleration in Cartesian earth-centered inertial coordinates. The algoridim 
does not handle stars, false alarms, or other forms of stationary clutter. All algoridims 
have been coded in FORTRAN and are being written in Ada. Test cases are currendy 
being run and analyzed. 

13. Organization: MITRE Corporation, Bedford 

Algorithm: The MTTRE Multi-Sensor, Multiple Target Tracker. 
Submitters and Developers: R. Varad and J. T. McKeman 
Description: This is a boost phase, individual object, three sensor, track 
initiation and track maintenance algoridmi diat associates two sets of stereo-associated data. 
Three sensors are divided into two pairs and observations are associated for each pair. 
Then the two pairs of associated observations are fused. Range is determined from the 
common sensor once target Usts from each pair are formed and dien associated based on 
hinge angles, in-plane angles, and estimates of die baseline ranges determined from each 
pair. A track is initiaUzed when association can be obtained in hinge angle, in-plane angle 
and range for data from two consecutive scans. After initialization, rates for hinge angle, 
in-plane angle, and range arc calculated for each track and predictions of target coordinates 
for die next scan are made. The track state consists of position and velocity estimates in 
diree dimensions. The algoridim is fidly designed, developed, and implemented in Pascal 
on a VAK/WUS system. 



14. Organization: Raydieon Company 

Algorithm: BMEWS Phased Array Radar Upgrade 
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Submitter: Fred Damn 

Description: This boost and midcourse single-sensor, cold-start track 
initiation and track maintenance algorithm processes phased array radar measurements of 
range, elevation, azimuth and target amplitude to estimate position and velocity vectors in a 
six-state Kalman filter, A multiple hypothesis track initiation scheme is used with pulse- 
pair track initiation and track maintenance. The observation-to-track association is 
performed by a nearest neighbor chi-square test 

1 5 • Organization: Raytheon Company 

Algorithm: Ground-Based Radar (GBR) 
Submitter: Fred Daum 

Description: This midcourse single- sensor cold or warm start track initiation 
and track maintenance algorithm processes phased array radar measurements of range, 
elevation, azimuth, target amplitude, and phase to estimate position, velocity, and higher 
order rotational dynamics (for discrimination). A multiple hypothesis track initiation 
scheme is used with pulse-pair track initiation and track maintenance. The observation-to- 
track association is performed by a nearest neighbor chi-square test Three Kalman filters 
are maintained: six state, seven state, and sixteen state. 

16. Organization: Space Computer Corporation 

Algorithm: Velocity Filter Algorithm for SDI Detection and Tracking 
Submitter: William J. Jacobi 

Description: Both boost and midcourse applications of this cold start track 
initiation algoritiun have been investigated. Operating with eitiier resolved or unresolved 
objects, the velocity algorithm performs a combination of signal-to-noise enhancement and 
scan-to-scan association functions utilizing a "track-before-detect" approach. A bank of 
filters matched to different vector velocities provides correlated object positions and 
velocities for track initiation. The vector velocities to which the filters are tuned are derived 
from cross-correlation of successive input image firames. When the data contains "velocity 
clusters," the algoritiim is inherendy robust against object proliferation, merging/crossing 
tracks, background clutter, and temporary loss of data. 
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1 7 • Organization: Systems Control Technology 



Algorithm: 



Dynamic Programming Algorithm (DPA) 



Submitter: 



Kenneth Kessler 



Developer: 



YairBamiv 



Description: Another "track-before-detect" approach, the dynamic 
programming algorithm is a practical and feasible alternative to replace exhaustive search 
techniques for detecting and locating entire target trajectories inside a sensor's field of view 
over some time interval. This is accomplished by batch processing data over a small 
number of frames through a bank of matched filters, where each filter represents a single 
possible two-dimensional trajectory. The algorithm produces simultaneous detection and 
two-dimensional tracking of targets because its output consists of detected targets and their 
associated hit strings. The analysis and software development has been ongoing for over 
seven years. 

18. Organization: TITAN Systems 

Algorithm: Knowledge-Based Sensor Fusion (KBSF) 
Submitter: Timothy E. Brockwell 

Description: This is a tracking algorithm to the extent that it is not entirely 
possible to decouple tracking from discrimination. The program was initiated to determine 
whetiier rule-based techniques could be applied to strategic sensor fusion, specifically to 
demonstrate a rule-based approach to multi-sensor discrimination. The goal is to host and 
evaluate competing algorithms, primarily Uiose designed for platform-to-platform 
association and track maintenance, by building a machine tiiat automatically selects die 
"best" algorithm for fusing multiple or single platform, multiple sensor track data. 
Experiments are being conducted with fuzzy techniques for determining "degrees of 
membership" in the set of valid tracks that is maintained by die system, and for uncertainty 
management in general. Currendy, a variation on die Sequential Probability Ratio Test for 
deteraiining "degrees of membership" in die set of valid tracks is employed. 

1 9 . Organization: TRW, Huntsville 

Algorithm: ADOP Scan-to-Scan and Track Algorithm Set 
Submitter: J. T. Lawson 
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Description: This is a midcoiirse, single sensor, individual object and cluster 
cold start track initiation and track maintenance algorithm that uses pattern matching in the 
data association and derives object velocity estimates from three-color correlation 
processing. In the track initiation phase, association is performed over six frames. The 
algorithm partitions the observations in a scan according to a distance separation threshold. 
For each partitioned set, the group centroid and group velocity are calculated, as are the 
predicted centroid for the next scan, a centroid gate, and a group azimuth-elevation extent 
gate. The grouped observations are associated by pattern matching between scan 1 and 
scan 2 and between scan 2 and scan 3; thereafter the association is performed in an 
independent, nearest-neighbor manner. Objects are assumed to follow non-maneuvering 
Keplerian ballistic trajectories. During the angle-only phase of tracking, a cubic-fit least 
squares filter in azimuth and elevation is used. A precision track mode uses a six-state 
Kalman filter. An iterative batch filter is used to transition between the modes: an initial 
state estimate is obtained by the Gauss algorithm for orbit determination given three angle- 
only sightings; the batch filter is then applied over the previous observation history to start 
the Kalman filter. A track is declared lost and no longer maintained after two consecutive 
empty gates. All clutter is assumed removed by the signal processor and not passed to the 
tracker. Versions of die algorithm's software (Pascal) are available at die USASDC 
Advanced Research Center for both the Honeywell and DEC VAX configurations. No 
significant improvements have been made since 1985. 
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APPENDIX A 
ALGORITHM SURVEY FORM 



TRACKING ALGORITHM SUMMARY 



Submitter's Name: 



Submitter's Company/Organization: 
Submitter's Phone: 



Address of Submitter 



Author's Company/Organization: 
Author's Phone: * 



Address of Author 



TITLE OF ALGORTIHM: 

SPONSOR: 

DEVELOPER: 



(Complete as appropriate) 
NOTES AND INSTRUCTIONS 

(The Summary/Abstract of the algorithm should be limited to half a page. The basic 
answers to the questions below should be limited to three pages. Additional information, 
refa-enced to the section numbers, should be included in a separate appendix that should be limited 
to four pages. Classified information should be included in a separate supplement. The total 
information including Abstract/Summary, basic answers, appendix, and classified supplement 
should not exceed eight pages. If some of the algorithms or details are proprietary, indicate what is 
proprietary and discuss only the nonproprietary aspects.) 

ABSTRACT/SUMMARY 

(Describe, in general terms, how the algorithm works and, if applicable, provide a high- 
level flow diagram.) 

1. CONTEXT 

(What processing chain characterizes your aljgorithm? [see Rgs. 1-4]. What functions 
within the chain are covered by your algorithm, e.g., Urack Initiation, Track Maintenance, etc? To 
what phase(s) of SDI is your algorithm applicable - Boost, Post-Boost, Midcourse, Terminal? 
What are your inputs, e.g:, single sensor or multiple sensors, clustered data or not clustered data, 
etc.? What are your outputs, e.g., 2-d tracks or 3-d tracks, launch parameters?) 

2. NOTABLE FEATURES 

(List any features that distinguish the algorithm.) 
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3. SENSOR ARCHITECTURE AND THREAT SCENARIO 

(Characterize the scenario that drives the algorithm design. What is the threat size and 
density? What are the assumptions about background clutter? Is all the clutter removed prior to 
tracking? What is the constellation size and orbit? What are the measurement errors/biases? What 
is the target kinematic model? Does the target maneuver? Is its signature related to aspect angle?) 

4. SENSOR MODEL/PROCESSING 

(What is the model used for sensor/signal processing? What is the probability of detection, 
false alarm? What is the clutter density and model for the clutter density, e.g., uniform density in 
regions with a Poisson model for clutter returns? What is the sensor? Is it a scanning or staring 
sensor? What is the frame time? Is it variable? How many wavebands of data are available? 
What is the measurement noise, bias, resolution? How are CSOs modeled? What are the 
attitude/navigation errors/bias assumed? What is the precision of your input data? Specify the 
interface between the sensor/signal processor and the tracker, e.g., [time, az, el, snr, extended 
objea intticator, covariance].) 

5 . TRACK INITIATION 

(Specify the method used to initiate tracks. How does it work? Are tracks initiated as 2-d 
tracks or 3-d tracks? Is cold start initiation performed, i.e., is handover data assumed for all 
tracks? Are tracks processed individually or in a batch or both? How is the initial precision track 
state estimate generated? Is track initiation performed sequentially or in a batch mode or both? 
How long does it take to initiate a track (scans or seconds)? How are cluster tracks initiated? If 
applicable, provide a high-level data flow. If Track Initiation and Track Maintenance are not 
separable in the algorithm design, so indicate and describe the Track Initiation and Maintenance 
algorithms in Section 5 and omit Section 6.) 

6.1 TRACK MAINTENANCE - DATA ASSOCIATION 

(Specify the data association approach for scan to scan, similar sensor platform to platform 
association and dissimilar sensors. Are multiple hypodieses generated? How is the combinatoric 
explosion controlled: pruning, merging, clustering? What is the branching factor in the hypothesis 
tree? How deep are the hypothesis trees? How are hypotheses scored: Bayesian, Likelihood, 
Heuristic, other? What search algorithm is employed to develop track hypotheses? Is an 
assignment algoritiim used: Munkres, Brogan-Lemay, otiier? Does track spawning occur? For 
satellitt;-to-satellite association, how is the resolution difference between the satellites accounted for 
in performing association? How are stars and false alarms [stationary clutter] handled?) 

6.2 TRACK MAINTENANCE - STATE ESTIMATION 

(Specify the approach used to estimate the state vector for the track. What is the coordinate 
system? What is the state vector? Whattypeof filter is used? What assumptions are made in the 
filter, e.g., what is the dynamical model, what dynamics are unmodeled and treated as noise, are 
noise measurements assuitied uncorrelated, is the measurement noise assumed to be constant? Is 
the filter nonlinear, iterated, batch, sequential? Are weighted sums of Gaussians used to generate 
the state vector estimate as in PDAF or JPDAF? How does the filter account for biases due to 
targets becoming reisolved [spawning] or become unresolved [crossing]?) 

6.3 TRACK MAINTENANCE - TRACK PROMOTION/DEMOTION 

(Specify the criteria for promoting, demoting and terminating tracks. Specify the method 
for maintaining the track status. [The track status indicates the degree of confidence that the track 
represents a valid target.] What is the criteria for pruning or promoting a track? What is the 
scoring method, e.g., Bayesian, Likelihood, heuristic measures such as m out of n?) 
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1. TRACK FILE MAINTENANCE 



(Specify the data that is maintained in the track file. Are extended objects, clusters, 
maneuvering objects, complete trajectories, threat corridors, etc. detected and maintained? What 
special data stracture for datafile management have been used?) . 

8. OUTPUT TO BM/C3 AND USERS 

(Specify the interface to users of the tracking data, esp. BM/C3. What is the data in the 
interface? How often is the data sent? What is the reporting criteria? How is the data computed?) 

9. COMPUTATIONAL REQUIREMENTS 

(Specify the throughput and memory requirements of the tracking algorithm. What 
scenario is used? How are throughput and memory measured? What machine/language is used? 
Are the results scalcable? Are the results empirical or theoretical? Are there performance bounds? 
What is the target machine? Is the target machine special or general purpose? Describe the degree 
of parallellism [e.g., 10 processors, 100 processors or more?] and the processor architecture. 
What are the sequential and parallel throughput requirements?) 

10. CURRENT STATUS 

(Describe the current status of the algorithm: conceptual development, design, 
coding/debugging, implementation, testing? Has the algorithm undergone performance 
optimization? Has hardware been optimized to execute the algorithm? What are the future plans? 
Has real data from current sensors been used?) 

11. PERFORMANCE MEASURES AND RESULTS 

(Specify the performance measures used to characterize the tracking algorithm. How are 
the performance measures defined? What scenario(s) were used to generate the performance 
values? What are the results? Are Uiere theoretical performance bounds, e.g., Cramer-Rao 
bounds? What are the performance limits? Under what conditions does the algorithm perform 
poorly?) 

12. REPORTS 

(List report(s) that are relevant to the description and performance of the algorithm: 
Additional data such as equation development should be referenced.) 
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GENERIC ALGORITHM PROCESSING CHAINS 



The panel on Critical Issues in Tracking has determined that 
there are four basic processing chains for track initiation and 
four basic processing chains for track maintenance. These chains 
are pictured below. 

The depictions of processing chains presented here are meant 
to be high level descriptions of the logical flow of information 
in algorithms. It is expected that any algorithm can be defined 
in terms of these chains. 

In order to remain at a sufficiently high level of 
generality so that these chains are applicable to all algorithms^ 
it is important to use terminology that is indicative of but not 
specific to techniques in particular processes. Hence, the terms 
association, filter and track are used in their most general 
sense. 

There are two major areas of detail lacking in these 
depictions. First, specific definitions and descriptions of 
techniques for each part of the process. Second, the 
distribution and location of computer resources. 
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APPENDIX B 
ALGORITHM SURVEY RESPONSES 



TRACKING ALQCRITHM SUmART 



SUBMITrED BY: Edward J. G. Goodchild 7»» February 1989 

Advanced System Architectures Ltd. 
Johnson House, 
73 " 79 P&rk Street, 
Camberley, 
Sxarrey, 

England, GUIS 3FE. 
Telephone Oil 44 276 682756 



AN GBJBCT-CRIENIED ABCHTTEOTIJRE FCR SEN5CR DATA FUSION/TRACKING 
IN DENSE TABGSr ENVUXMlBns 



SPONSORS: ADVANCED SYSTH1 ARCHITECTURES, 

SDIO BM/C , Captain Johnson 

DEVEIX)PER: Advanced System Architecttires» Data Flision Department 
AmHOR: Edward J. G. Goodchild 



SIMIART 



The algorithm" is a target-oriented system desijtned to perform multi-sensor, 
multi-object data fusion and three-dimensional tracking in a comoutationaUy 
practical manner in dense target environments. It has been designed to ocerate 
with data from any mix of bearing-only and 3-D sensors, and tracks both 
individual targets and target clusters of arbitrary size. 

•nie algorithm is based upon an object-oriented system architecture, emheddintf 
the tracking filters and data association gatiwt functions within renlinabl« 
logical objects. Oie such object, or process, is assigned to track eanh 
perceived target, thus exploiting to the maximum the inherent parallelism in the 
data fusion task. A tracking process has the sole task of increasing its 
knowledge about the body it is modelling by examination of sensor data, and 
using any such reports that are relevant to that body. 
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!• CONTEXT 



!•! PKXESSING CHAINS 

The processing chain characterising the data fusion algorithm differs depending 
on the type of data provided by the sensor. The systan has been designed to 
accept data in the. form of hits or-scan-to-scan correlated tracks in either two 
or three diiDensicns. 

Track initiation from uncorrelated 2-D hits follows the processing chain shown 
in figure 1, which indicates how new track hypotheses are formed fron hits from 
all available sensors without an intermediate stage of track formation by 
individual sensors. Track initiation from sensor tracks or with three-dimension- 
al hits follow variants of the "Centralised Sensor by Sensor" processing chain 
shown in figure 2^ . Track maintenance follows variants of the "Centralised 
Sensor by Sensor*" processing chain for all types of input. The principed 
variations from figure 2 concern the type of gating applied to the sensor 
inputs; this is dependant on the sensor data type. 

The algorithm covers the 3-D tracking functions and the gating and association 
of incoming data with the 3-D track projections, but excludes any single-sensor, 
scan-to-scan correlation and track forming functions. 

1.2 PHASES OF TOE SDI BATTLE 

The algorithm has been designed for birth-to-death tracking of ballistic 
missiles, from launch all the way through to the terminal phase. Effort has been 
concentrated, thus far, on the transition from boost to mid-course phases, 
including the period of target proliferation during the post-boost phase. 

1.3 INPUTS 

The algorithm has been designed to accept any available sensor data with no 
prior assumptions on the number or configuration of sensors. Sensor data may be 
clustered or otherwise, and may be scan-to-scan correlated or otherwise. 

1.4 OUTPCTS 

The fidgorithm produces full three-dimensional tracks of individual targets, of 
clusters of targets, and of individual members of clusters where sufficient data 
has been received to resolve them. Extensions are planned, but not currently 
supported, to include target discrimination functions, and to provide extra- 
polation of the tracks both forward and backward to yield aim-points and launch 
parameters . 

2 NOTABLE FBAIIIRES 

- the ability to handle any mix of sensor data, as described above. 



i Figure 1 in attachment 4 of IDA Memorandum dated December 27*^ 
from Gabriel Frenkel. 
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the ability to initiate full 3-D tracks from three uncorrelated 2-D 
observaticxis on a target from any sensors (even a single sensor provided 
that its movement relative to the target is sufficient), and to prune 
swiftly the rapid proliferation of false track hypotheses. 

the scalability of the architecture to any conceivable threat size, with 
predictable growth with growth of the threat environment. 

3 SBNSCR ABC&nSCnBS & TBSBftT SCENARIO 

The Edgorithm, although originally concieved as a generic data fiision architect- 
ure, has been developed in the context of the SDI. The scenario driving the 
design has been that of a mass ballistic missile attack including many launchers 
producing vast nunnbers of mid-course bodies (>10« ) . These have been assuned to 
travel in cliasters of up to a few hundred objects each. 

The system has been designed to operate with any likely configuration and nunber 
of sensors. The algorithm has been designed to take account of measurem^t 
errors and biasses (including sensor own-position reporting, pointing accuracy* 
and bearing resolution) but performance measiirements have not been made* 

Three kinematic models have been assumed for targets, a pure ballistic model for 
the mid-course {iiase, a ballistic model with acceleration along the trajectory 
for the boost phase, and a ballistic model with variable acceleration, pre- 
doniinantly along the trajectory for the post-boost vehicle. This variatble 
acceleration is included as a noise term in the kinematic model. 

At present, tracking is on target position and motion only. Versions of the 
design exist which include photometric parameters in the target state vector • 
and it is planned to make use of these for target discrimination purposes in the 
future. 

No specific assusoptions have been made, to date, as to the nature. Quantity, or 
pre-filtering of background clutter. 

4 SENSOR IfXfEL/HOCBSSTHa 

In principed, the algorithm is not limited to operation with particular sensor 
types. The evolution of the design in the SDI context has directed concentration 
towards operaticm with passive, focal-plane, IR sensors. It will operate with 
both scanning and staring sensors , producing either geometric data only or both 
geometric and photcmetric data. Scan rates and frame rates may be fixed or 
variable. 

The algorithm requires data to be supplied on at least: observation time, sensor 
position, bearing of observed object in two orthogcmal planes. No other limits 
have been aiasumed on sensor signal processing capabilities. 

5 TRAGE miTIATIGN 

5.1 TARGET TRACKS 

Tracks are initiated directly in three-dimensional form even from un-correlated 
2-D hits. They are initialised by association of any sensor reports not 
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associated with existing tracks to form new-trac* hypotheses. These are tested 
for realism and likelihood at every stage, with hypotheses failing either test- 
being deleted* Association of three hits in a new-track hypothesis is suffici- 
ent, in general, to allow its promotion to a fiill track. 

There is no necessity to initialise the system with tracks handed over from 
elsewhere following a cold-start. Provision for svch hand-over could be added to 
the algorithm without difficulty, and could reduce the settling-time of the 
system under cold-start conditions. 

5.2 CLUSTER TTZACKS 

All launch v^iclea are assxmied to produce clusters of mid-course bodies, thus 
cluster tracks are generated for all boost-phase tracks. At this stage they have 
a membership of one target, i.e. the booster itself. New tracks formed by 
splitting during the post-boost phase are assigned as meinbers of the cluster. 

Cluster tracks are also formed later by association of a mmiber of non-cluster- 
member tracks all having near identical trajectory behaviour. 

6 TRACK MAXyrSNANCE 

6 . 1 DATA ASSOCIATION 

A multiple gating approach is used for data associaticxi. Firstly, new data is 
gated with the cluster tracks, using a simple rectangular gating strategy. A 
two-level gating strategy is then applied to the data by target tracka belonging 
to clusters with which the data is successfully associated. The first level 
comprises a rectangular gate; svccessful association at this level leads to a 
more stringent, ellipsoidal gate. Only data passing all stages succesfully is 
used for track updating. 

6.2 STATE ESTIMATION 

A Kalman filter technique is onployed for target track state estimation. 
Currently, this comprises a six-degree-of -freedom filter, tracking positicm and 
velocity. Plans exist to extend the state vector to include aim-point, launch 
point, and photometric parameters with appropriate filters. 

A similar six-axis Kalman filter is used to maintain the cluster track 
trajectory, and cluster extent. The cluster track also maintains its target- 
track meinbership list. 

6.3 TEIACK FBOMDTION/DEMQTION 

Tracks are continiKxisly monitored for validity. Pruning testa implemmted and 
planned include: consistent aiKl sensible track behaviour, behaviour of the error 
co-variance with time, and continued updating with new sensor data. Conaistait 
failure of any of these tests would cause deletion of the track. 

7 TRACK FILE MAINTENANCE 

Currently not implemented. 
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8 OUIFUr TO m/<P USERS 



Currently, the output of the algorithm is all the track state vectors, as they 
are updated* 

9 OGHEVEATIGKAL RBQUHZSIBm 

The algorithm has been designed to operate on a large, concurrent, message- 
passing processing system, with the ability to distribute the system over a 
number of platforms. The algorithm has been designed explicitly to maximise 
processing parallelism. The ntsnber of processors required in such a machine 
would probably exceed 1000, in a loosely-coupled architecture. 

The edgorithm has been developed with the Auto-G CASE tool, vising the G 
notation, with automatic code generation of Ada and the Ada-based SAIKT language 
developed by the IDA. Simulations have been run of jpart of the design coded in 
SAIMT on SUN 3 workstations. 

10 CURRBMT STATUS 

The algorithm is in a concept proving stage t with the core part of the design 
having been completed in Auto-G. A portion of it, the target tracking process 
has been converted to SADMT code. 

Only limited performance optimisation has only been carried out. 

11 FERFCBMANCB MEASURES AND RESULTS 

The target tracking process coded in SAD^f^ has been run successfully in 
simulations using the IDA SAD^f^ Simulation Framework, demonstrating track 
initiation and maintenance. Performance metrics have not been produced. 

12 REPGRTS 

ORIGINAL PRESENTATION OF CONCEPT:- "An Object Oriented Approach To Data 

Fusion"; Chapter 3, Section 3.5 of "Application Of Artificial Intelligence To 
Command And Control Systems"; C. J. Harris (Ed); Peregrinus Press for the 
Institution of Electrical Engineers. 

INITIAL DESIGN OF ALQORITTM:- ASA Technical Report T88/001, dated 15* ^ 

hfarch 1988, produced for the SDIO as the final reix>rt of the original DoD SDIO 
contract SDIO84-87-C-0040, entitled "Demonstration Of Specification Methodology 
For SDI Data Fusion". 

DETAILED DESCRIPTION OF ALQORITHM:- ASA Technical Report T88/009, dated 15* 
June 1988, produced for the SDIO as the First Phase reix>rt under the extension 
to DoD SDIO contract SDIO84-87-C-0040, entitled '*Data Fusion Architecture 
Refinement And Simulation". 

DEVELOPMENT OF SIMULATION AND PRESENTATION OF SIMULATION RESULTS:- ASA 
Technical keport T88/014, dated 15«^ December 1988, produced for the SDIO as the 
Final Report of the extension to DoD SDIO contract SDIO84-87'-C--0040. 
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AUTHOR'S COMPANY/ORGANIZATION: ALPHATECH, Inc. 
AUTHOR'S PHONE: (617) 273-3388 

ADDRESS OF AUTHOR: 1 1 1 Middlesex Turnpike, BurHngton, MA 01803 

TITLE OF ALGORITHM: Multiple Information Set Tracking Correlator (MISTC) 
SPONSOR: U.S. Army Strategic Defense Command CSSD-H-SBY 
DEVELOPER: ALPHATECH, Inc. 

1. DESIGN SCENARIO 

The algorithm design was driven by midcourse scenarios after RV deployment and 
before re-entry. Total threat size before daui partitioning was assumed of the order of 10^ 
objects, group size after data partitioning was assumed of the order of 100 objects. A range of 
target densities were used for simulated scenarios in evaluation of group tracking and 
correlation. The threat SDC-M was used for evaluation of data partitioning. All objects were 
assumed to fly Keplerian trajectories without maneuvering. Models for space-based, airborne, 
and ground-based sensors were used to allow representation of SSTS, GSTS, AOS, and 
MGBR. Handover track data could also be included. A variety of sensor architectures were 
used. Registration, calibration, etc. was assumed to be performed prior to tracking; only 
sensor noise errors were modeled. 

2. SENSOR MODEL/PROCESSING 

The sensor models included effects for detection (constant probability of detection 
within the sensor FOV, Poisson uniform false alarms in the FOV), measurement resolution 
(measurements within specified sensor resolution were combined — however, no algorithms 
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were developed to handle this during the MISTC project), measurement accuracy (additive 
Gaussian noise), and single sensor data processing (tracking in the single sensor's 
measurement coordinates — e.g., angle scan-scan correlation and tracking). Attitude, 
navigation, etc, errors and biases were not modeled. The sensor output to the tracking 
algorithm included time, CSO indicator (however, not used in algorithms), measurement track, 
measurement track error cdvariance, measurement track ID, measurement track probability 
distribution of target class, nieasurement track rime. The measurement track format varied 
depending on the type of sensor (e.g., four dimensional angle and angle rate state for angle- 
only sensors). 

3. TRACK INITIATION 

Tracks are initialized as 2-D tracks until multiple sensor correlations are made (then a 3- 
D track is initialized). Handover data can be used if available, but is not necessary. Tracks are 
processed individually and sequentially. 

4. DATA ASSOCIATION 

Data association was accomplished in two stages. The first stage partitioned one scan 
of LOS data from multiple sensors (two or three) into spatially separated groups. Then each 
group of data was processed by the tracking/correlation algorithm. Within the group eight 
different data association algorithms were investigated, including three zero-scan, pairwise 
approaches (row-column, row-column with backtracking, and optimal RELAX assignment) 
and five N-scan, multiple hypothesis approaches (multidimensional, maximum, marginal 
xetum, branch and bound, and three new algorithms by Tsaknakis). These algorithms are 
described in [2]. The same correlation algorithm performed scan- scan and sensor-sensor 
correlation. Hypotheses were scored by likelihood ratios. Ten branches were allowed per 
target (five appeared adequate). The initial data partitioning was able to keep the group sizes 
fairly small so that gating and pruning were sufficient to keep the number of hypotheses firom 
growing. p 



5 . STATE ESTIMATION 

State estimation was accomplished using extended Kalman filters assuming Keplerian 
motion between measurement times. The algorithms developed did not handle track spawning. 

6 . TRACK FILE MAINTENANCE 

The track file stored track time, track state type (e.g., 2-D or 3-D), mean, error 
covariance, ID, probability distribution of target class (e.g., RV, balloon, light replica), list of 
sensors correlated for that track and corresponding sensor track IDs. 

1 . OUTPUT TO BM/C3 AND USERS 

Not applicable. Entire track file was assumed to be available to user systems. 

8. PERFORMANCE MEASURES 

Performance measures included number of true tracks generated, number of false tracks 
generated, trae track accuracy, track time to initialize, and life-time. These are detailed in the 
report [1]. Performance bounds were not computed. 

9. COMPUTATIONAL REQUIREMENTS 

Computational requirements were obtained by partitioning realistic threat tapes to obtain 
a distribution for group size and then running the tracking/correlation algorithm on groups of 
different size to obtain empirical estimates of run-time. Run-times on a VAX 1 1/750 (using 
VMS FORTRAN) were translated to MEAS (million equivalent floating point additions per 
second) and scaled to larger threat sizes. The details and results of this computational 
evaluation are contained in [3]. The rough requirements were 90-120 MEAS per 10000 objects 
for average throughput with about 80 MBytes for total memory. 

10. NOTABLE FEATURES 

The data partitioning algorithm proved effective in decomposing the sensor data into 
groups small enough for uniprocessor computation. The Tsaknakis algorithms exhibited near- 
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optimal performance with processing not much greater than the simplest algorithm. 
Description of these algorithms and their tracking performance are contained in [2]. 



11. CURRENT STATUS 

Algorithms have been implemented for sequential, off-line processing in FORTRAN 
(VMS)- Code was not optimized Plans are to incorporate CSO tracking and resolution into 
algorithms and to implement on different parallel processors. 



12. REPORTS 

All reports are unclassified with the usual SDI limitation on distribution and are available from 
the U.S. Army Strategic Defense Command. 

1 . Catiin, R. A,, and R.B. Washburn, Evaluation Methodology for Multiple 
Information Set Tracking Correlator (MISTC), CDRL Item A009, ALPHATECH 
Technical Report TR-385, Burlington, Massachusetts, March 1988, 

2. Allen, T.G., L.B. Feinberg, R.O. LaMaire, K.R. Pattipati, H. Tsaknakis, R.B. 
Washburn, W. Wren, P. Patterson, and T. Dobbins, Mw/rip/e Information Set 
Tracking Correlator (MISTC) Final Report CDRL Item A007, ALPHATECH 
Technical Report TR-406, Burlington, Massachusetts, September 1988. 

3 . Washburn, R.B., Multiple Information Set Tracking Correlator (MISTC) 
Processing Requirements, CDRL Item A006, ALPHATECH Technical Report TR- 
348-2, Burlington, Massachusetts, September 1988. 
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ADDRESS OF SUEMTTIER: 1100 Wilson Blvd. Siaite 1710, Arlington, Va. 22209 
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TITLE OF AIX3QRITE1M: SDI Mldoourse Tracker/Corcelator Algorithm (TRC) 
SPONSOR: Nacval Research laboratory (Dr. Kurt Askin, Mr; Steve McsBumett) 
DEVEIORER: Ball Systems Engineering Division and Daniel H. Wagner, Associates 



The follcvring responses are offered in response to questions posed in the 
summary related to the NRL Tracker/Correlator Algorithm under develcpoent by 
Ball Systems Engineering Division (VERAC, Inoorpbrated) and Daniel H. Wagner, 
Associates. Since scsob responses could result in lengthy algorithm 
descriptions that exist in our Algorithm Design Document, a reference to a 
pertinent section or sections of the docunent is provided in lieu of a xtore 
definitive answer. A copy of the document is enclosed for your convenience. 

1. DESIQJ SCENARIO 

A. Ihe algorithm design is. essentially scenario independent except that it 
takes advantage of the fact that objects tend to cluster into spatially 
inseparable groups in order to reduce the conibinatoriad e3?)losion that results 
from the hi^ target densily and potentially high sensor observation rate. 

B. No assus^stions about the threat size and density were made in the design, 
however the algorithm's implementation on the SUN ocaatputer provides upper 
practical limits due to processor speed and memory avadlability. A future 
inplementatidn, currently in p r ogress, on a parallel processor will extend the 
utility of the algorithm to greater numbers of objects. 

C. Ihe algorithm assumes some p r^roc e ssing by the sensors to remove clutter. 
It does hcwever treat the case of false reports in its scoring formulae. 

D. Ihe eilgorithm design is entirely independent of the sensor configurations. 

E. No fixed values are assumed for sensor measurement e r r or s and biaises. 
Algorithm parameters may be adjxasted to permit use of a full rai^e of 
postulated sensor errors. 

F. Post Boost Vehicles (PB7s) are assumed to undergo random accelerations 
during MIRVing. All Other objects are assamied to be in ballistic K^lerian 
orbits thixu^ midccurse [See sections 3.4 and 3.5.2 and i^pendix A] « 

EMCDDSURE 1 
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G. We have selected attrlhutes vdiich are independent of aspect for association 
purposes. 

2. SENSOR MDDEL/FRDCE5SING 

H. Sensor configurations and parameters are established b/ the MRL systems 
engineer. Hie algorithm naJces no assunptions except those d e scribed in section 

I. 2. 

B. TbB TRC receives the following informtion from the sensors: 

tjjm, sensor position and velocity, azimuth and elevation and errors 
IR only: irradiances at three wavelengths, S/N 
Radar only: range, radar cross section 

3. TRACK INinATIOT 

A. Currently tracks are initiated by a boost-phase handover report a ssu m ed to 
provide the 6-state estimate of an c±>ject (presumably a PB7) plus an 18-element 
covariance matrix. Enhancements to the algorithm, presently i:inder design, will 
initiate tracks from 2-D and 3-D sensor r^rts. 

B. Sensor reports are batch processed as indicated in Section 3.2.2. 

4. DftTA ASSCXnftTION , 

A. Reports from each sensor scan are associated against omrent target space, 
i.e. the algorithm is a r^rt-to-tradc correlator. 

B. Ihe adgorithm takes a multi-hypothesis approach. 

C. Coiibinatoriad explosion is controlled in several ways: 

1. local processors and clusters [Section 3.3.1] 

2. within a cluster, r^rt-to-map associations must pass a retention 
test to be considered for inclusion in a scene (hypothesis) [Section 
3.5.5] 

3. only hypotheses which pass a retention test are passed on for future 
i2se. Others are pruned. [Section 3.8.4.4] 

D&E. We do. not have a fixed d^Jth hypothesis algorithm. Hypotheses are deleted 
if they are below (in score) an input percentage of the score of the best 
hypothesis. There is also a maximum number of hypotheses which can be saved. 

F. Hypotheses are scored usinj a Bayesian approach and taking into acocunt a 
prcbability of detection, probability of false alarm, false and new target 
densities, and- r^)ort-to-map association scores. [Section 3.8.4.2 and. 3.8.5] 

G. A heuristic search algorithm is \ased to formulate the hypotheses. The 
algorithm enplqys knowledge about target space and a key assumption that no two 
reports in a sensor scan were generated fay the same target. [Section 3.8.4.3] 
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H. Tracks are spawned as described in Sectican 3. 8. 4.3 figure 3,8-2. T&rget 
maps are merged as described in Section 3.11. 

5. STAIE ESTIMftTION. 

A. The approadi for state vector estimation for a track is desoribedxrin- 
Section 3.10 and Appendix D. 

Bie state vector is a 6-dimensional position/velocity vector plxis 
covariance matrix plx:is the target tenpeafiature estimate and its. associated 
error. 

C. Ihe algorithm xsses a Kahlmah filter with equations of motion based on the 
modified Euler method. Ihe dynamical model is described in Appendix A. 

D. For information en spawning, track scoring and pruning see Section 3.14 as 
well as previously referenced sections on report*to-map asisociaticn and. scene 
scoring. 

6. TE^CK FILE MAINIENAWCE 

As above, see Section 3.14. At each periodic snapshot, constructed track 
ccsrplex data are out^xit to the track file to v^date the state estimates, for all 
tracks. CTCs vdiich are not output have been deleted frcm the track file. In 
the NRL TSS, the track file data base will be maintained by a separate CSCI. 
Output crcs will include a flag which indicates if the object is a FB7. 

7. OUTHTT TO EM/C3 AND USERS 

A. Snapshots of the current state estimate of objects which exceed a 
confide n ce (CTC wei^t) threshold are output to the track data base. [Section 
3.15.3] 

B. Data are output periodically, (noninally every 12 seconds) . Ihe snapshot 
interval can be varied by changing an iiput parameter. 

8. F£kl^Oi<MANCE MEASURES 

A. The algorithm is measured in four areas: tracking, target count, 
correlation and execution time. Please refer to the Test Bed Design Document 
(Enclosed) for a detailed description of all measures of performance. 

B. We are not using maximum likelihood estimation and do not have a 
theoretical performance measure. 

9. OCMrorftnONAL rbcoiremewis 

A. Throughput and memory requirements -* none have been established 

B. 'Ihroui^:$)ut is measured as number of reports processed per second, where a 
report is a single observation of a target or unresolvable groap of targets by 
a sensor. 
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C« The algorithm is isplemented on a SUN oarputer in irUKTHAN. 

O. Results can be considered scadable on a global, but not a locaU. basis. 
That is, objects frcm each missile are considered and processed independently 
of objects frcm other missiles. Within the o uuip u nent cl\:isters of a local 
processor, the problem scales non-linearly with the density of the targets. 
Problem should scale linearly with the number of missiles provdded eadi missile 
produces similar target densities. 

E« ThB results to date are enpirical based on sinple target scenarios. 
F. No ccaoiputational performance bounds have been isposed at this tmB. 

10. NOIABIE FEMURES 

A. The combinatorial explosion resulting from high target densities and report 
rates is limited by clustering and use of retention thresholds. 

B. Constructed Track Coirplexes give a useful output to other battle manager 
functions. These are a set of maps judged to be edtemate r^)resentations 
(under different report-to-track correlation hypotheses) of the sane target. 

C. Hierarchical algorithm architecture is designed for eventual parallel 
processing. 

D. Target non-kinematic attributes are processed and associated. Currently 
the algorithm uses tenperature derived from irradiance measureznents. 

11. CURRENT STRTDS 

A. A second generation algorithm has been conpleted and is in the testing 
phase. It has been installed at NRL and lANL. 

B. Performance optimization is currently being investigated. 

C. Hardware has not been qptimized for this algorithm. 

D. Future plans include: 

o esdiaustive testing and behavior analysis under stressing scenarios 
o continued enhancements underway or pr opo s ed: 

- track initiation without handover 

- inproved PB7 motion models 

- vise of color indices as association attributes 

- c^namic thresholding 

- dynamic new and false target density VEdues for scene scoring 

- inproved report-to-map assignment scheme for accelerated processing 
(monotonic logical grid) 

- - improved scene processor 
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- i api 'o v Bd scene processor 

- parallel iiiplemerttation on Butterfly 

12. KKRJKIS 

Relevant algorithm description documents are enclosed. Performa n ce testing 
reports are pending ocaipletion of testing and will be available frcm Dr. Hart 
Askin at MRL or Mr. Larzy Filippelli at Ball Systems Engineering Division. 
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(Complete as appropriate.) 



NOTES ft INSTRUCTIONS 

(The Siunmary/Abstract of the algorithm should be limited to a half 
of a page* The basic answers to the questions below should be 
limited to 3 pages. Additional information, referenced to the 
section numbers, whould be included in a separate appendix that 
should be limited to four pages. Classified information should be 
included in a separate supplement. The total information including 
Abstract/Summary, basic answeres, appendix and classified supplement 
shoulld not exceed eight pages. If some of the algorithms or 
details are proprietary, indicate what is proprietary and discuss 
only the non-proprietary aspects.) 

ABSTRACT luu APV 

(Describe, in general terms, how the algorithm works and if 
applicable provide a high level flow diagram.) 

1. CONTEXT 

J[ yfl^sit processing chain characterize your algorithm? (See 

Figs. 1-4.) What functions within the chain are covered by your 
algorithm, e.g.. Track Initiation, Track Maintainence , etc? To what 
phase(s) of SDI is your algorithm applicable - Boost, Post-Boost, 
Uidcourse, Terminal? What are your inputs, e.g. single sensor or 
multiple sensors, clustered data or not clustered data, etc? What 
are your outputs, e.g. 2-d tracks or 3-d tracks, launch parameters? 

2. NOTABLE FEATURES 

(List any features that distinguish the algorithm.) 

3. SENSOR ARCHITECTURE ft THREAT SCENARIO 

( Characterize the scenario that drives the algorithm design. What 
is the threat size an^ density? Whar are the assumptions about 
background clutter? Is all the clutter removed prior to tracking? 
What is the constellation size and orbit? What are the measurement 
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SimuIalioii88 Trackuig Mudt:! 

Thomas D. Cottschalk 

California lostilule Of Teciiaology 
356^8 

Pasadeoa, CA 91125, 818-356-6671 
ABSTRACT 

This note describes ttie tracking module for SlmulatioaSS - a boost/post-boost 
simulation of SDI functions for an SS18/ASAT direai scenario performed on the Mar- 
klU bypexcube. The simulation involves a number of separate tasks (enviromnent gen- 
eration, s e ns in g, tracking and battle planning) each running on separate subcubes of 
the Maxkm bypercube. WiUiin a given subcube, the particular SDI task is done in a 
concurrent fashion. Communications among subcubes are done asynchronously. Tlie 
Sim88 tracking task involves separate subcubes for sensing and tracking, wilh due 
tracking task of an individual subcubc/tracker in turn divided into 20 mono tracking 
(using only data from a single sensor) and 3D precision tracking (using 2D mono 
tracks from itself and an addidonal tracker/subcube). While the details of the present 
tracker are rather tlghdy ded to specifics of the Slm88 threat model, the overall con- 
current tracking prescripuon can be generalized to deal with aibitraxy threat scenarios, 
and modifications along these lines have been begim. 
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Couiext 

The traddog model for SimulatiooSS involves eight separate subcubes of the MoxiHI hypercube, 
perfonuiug MEO and GEO seosiag and tiacking ta&ks for boost aaJ post-boost phases oa uiulU*iaigct. 
sceoanos. lu barest terms, the task of the tracker is to process pixels from tlie semors aoJ ulilfuately 
produce trajectory parameter estimates for iodividual targets which are then passed oa to tl^ baiUe-. 
plamiiag phase of the iutegrated SimuIatiooSS demoostnitioa The overall processiog chaia is iUusiraie J 
io Fig.(l)« Al the GEO level, each seasor develops pixel sets at iadividual scans and passes them to aa. 
asso ci at ed tracking function. The tracken develop 2D kinematic mono tracks from these data. Mi>no 
track files ate exchanged between stereo partners and associated to fonu 3D kiueuiatic tracks (EQ posi* 
tion and velocity estimates) which are then used to solve for launch parameters according to a powered 
flight model Tracking at the MEO level is concepmully identical to thai of the GEO system, with the 
addition of direct initiation of precision tracks from the GEO tracker. 

Figure (2) shows the general processing strategy used for both the GEO and MEO trackers during 
a single scan.. The 20 tracking is done using simple, 2D kinematic Kalman iillers for each of the two 
projections measured by the sensor. Tracks are initiated using a 3-scan batch processor and extended 
using a track-splitting formalism. Focal plane tracks of sufficient age (typically 4 or 5 total scans) are 
propagated to a common reference time and collected imo a 'Tocal Plane Report'* which is exchanged 
with the stereo partner. The association of 2D reports into 3D uracks is done in terms of projected posi- 
tions and velocities along a common reference line (the *iuiersecUon* of the two focal planes). The pre- 
cision tracking task actually consists of two distina paru. For precision track initiation, the EQ posi- 
tion and velocity vectors from \bc stereo association are used to determine initial values for tiie launch 
|)arameters. Once initialized, parameters are updated on subsequent scans by means of extended Kalnian 
filters. 

The focus of the SimuiationSS effoit was the establisbmem of a flexible, integrated fnunework for 
boost/post-boost SDI functions, and the demonstration that the computations for tliese tasks could be 
done concunentiy in an effideot manner. In ihc interest of meeting somewhat ambitious dcmoostration 
deadlines for the entire simulation, a number of 'shortcuts' la \bc tracking formalism have been taken - 
most notably, the decision to design a tracker which is intimately tied to a specific powered flight 
model and the use of extremely simplistic sensor models. Present work is aimed ut removing these 
defects, witii the intent of producing a generalization of die SimuiationSS tracker wluch can process 
threat tapes. 

Architecture And Threat Scenario 

The generic sensor architecture for SimuiationSS consists of a ring of LWIR sensors in geosta- 
tionary oibits and two rings of SWIR sensors in circular, polar orbits. The number and nature of sensor 
rings as well as the number of sensors in an individual ring are specified at run time through an archi- 
tecture definition file. 
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While the eaviromhent generaiion phase of SimSS aciuiJly 'llie:^* ih^ enure seming cao^itcllauoi], 
only subsets (typically 2 GEO sensor/tracker paiis and 2 MEO paiis) afe siiiiuiatcd ia high fidelity as 
separate subcubes within the hypercube. The choice as to which particular sensor/truckers are actually 
simulated is detennined at run time-within the overall simulation Iramcwork. 

The threat scenano for SimSS consists of a disuibuted launch of 'generic* SS18*s togedier wiUi 
some number of ASATs. Trajectories for individual SS18*s and ASATs are described usuig a simple, 
analytic powered flight model (TowerSS*) which was developed spcciucally for the SiiuSS taJc An 
individual launch is described by the longitude and ladtude of the launch site, the launch time, the iui* 
tial launch azimudi and a final thrust angle for each stage. The SSI 8 model has two stages (6 panuue- 
ters) while the ASAT model is single stage (5 parameters). The SS18 thieat typically has about 40 
boostera from each of six launch sites launched wiUim a two minute wave. The exact nature of the 
threat (including targeting to specific sites) is run time recoofigunble. The ASAT dueat is generated 
dynamically at die start of the simulation, with the actual ASAT launches detemiined by the nature and 
kinemadc properties of the CV architecture. 

' Sensor Modeling 

« 

Sensors for die current program are embarrassingly siniplisdc : 3D target posidons are projected 
onto sensor focal planes giving a 2D data point. For the GEO sensors, (be focal plane is divided into a 
number of rectangular pixels, and die reported observation is simply iLe center of an active pixel. The 
GEO sensor model does not have muld-hil disciiuunation within iDdividuuI pIxeL. For die MEO sens- 
ing model, the truth posiuon on the focal plane is simply smeared u^ung Gaussian distribudons. At 
present, Pq and diere is no noise. In typical runs, the scan time for die GEO system is S seconds 
and that for the MEO sensor is 10 seconds. 

In the integrated SunSS framework of Fig.(l), an endre subcube (typically 1 node) is used for 
each sensor. This is sufficient computing power to do a more realisdc description of die sensing pro« 
cess, and this task will be included at some level m die next genen&don of the simulaUou. Tlie actual 
fidelity of sensor simuladon remains to be deteniuned. 

Track IidtL*iioa 

Two types of track initiation occur within SimSS: 2D focal plane and 3D precision uacks. The 
inidadon algorithms are quite diffcrem. 

Fbr focal plane tracking, iniUation is done using a 3-scan batch processor. Al any time, die full 
sensor pixel lists for the present and two previous scans are maintained in memory, and the batch uiitia- 
tor constandy searches for provisional new tracks subject to essentially only three constmints: 

1) The three data used in the track must not coincide widi the last duee data in any exisung 
track. 
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2) The last datum ia the tnck must DOt coiacide wtdi die la^t datum to any esiublLiUed uuok 
(defined by Nsokus ^ ^o> typical cutofi Nq s 7. 

3) The 3-scan segmecu must be suffidendy straiglit, as deliued by a uuuiber of suuple. Iieuris- 
tic cuts 00 the appaiem acceleratioo. - _ _ "f 

For \bc 3D tracker* a precisioa track ioidatioo is allempled ^beoever tije ECI posilioa vector dcicr* 
mined by Ibe kinemalic stereo algorithm caiux>t be associated widi any existing preci^ou track. In such 
cases, the ECI posidon and velocity vectors of the target (together with covariai^e estimates) are first 
determii^d by a straightforward geometric aigohihmt and these values are dien used to detemune inidal 
launch parameter and covariance values (witliin the PowerSS framework) by a Newton-Rapbsoo iuver* 
sion of the equations of motion. 

The coiieladon of focal plane reports to form 3D EQ posidon and velocity vectors is a central 
pan of the precision inidation process. As noted previously, the focal plane tracks from the two mono 
trackers are first propagated to common reference times, atkl values of die focal plane posidons arul 
veloddes projected onto a common reference axis are computed. The two focal plane reports are then 
associated by forming associations between these projecdons according to a modified nearest-neighbor 
algorithm. 

Truck Maiulenance 

For the focal plane trackers, track extensions are done using a track-splitting formalism. 2D tracks 
are extended to predicted posidons using the system/measurement models and association regions based 
on the posidon covariance are formed. Each observation widun the association gate results in a separate 
extension of the iruUal track. Tracks are presendy deleted on a single missed detecdon (this can/will be 
generalized to accormnodate Pp ^ 1). Track pruning is done by a simple 'common history* algoridun: 
tracks are deemed equivalem if tl>ey use the same sensor observauons over the past N£QUiY scans, widi 
typical value Nsg^iy s 4. Since tiie kinemadc filler is high gain, die choice of whidi lra<:k is kept of an 
equivalent pair is inconsequentiaL 

For precision tracks, the associadon of tracks widi data is done using a Nearest-Neighbor global 
assodadon scheme. Specifically, predicted data positions are evaluated for each sensor, aikl a global 
associadon of predicted and actual data sets is done using a modified nearest-neighbor scheme. The pre* 
cision track is then updated using an Extended Kalman Filter, with each projection processed separately. 

Present Status and Future Geiteralizatiuus 

The tracking model based on the prevision descriptions has been completed and integrated widiiu 
the entire SimidadonSS frameworlc, and was pan of the successful fiiil SimSS demonstradon for ESD in 
January 1989. While the achievement of this milestone was a non-trivial accomplishment in itself, there 
are a number of shortcomings in the present tracker. Specifically 
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1) The present tracker is tied loo iuUiuately to tlie assumed PowexSS ibieut model. 

2) The tracker does not do as good a job as it sboidd in estimatiug target velocities duxiog 
post-boost maneuvers. 

Both problems have a coimnoa origia : the design of the precision tracker is based entirety on 
parameterized trajectories, with all updates of existing precision tracks done using Extended Kalman 
Filters for assumed trajectory models. For boost phase, tliere is no reason to believe that the PowerSS 
model might describe a get^ric threat During post-boost, the 'noisy*Keplerlan' system model has been 
found to be inadequate* 

The next generation tracker will attempt to 'solve' both these problems by using a 3D ECI 
kinematic sute vector formalism as tiie basic part of the precision track file (for purposes of track-hit 
assodadons), with parameter estimations according to arbitrary flight models entering in a maimer 
which is well-removed from the multi-target tracking logic. The essential open question for the new 
tracking scheme concerns the output of the tracker and the imeipreiation of this output by the Battle 
Plaxming module of Sim88. The present tracker produces trajectory parameters which are ti^o used by 
the Battle Planner to predict positions and select engagements. The assumed 'complete knowledge' 
predict ahead used in Sim88 will r^ed to be replaced by something more fuzzy/realistic. 
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Sim88 Tracking Formalism 
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TRACKING ALGORITHM FOR PROJECT SWAT 

Status: August 1988 
Submitted by Jack Liu 

Project SWAT is sponsored by DARPA with funding from SDIO. Phase 1 of this program 
is an 18 month effort to evaluate the applicability of applying DARPA's MOSAIC architecture, 
developed by ESL in supporting the development and implementation of algorithms required for 
tracking and battle management for SDI. To demonstrate feasibility, the phase 1 project scope is 
limited to the problem of midcourse track maintenance with the possibility of evolving CSO 
resolvability and low velocity object deployment One potential mechanization was developed 
based on the classical concepts of gating, track splitting, multiple hypothesis track spawning and 
pruning, and extended Kalman filter prediction and correction. To date (August 1988), scenarios 
consisting of hundreds of tracks have been successfully run with the algorithm. Single platform 
tracking has been demonstrated; multiple platform scenarios will be addressed in the near future. 
No claim is made regarding the dptimality of the algorithm or its computational efficiency. This 
algorithm has been selected on the basis of being representative of the class of algorithms that 
would need to be considered for development 

The key feature of the MOSAIC architecture is its ability to link heterogeneous (different) 
processors to address a common problem. MOSAIC provides the capability for extremely high 
data rates (80 Mbytes/sec in, 80 Mbytes/sec out) between each processor. This is an important 
feature because of the large number of objects that must be tracked. Sixteen processors can be 
interconnected in a tightiy coupled configuratioiilising MOS AIC. The processors can be parallel 
machines, algorithmically specialized processors, or general purpose computers and workstations. 
MOSAIC provides the capability to utilize the power of parallel machines and algorithmically 
specialized processors more effectively. MOSAIC has a programmable scheduler to control the 
data flows between processcn: nodes. The schedules are precompiled 

The SWAT algorithm for midcourse tracking selected for the project utilizes classical 
concepts. There is a data source that generates the scenarios consisting of the objects to be 
tracked, the sensor platforms, uncompensated stars, and system error sources. The data/track 
associator allocates the sensed data to existing tracks and spawns and prunes tracks when 
appropriate. The extended Kalman filter updates tracks and generates gates for use by the 
data/track associator. The track file maintains the databases of true object tracks and algorithm 
generated tracks, emulates the process of field of view extraction, and produces indicators of 
algorithm performance. These components are being allocated to separate computen connected to 
MOSAIC. The data source is being implemented on a Sun 4; the data/track associator, on a 
Convex; the trajectoiy estimator, on a Warp; and the track file, on an Encore. These computers 
were GFE for the project This configuration was chosen for convenient attachment to MOSAIC; it 
is not necessarily an optimal configuration for an algorithm development test bed. 

The data source consists of a scenario generator, a propagation model, and a sensor model. 
The scenario gen^ator specifies the objects and their flight parameters and the sensor platform 
configuration.. The propagation model provides selectable dynamics propagation in a non- 
spherical earth gravity field of up to J6 . TTie sensor model has the capability to implement an ''all- 
seeing" sensor, a staring sensor, and a coolie hat sensor. The data can be provided to the tracking 
algorithm in tenns of either sensor focal plane azimuth/elevation coordinates or platform centered 
inertial right ascension/decUnation coordinates. 

The data/track association algorithm selected utilizes a multiple hypothesis track spawning 
and pruning approach for midcourse track maintenance. Therefore, it is a multiple scan approach, 
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the scan level being variable and affected by a redundancy parameter and the density of objects. 
Each track gate is detennined by the predicted track uncertainty projected on the sensor focal plane. 
The closest detection within the gate is processed for updating the track. Other detections within the 
gate are used to spawn new tracks. Detections within overlapping gates are utilized by all affected 
tracks. Rectangular gates are used Tracks arc scored based on the log likelihood function of the 
detection relative to the prediction plus a maintenance bias. The scores are cumulative. Tracks that 
have no detections are degraded. When a set track redundancy level is reached, the low score ~ 
tracks are deleted Tracks that are very close in terms of estimated object state and uncertainty are 
merged probabilistically based on their track scores. 

The track estimator is a standard extended Kalnian filter prediaor/corrector. The state 
dimension is 6 (3 positions and 3 velocities). Predicted state estimates for determing the gate 
position are generated by Runge-Kutta solution of the differential equations of motion. The earth 
gravity model is selectable up to J6- The gate size is determined by propagatiung the state estimation 
errors using the filter dynamics partials, projecting this on the sensor plane, and adding filter 
assumed sensor jitter. 

The track file emulates two functions. The first is to perform trackfile database management 
for the tracking algorithm, and the second is to evaluate tracking perforaiancc by relating the 
algorithni generated tracks to the objects produced by the scenario generator. Algorithm generated 
tracks are kept in a geometrically organized database to enhance access of tracks in the sensor's 
field of view. Various tracking algorithm evaluation functions are being developed with a 
preliminiary set currentiy in place. 
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Abstract /Summary 

These algorithms provide a means to track multiple clusters of 
objects and take advantage of data from multiple^ sensors. 
Multiple sensor data is especially useful with passive sensors 
because clusters can appear to cross or overlap in the field of 
view of. a single sensor. 

The cluster tracks can be used to initiate tracking of; the 
individual objects when the closely spaced objects are resolved. 
These algorithms use only a fraction of the computing resources 
normally required to initiate and maintain tracks on individual 
objects after deployment. 

These cluster tracking algorithms were conceived to overcome the 
difficulties and meet the unique needs posed by the SDI target 
threat. Cluster tracking is challenging because a cluster, as- 
seen from different platforms will differ significantly. The 
apparent size , and shape of a cluster varies from sensor to 
sensor. The apparent number of objects in a cluster can also be 
different for each sensor. 

These algorithms are based on new concepts that permit the 
estimation of not only the cluster centroid position and velocity 
in inertial space but also the cluster "extent." The extent 
establishes the cluster size and shape in inertial space, not- 
just relative to a particular sensor. 

A major feature of these algorithms is that data from multiple 
sensors is combined in a natural way. As a consequence, the 
approach is quite general. It is applicable to passive or active 
sensors or a combination of both. This approach is also 
applicable to tracking extended objects, such as plumes, and may 
be useful in other applications such as tracking a salvo launch. 
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Multiple Sensor Cluster Tracking Algorithms 



1. CONTEXT 

These algorithms were designed to provide a means to track 
clusters during the early midcourse stage of SDI, starting with 
deployment of clusters. In this context^ a cluster might consist 
of a reentry vehicle and attendant objects, including decoys.* 
Immediately after deployment of a cluster, many objects will not 
be resolved because of sensor resolution limitations. Thus 
initially, the cluster might appear from the sensor data to be a 
single "clump" (unresolved closely spaced objects) or a small 
group, namely, a collection of clumps and individual targets. 

It would be difficult and too demanding of a processor to track 
each individual target or clump during the early midcourse phase 
of SDI after deployment of a cluster. After a cluster is 
deployed, many new objects will continue to appear as more and 
more objects become resolved. As a consequence, the number, type 
and size of resolved objects of a cluster can change rapidly 
during this early phase. To compound the problem, the number and 
relative location of objects in a cluster would vary from sensor 
to sensor. In addition, the targets are initially so close that 
the high density can make individual target tracking impractical. 
See reference 1 or 2 for a discussion of these difficulties. 

These cluster tracking algorithms permit efficient tracking of 
cluster position, size and shape in spite of crossing and 
overlapping clusters and the other difficulties outlined. This 
permits a smooth transition from post boost vehicle (PBV) to 
individual target tracking without overloading the processors. 
The algorithms are intended to be compatible with birth-to-death 
tracking, which has been conceived as a practical solution for 
the computationally intensive SDI midcourse tracking task (see 
reference 1 or 2) . However, these algorithms should be useful 
with minor modification for other applications. 

While this approach could be used for independent sensor tracking 
as in the Type I processing chain, the best use is for multiple 
sensor tracking. This approach can be used effectively in any of 
the multiple sensor processing chains, namely, Types II, III and 
IV. The four types of processing chains are defined in the 
appendix. 

This approach should be most useful in the early midcourse phase 
of SDI, before almost all the targets are resolved. In that SDI 
phase this approach would be useful in three ways: 1) to provide 
cluster handoff information from one sensor system to another, 2) 
to combine cluster data from different systems, and 3) to combine 
data from multiple sensors within a specific sensor system. Thus 
SSTS, GSTS, GBR and AOS are applicable sensor systems. 
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Multiple Sensor Cluster Tracking Algorithms 



The inputs to the tracking algorithms are assumed to be the 
observations as provided by the sensor signal processor. Either 
or both active and passive sensor signal processing data can be 
used effectively. For a passive sensor^ the tracker could use 
simply angle location data (azimuth and elevation) for each 
detected signal (observation) without any' indication of whether 
it is a resolved object , extended object or unresolved closely 
space object. 

On the other handr additional information such as object 
intensity, multiple wavelength intensities , and object type and 
size could be used effectively and should improve tracking 
performance. The outputs (during most of the midcourse cluster 
tracking phase) include the estimated position (centroid 
location) , velocity and extent of the cluster (group) in inertial 
coordinates. 

Certain aspects, phases and outputs of this approach are Hughes 
proprietary so some details have been omitted from this Tracking 
Algorithm Survey. Also, this is not the only approach being 
developed at Hughes. Because of the importance of cluster 
tracking to SDI, other alternative cluster tracking methods are 
also being explored. 



2. NOTABLE FEATURES 

The critical issue that this approach addresses is the tracking 
of multiple clusters with multiple sensors. The use of multiple 
sensor data greatly improves the tracking accuracy and the 
information that describes or characterizes a cluster. 

This cluster tracking is designed to provide an effective and 
efficient transition from cluster deployment to individual target 
tracking. Through cluster tracking, the onboard processor size, 
weight and power can be substantially reduced. 

These algorithms accommodate the differences in the data obtained 
for a cluster by sensors at different locations. This approach 
is quite general and can be adapted to various types of 
applications and sensors. It is designed to process multiple 
clusters that overlap or cross in the field of view of passive 
sensors. It has been modified and tested for tracking extended 
objects, i.e. small targets larger than a point source. 

The way the extent is characterized, parameterized and processed 
in this approach is believed to be unique and novel. The extent 
estimate as formulated greatly facilitates both the processing 
and the usefulness of the output. 



B-33 



Multiple Sensor Cluster Tracking Algorithms 



3. SENSOR ARCHITECTURE AND THREAT SCENARIO 

The initial design of the algorithms was for SDI midcourse 
tracking. For initial feasibility testing, the simulation 
modeled an above the horizon infrared sensor system akin to SSTS 
with two platforms. The platform altitude was in the range of 
1000 to 5000 km, the sensor resolution was in the range of .001 
to .005 degrees and the angle measurement error standard 
deviation was from one-fourth to one-half the resolution. 

A typical cluster in the initial simulations was 40 targets (RV's 
and decoys) with a maximum dispersion rate of 2 meters/sec, 
nominally. The targets trajectories were generated using a 
ballistic model including the J2 term. See reference 4 for 
further details. Further testing with a range of simulation 
parameters in underway with limited funding. 



4. SENSOR MODEL/SIGNAL PROCESSING 

This approach can accommodate false and missing observations. 
For simplicity in initial feasibility testing, clutter and false 
signals were not included. A simplified signal processing 
simulation was used that included missing observations -due to 
unresolved closely spaced objects caused by finite —sensor 
resolution. Gaussian errors were first added to the true 
position for each target in the field of view for a sensor. The 
resulting clumps and resolved objects were passed to the track 
processor. The worst case was modeled by not providing 
observation intensity information to the track processor, which 
received only azimuth and elevation angles and time for an 
observation. The range of frame times was from five to 100 
seconds. 



5. TRACK INITIATION 

For initial testing the cluster tracks were initiated by track 
spawning based on the estimated PBV track. The centroid 
covariance matrix was larger than that of the booster state and 
the time of cluster deployment was assumed known. The covariance 
matrix for the extent was set to the identity matrix multiplied 
by a large scalar value. These initial conditions were then used 
for track maintenance processing. A more advanced track 
initiation approach has been designed but requires further 
development . 
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Multiple Sensor Cluster Tracking Algorithms 



6« TBACK MAJENTEMANCE 

In this approach^ the cluster is modeled in inertial space as an 
ellipsoid. This is modeled in two parts ^ the centroid and the 
extent. The., elements of- the centroid. state vector . are the 
inertial position and velocity. In this approach the extent is 
the second central moment in inertial space of the objects in the 
group. 

The extent provides information on the size and shape of the 
cluster relative to the centroid. The extent not only describes 
the characteristics of the cluster^ but it has a very important 
function in the track processing. The extent can also play a 
vital role when making a transition from cluster to individual 
target tracking. 



6.1 TRACK MAINTENANCE - Data Association 

The extent is used to determine which observation (signal/object 
detected by a sensor) belongs to which group. This is a vital 
function for tracking multiple clusters^ especially for crossing 
or overlapping clusters. Based on the projection jDf. the 
predicted extent on to the field of view of a sensor^ ^ gate 
(validation window) is computed for a cluster. 

For an isolated cluster, the first and second central moments in 
2-dimensions of the observations in the gate are computed. For a 
Type I, II or IV processing chain, these moments are passed to 
the filter described in Section 6.2 to update the established 
tracks. The four types of processing chains are defined in the 
appendix. For a Type III processing chain, the moments computed 
from multiple sensors must first be combined and then the 
resulting estimated moments in inertial space are passed to the 
filter. 

For clusters that are overlapping or crossing in the field of 
view of a sensor, association processing is needed to resolve 
the ambiguities. Association algorithms have been designed but 
require further development. 



6.2 TRACK MAINTENANCE - State Estimation (Filters^ 

This approach employs two filters, one for the state of the 
centroid and another for the extent. The filter for the 6-D 
centroid state is a simplified extended Kalman filter. The 
extent state has been simplified to 6 elements, which greatly 
reduces the processing load. The filter for the extent is : a 
pseudo-linear filter. The details and unique approaches used in 
this filter are provided in reference 4. 
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Multiple Sensor Cluster Tracking Algorithms 



7. TRJICK FILE M21ZNTENANCE 

The track files in the initial tests contained the two states and 
the corresponding covariance matrices. The advanced versions of 
this approach retain some additional information for purposes of 
association* 



8. OUTPUT TO USERS 

During early cluster tracking, the available data is the number 
of clusters, the number of objects in each cluster^ the state 
estimates (centroid and extent) and corresponding covariance 
matrices. After many of the objects in a cluster are resolved, 
it is practical to transition from cluster to individual target 
tracking. At that time the available information includes the 
state vector and corresponding covariance matrix for each target 
and additional information as needed, such as predictions and a 
target object map. 



9. COMPUTATIONAL REQUIREMENTS 

With isolated clusters the reduction in throughput and memory is 
approximately a factor of 2/N, where N is the average number of 
targets per cluster. For example, with an average of only 20 
targets per cluster the savings in track processing would be 
roughly a factor of 10. With overlapping clusters the savings is 
not as much and would depend on the processing accuracy required. 
Further detailed analysis is required to determine the processing 
load for a specific application. 



10. CURRENT STATUS 

Feasibility tests have been conducted successfully for cluster 
tracking under realistic conditions using practical values for 
sensor resolution and measurement errors. Earlier testing 
results with finite sensor resolution but without measurement 
errors are reported in reference 4. Further limited development 
and feasibility testing is underway on Hughes IR&D. 

While single sensor cluster tracking is relatively mature (see 
reference 3), multiple sensor cluster tracking has only recently 
been addressed* Accordingly, new ground is being broken, even in 
how to evaluate performance as well as how to track clusters. 
Substantial algorithm development and testing remains to evaluate 
the various phases and aspects of the approach and then adapt- it 
to a specific sensor application. 
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11. PERFORMANCE MEASURES 

Measures of performance for the cluster centroid can be similar 
to the usual measures for an individual target. However as 
formulated, the extent is new and therefore new meaningful- 
measures of performance must be established. One identified 
measure of extent estimation performance is the ratio of the 
estimated to actual volume of a cluster. Another measure is the 
eigenvalues of the estimated relative to the actual cluster 
second central moment. These measures and some early test 
results are discussed in reference 4. 
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13. APPENDIX - Types of Processing Chains 

Four types of processing chains have been adopted by the E)A SDI Tracking Panels for ease of 
reference. Reference 2 discusses these generic processing chains and summarizes the four types 
as follows. 

Type I: Independent Sensor Processing - Tracks are processed for each sensor independently of 
the data from the other sensors. 

Type n: Hierarchical Processing - Sensor level processing is followed by track fusion. 
Frame-to-frame association and filtering are followed by sensor-to-sensor processing. 

Type ni: Observation Fusion - Multiple sensor, observation association and target position 
estimation is followed by firame-to-frame association and filtering. Sensor-to-sensor processing 
precedes frame-to-framc processing. 

Type IV: Centralized Pnx:essing - Observation-to-track association is followed by multiple 
sensor filtering. Association and filtering is performed on each frame of data from each sensor 
as it becomes available. 
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ABSTRACT/SUMMARY 

The algoritha processes measurements from a constellation of passive 
LWIR sensors to generate precision cartesian state vectors on 
midcourse objects, clumps or groups. The algorithm performs birth 
to death tracking on RVs , penaids and PBVs . The algorithm also 
performs group tc object tracking. 

1 , CONTEXT 

Track Initiation is performed using the Mono- theu_ Stereo 
architecture. After Stereo tracks are initiated, track maintainence 
is carried out using the Centralized architecture. The algorithm is 
applicable to Midcourse. The inputs are line of sight measurements 
from multiple LWIR sensors. The line of sights may reference 
distinct objects or clumps of objects (CSOs) . 

2 > NOTABLE FEATURES 

The algorithm uses a Bayesian network architecture to carry out 
state estimation and association. The algorithm decomposes the 
association probl<»m into scenes. Each scene is a distinct 
association problem, either a pontact-to-track or track-to-track 
problem. A variant of the A* sfearch algorithm is used to generate 
feasible scene hypotheses. The algorithm performs birth to death 
and group to object tracking. 

3, SENSOR ARCHITECTURE ft THREAT SCENARIO 

The threat scenario driving the algorithm design is the TSGBl 
threat. Background clutter is lassumed to consist of stars, RSOs , 
structured background and nuplear redout. All clutter is uot 
assumed to be removed prior to tiracking. The constellation consists 
of 18 satellites inclined to 90 degrees. The RVs and penaids are 
assumed to be in free fall. Tfhe PBVs execute typical laydown 
trajectories. 
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4, SENSOR MODE /PROCESSING 



The sensor is a scanning, 3 color LWTR sensor. The frame time 
varies from 5 sec to 30 sec. The scan 5ipeed can also vary. CSOs 
arc lEodeiied usir.g a :i:nctional model. The interl'acc to th*- tracker 
consists* of [tine, az, el, snr, covariance, extended object 
indicator, extent parauneters] . 

5, TRACK INITIATION 

Tracks are initiated as 2-d tracks. After 4 to 6 updates a passive 
ranging algorithn is executed to initiate a coarse cartesian track. 
Cold start initiation and warm start intitiation are performed* 
Tracks are processed both individually and in a batch. False alarm 
and new track densities as well as chi-square scores are used to 
score new tracks. Tracks satisfying the firm track criteria are 
promoted; otherwise they are deleted. 

From a data association standpoint, the track initiation and track 
maintainence phases are not clearly separable since a variant of a 
multiple hypowhese algorithm is used. 

6 . 1 TRACK MAINTAINENCE - DATA ASSOCIATION 

A track is considered to have completed initiation when a precision 
cartesian state vector is created. This requires that individual 
mono tracks be formed. Passive ranging creates coarse cartesian 
tracks and track-to-track association forms precision cartesian 
tracks completing initiation. — 

Pruning, merging and clustering are used to control the combinatoric 
explosion. In addition, scenes are managed to assure that they do 
not get too large. The hypotheses are scored using a Bayesian 
approach. A variant of the A** algorithm is used to search the 
hypothesis tree. Track spawning is managed. Resolution is 
accounted for in, track-to- track association by allowing a contact 
to be shared by multiple tracks. Stars are eliminated using a star 
catalog. 

6.2 TRACK MAINTAINENCE - STATE ESTIMATION 

Triangulation is performed to initiate the precision cartesian state 
vector. ECI coordinates are used. State vectors are combined when 
a merge action is declared and the covariance matrix is adjusted 
accordingly in a manner similar to PDAF. The filter is an extended 
Kalman filter. The passive ranging algorithm is an iterated maximum 
likelihood algorithm. Biases ja.re accounted for by adjusting the 
system noise matrix or measurement matrix. 

6.3 TRACK MAINTAINENCE - TRACK PROMOTION /DEMOTION 

A Bayesian multiple hypothesis approach is used in data association. 
The hypothesis scores are used to generate track quality scores for 
each track. The track quality scores account for clutter density, 
new track density, miss detection rates and cumulative chi square 
scores. Thresholds are used to ^promote or prune tracks. 
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7. TRACK FILB MAINTAINENOE 

Essentially all zrack data is maintained in the Bayesian network. 
The network icaintains singleton object tracks, singleton CSO tracks, 
singlctcn gro^ui* tracks and forr.^^tious cf these tracks. Critical 
deployment points are also maintained. 

8> OUTPUT TO BM/C3 AND USERS 

State vectors for threatening targets are output to BM/C3. 
fl. COMPUTATIONAL REQUIREMENTS 

Computational requirements are discussed in Reference 2. 
10, CURRENT STATUS 

The algorithm is currently being implemented, in ADA, in preparation 
for an algorithm demo in the first half of 1989. The algorithm to 
be demonstrated will incorporate f rame*to-f rame and track-- to- track 
association, birth to death tracking and group to object tracking on 
RVs ajid penaids. Future demos will include PBV tracking. 

11 PERFORMANCE MEASURE k RESULTS 

The primary performance measures to be used in the algorLthm demo 
are track accuracy, track resolution, track continuity, track time 
to firm, track purity and track handover performance • 

12 REPORTS 

Pertinent information and performance data can be found in 

1. Technical Operating Report. Processor Development Report* CDRL 
008A9, Volume II. Processor Preliminary Development. Section III. 
Critical Algorithr Development. 31 July 1987. LMSC. Contract 
F04701-84-C-0102. (Secret) 

2. Technical Operating Report. Processor Development Report. CDRL 
0G8A9. Volume II. Processor Preliminary Development. Section I. 
Processor Architecture Development. 31 July 1987. LMSC. Contract 
F04701-84-C-0102. (Secret) 

3. Design Review Data Packaged SSTS Demonstration Design Review 
Data Package. Tracking Algorithms.. 12 October 1988. LMSC. 
Contract F04701-87-C-0093. (Unclassified) 
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ABSTRACT 

The algorithm presented Is intended to perform fxmctions of track 
initiation and track maintenance in the exoatmospheric ballistic missile 
defense scenario, assuming LWIR sensors aboard a set of mid-altitude 
satellites. The target density is effectively reduced by calling upon an edge 
tracker initially, which forms track files of edges of clusters. Track, 
initiation is then accomplished by referencing to those edge track files, 
assuming that targets in the same clxister travel in parallel. The 
iterative-least-square (ILS) filter employed in the track initiation process 
is itself initialized in a special manner in order to handle conditions of 
high target density and poor observability. Tracks initiated by two sensors 
are merged to provide precise state estimates to the extended Kalman filter 
that is used to carry out the track maintenance task. 
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1. CONTEXT 

The algorithm Is applicable to the mld-coiirse tracking problem. It 
contains two major components: track Initiation and track maintenance. The 
track Initiation portion basically follows the ftinctlonal chain of ~ ; - ~ 

monO' then' stereo and the track maintenance portion follows the centralized 
sensor - by- sensor approach • 

Inputs to the algorithm are angle measurements from either a single or 
multiple passive optical sensors. The outputs are track files and associated 
state estimates. 

2. NOTABLE FEATURES 

In the track initiation stage, track files are formed based on the 
assumption that targets within the same cluster travel in parallel. Reference 
trajectories are first established for edges of clusters using an edge 
tracking algorithm. The state associated with a track file is estimated using 
an iterative -least- square filter which is initialized in a special maimer. 

In the track maintenance stage , a centralized extended Kalman fil;:ter is 
used to continue tracks based on angular measurements from multiple sensors. 

3. SENSOR ARCHITECTURE & THREAT SCENARIO 

The algorithm can handle satellite •borne or probe-borne sensors. Sensor 
parameters such as resolution, sensitivity, measurement precision, bias, and 
detection threshold, are adjustable. 

Structured backgroimd clutter and background stars are assumed being 
removed prior to tracking. The only noise effect considered is the additive 
Gaussian noise. 

CSOs exist In the threat. Targets in the threat are all ball£stlc and 
non-manetjverable . 

4. SENSOR MODEL/PROCESSING 

A sensor functional model is used to generate angular measurements from 
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a scanning passive optical sensor. Included in the fimctional model are 
effects of random noise, bias, and CSOs. The CSO modelling is particularly 
elaborate. The scanning mechanism is also simulated. 

5. TRACK INITIATION 

Tracks are initiated from a cold start. An edge tracker is first used 
to establish track files for edges of clxxsters. All targets are then 
initiated in either a closely- coupled or a loosely- coupled parallel search 
mode, asstuoaing targets within the same cluster travel in parallel. Five 
frames of data are \xsually processed in a batch mode. Tracks are split 
initially and then pruned based on the nearest neighbor rule or the criterion 
of minimum residual chi-square. A complete trajectory is generated for each 
track file by calling upon an angle -only tracking algorithm. 

6.1 TRACK MAINTENANCE - DATA ASSOCIATION 

A simple nearest neighbor rule is used for data association during the 
stage of track maintenance. 

6.2 TRACK MAINTENANCE - STATE ESTIMATION 

An EKF is used to maintain tracks. The filter recursively updates the 
state estimate of each track file using angle measurements from multiple 
sensors. The ECI coordinate is adopted. Only white Gaussian noise is 
considered for the filter. 

7. TRACK FILE MAINTENANCE 

A track file maintains three categories of data: (1) the current state 
estimate, its covariance matrix, and the chi-square value, (2) status of: 
CSO's, and (3) intensity measurements from three color bands. 

8 . OUTPUT TO USERS 

The current algorithm does not explicitly provide outputs to tisers. 
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But, signature data associated with each established track file could be sent 
to discrimination algorithms, and metric data of all objects could be sent to 
BM/C^. 

9. COMPUTATIOHAL REQUIREMENTS 

The algorithm was Implemented in an IBK mainframe computer. The 
computational requirement depends on the number of targets and the number of 
frames to be processed. 

10. CDRRENT STATUS 

The algorithm has been implemented, tested, and run in a number of 
simulated threat/sensor scenarios. Currently, It is being Integrated with 
discrimination algorithms and radar tracking fxmctions. 

11. PERFORMANCE MEASURES & RESULTS 

Two scoring methods, the target oriented measure and the track- file 
oriented measure, have been used to evaluate the performance of the algorithm. 
The most critical factors which Influence the performance are the target 
density and the sensor resolution. 

12. REPORTS 

1. C. B. Chang and L. C. Youens, "An algorithm for multiple target tracking 
and data association," Technical Report 643, Lincoln Laboratory, M.I.T., 
June 1983, DTIC AD-A1131313. 

2. C. B. Chang, K. P. Dunn and L. C. Youens, "A tracking algorithm for dense 
target environments," Proceedings of American Control Conference, San 
Diego, CA, Jirne 1984. 

3. M. J. Tsai, L. C. Youens, and K. P. Dunn, "Track Initiation in a dense 
target environment using multiple sensors," Proceedings of SPIE 
conference on Digital Signal Processing, Association and Tracking of 
Point Source, Small and Cluster Targets, Orlando, Florida, March 1989. 
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Abstract/Summary 

The algorithm.is applicable to two passive sensors tracking individual objects in the 
late midcourse phase. It con-elates observations between sensors "up front" -before 
attempting track initialization. 



1. Context 

The processing chain is Type Iv (?). The algorithm assumes a cold start. It 
includes track initialization and maintenance. It is applicable to the late 
midcourse phase, when objets are resolvable but before any atmospheric 
influence at reentry. It assumes two passive sensors in exoatmospheric ballistic 
trajectories are tracking the target objects during the same time-period. The 
outputs are 3-D target track files. 

2. Notable Features 

The algorithm is much faster than altematives we have pursued or heard about. It 
initiates track files relatively rapidly, and its conreiation mistake ratios are low. It 
also con-ects for bias and navigation errors in the observers. 
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3. Sensor Architecture and Threat Scenario 

The threat consists of thousands of objects in ballistic trajectories. The inddence 
of unresolved objects is low. The sensors are passive scanning infrared sensors 
deployed in exoatmospheric ballistic trajectories. — 

4. Sensor Model/Processing 

The Lincoln Laboratory sensor functional modeli is used. The sensor is a 
generic scanning IR sensor. The frame time is variable. One wave band is used. 
The sensor-signal processor supplies line-of-sight, LOS-rate, observation time, 
and estimated measurement accuracy to the data processing function. Attitude 
bias errors arid navigated position and velocity errors are included. Detection is 
based on signal to noise ratio. 

5. Track initiation 

Cold-Start Track initiation is performed. The track files initiated are 3-D. A scan- 
to-scan congelation is done on two frame cycles of data. These paired 
observations are interpolated to a common time. The interpolated angles are 
matched between sensors with a "pseudo-elevation" angle^. With this correlation 
four observations are now correlated. These are used to initiate a track file with a 
square-root information filter. A Householder algorithm is used to solve for the 
state vector. Multiple hypotheses are not generated. On the first two scan cycles 
scan-to-scan correlation is performed by estimating line-of-sight rate from the 
average LOS rate taken from a sector by the signal processing. Once one or 
more track files have been established the line-of-sight change between two 
frame cycles is estimated by assuming the object in track con'esponding to the 
Observation nearest to the obsen^atioh in question is at the same range and has 
the same velocity as the object corresponding to the observation to be correlated. 
Sensor-to-sensor correlation is performed by matching the pseudo elevation 
angle generated from the interpolated line-of-sight measurements taken from two 
frame cycles of data. 

1 C. B Chang and K. P. Dunn, "A Funptional Model for the Closely Spaced Object 

Resolution Process". Lincoln Laboratiory technical Report 611, 20 May 1982. 

2 Blackburn, T. R., "A Practical Correlation Test for Cooperative Passive Optical 
Sensors", AIAA Journal of Guidance, Control and Dynamics. Vol. 6,. No. 1, Jan-Feb. 
1983. p. 62. 
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6. Track Maintenance 

6.1 Track Maintenance 

Data association in track maintenance relies on the apriori prediction of the next 
obiservatlon inherent in the stmcUire of recursive (Kalnrian) filters. New 
observations "nearest" to the prediction are chosen for updating the track file. 
The nearest distance measure is weighted with the filter enor residual 
covariance. 

6v2 Track Maintenance ■ State Estimation 

A U-D mechanization of the Kalman recurelve filter is used for state estimation in 
the track continuation phase. Coordinates are earth-centered inertia!. The filter 
utilizes estimates of the measurement noise level generated in the signal 
processing. The object tracking filter utilizes a 6-element state vector, containing 
object position and velocity. A separate estimator is used to estimate the relative 
elevation angle bias between the two sensore. 

6.3 Track Maintenapce - Track Promotion/Demotion L 
If observations fail outside 3-sigma bounds on the error residual four times in a 
row the track nie is dropped. The dropped file is expunged from the inventory of 
track files. 

7. Track File Maintenance 

individual or assumed-individual objects 7-D state vectors are maintained in track 
files, along with estimation uncertainty covariance matrices stored in upper 
triangular square-root format. 

8. - Output to BM/C3 and Users 

These questions havent been addressed. The data of item 7 above is what's 
available now. 

9. Computational Requirements 

WE have done most development woric on a CDC-860 machine which is rated at 
about 2.1 mega-flops. It took 4.28X10-3 sec per object per scan cycle. The 
computing load increases lineariy with the number of objects and frame cycles. 
This speed is over tOO times faster than the type -11 scheme we developed eariy 
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in the decade. Recently we have rehosted it on a VAX computer. It functions weir 
almost entirely in single precision on this 32-bit machine. 

10. Current Status 

The algorithm is in the late conceptual development stage. There are still areas 
for refinement. It has been tested and debugged mnning against threats with 
slightly over 600 objects in them. It has not been integrated with bulk filtering for 
star and false alarm rejection, etc. 

Performance Measures - 
The performance measure concentrated on has been correct obsen^ation 
assignment. The figure of merit is the difference between the number of 
observation assignments to a file and the number of times the source object 
assigned to a file the most was assigned to it. The difference between these two 
numbers is the number of mistakes made. Our scores are running better than 
99% against threats fumished us. The tracking filter estimation accuracies reflect 
Cramer-Rao bounds when files are free of assignment mistakes. j_ 

12. Reports 

IRAD Program Description 1-221, Optical Sensor Technologies, Appendix-A, 
MDAC Report MDCQ0931-1, 1989 Independent Research and Development, 
McDonnell Douglas Astronautics Co., Huntington Beach Division, March 15, 
1989 
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Abstract and Summary 

The Cluster Map Tracking approach is designed to treat the problenis of muld-sensor, multi- 
precision, and multi-coordinate system trac&g within a ba^e management system. This problem 
is resolved beginnm with the splutioh to the transition of thnisted obje^ into clusters, and their 
subsequent transition into single object track files. The inx)cess by which this is accomplished is 
C2iiltd clustetniapping^ }^ the details of this concept are presented in the paper Cluster Map 
Tracking in a Battle Management System produced by this author for the battie management 
IAbgram*s design notes. 

The crux of the idea is tiiat there are patterns within a threat that are naturally fomied by die objects 
dispersed from the same PBV and which are heading toward the same target This collection of 
cpmoving objects is called a c^ter within this discussion and will serve as the keystone to the 
understanding of the data reduction and parallelization that can be gained from the clusto* mapping 
methodology. 

When angle measurecqents are inade of the objeas in a cluster, the patterns formed can be traced 
from one scan of the associated sensor to another by using the previous scan's two-dunensional 
map of the cluster as a pattern for the next scan's data associations. This single technique leads to 
a battle system that is capable of learning as it evolves from single ppiiit source clusters, just 
b&thed mm a PIB V, to frilly developed clusters with all objects resolved. This learning 
phenonomena makes the battie manager trackihg system adaptive by nature and robust under a 
variety of different threat scenarios. 

These patterns of angle data measurements are stored from one scan to the next as cluster maps. 
Each clustier has a list of the associated cluster maps, one for each sensor that views die cluster. 
Notice here thiat die explicit problems of miilti-sensor correlation are not tackled until the cluster 
maps have been extracted. This interlocking of the cluster to its source data permits each map to be 
used as a filter to remove the cluster's new n^asurements from the freldrof-view for that sensor. 
The extracted data set, if determined to be valid, can then be used to replace the existing set as the 
new cluster map. This new map continues into die extraction process on the next scan, with each 
new inap definmg more of die cluster's detailed structure. Eax:h time a cluster map is updated, die 
associat^ cluster centroid state is updated via extended Kalman filtering. Where the data used in 
the iipdate is a calculated pseudo-measurement created from the collection of angle measurements 
that define the cluster map. These updates provide a mechahism by which the cluster centroid state 
is improved cohtihuously from the many different sensor views. 

By using cluster maps in this manner, the multi-sensor correlation problem is first indiiecdy solved 
at the cluster level, aiid then later at the object level As each cluster map is updated, and if that 
cluster is deemed by the batde system to be subject to object tracking, the estimated measurement 
data is passed on to the objea trackihg algorithni. Because the cluster mapping process has already 
divided the threat into clusters, this process can be accomplished in parallel witii specialized . 
hardware. If it is the first time that data is being passed to the object tracldng system, then an initial 
state vector and covariahce is formed for each object within the clustv. This is done through the 
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cmtrpid state* the angular separation of the objects within the cluster from the centroid, the cluster 
spread and si^:ead rate, and the angular measurement accuracy. 

Once an initial state has been sent, the object tracking algorithm uses that state as the initial track for 
that object, lifKlatihg it as tiine progresses widi the new measurements as tiiey appear. Splitting, 
merging, and other object track phenomena are handled within the object tracking environment, but 
are always restricted to die domain of the clustcrrOh subsequent updates of the cluster map, the 
associated object angle measurement data is sent to the object trackmg algorithm, where it is 
correlated with the hmited list of tracks associated with the cluster in question. This process, as it 
was conixived, will automatically generate multi-sensor cbrrelated tracks for each objea within the 
cluster without actually performmg multi-sensor correlatioiL 

1.0 Context 

The cluster mapping approach is designed to cover the entire domain of a batde scenario, from the 
details of thrust^ toost-phase tracking to the sp^niics of engagement-level tracking and the 
reporting of the needed in-flight guidance updates to the kill vehicles. Because the concept is 
applied over the entire scenario of the battle, from the birth of objects in die sensor view to their 
death fiom a kill operaticwi, it is necessary to specify the details of each phase to appreciate the 
context in which the method can be applied 

For die boost-phase, the approach requires that the thrusted objects (boosters and PB Vs) be 
tracked with a minimum of 9-12 state elements being developed maintained the PBV. This is 
necessary so that the deployment seed points for the lethal objects can be calculated via thrusted 
flight algorithms. TTie main requirement to track these thrusted objects is two-fold. First, by 
knowing die details of die PBV's track, it is possible to determine die approximate location at 
which a cluster of objects was released. This then gives the cluster a tie-in to die PBV dia^ 
*birthed* it Second, with knowledge of die booster track, die launch point can be determmed. 
This allows intelligence (a priori) iirformation to be used in die determination of the vehicle type 
(SS18, etc.). Li addition, widi explicit knowledge of die booster track, it is possible to make an 
association between die PiBV and die booster, thus providing a relationship between all deployed 
objects in b. cluster to their parent bboster and suspected payload. 

The details of the boost-phase operation are more critical to the batde planning, engagement 
planning, and threat assessment aspects of a batde manager than diey are to die tracking systenL 
TTiis is because it is still possible to create an initial track for a cluster without regard to its parent 
booster. This does imply a much more complex system than is currendy envisioned for a batde 
manager. For this reason, the details of how the a priori information should be used in a batde 
manageinent system are not discussed in diis paper. ; 

The contribution to tracking given by the boost-phase is that a cluster's centroid state can be 
initially estimated based on die PBV state at the time of deployment This warm start approach 
greatiy simplifies the process of track initialization. If PBV association is not available, then die 
standard method of track iiutiatioh is pCTfomied, whe^ sev(Eral frames of sense* measurement data 
must be gadiercd before die track can be initialized. In doing dus however, the information 
spedifying die lethality content of die cluster is unknown. 

In die post-boost phase, the boosters, PB Vs, and clusters are in a state of transition, widi die 
boosters eidier completely or nearly buriied but, the PBYs in the process of deploying dieir 
]^aylbads, aiid the clusters in the process of being forrnoL In this phase, the battle management 
system must begin to deal with ^e multi-sensor correlation problem. First, all tracking performed 
on the dxrasted objects must eidier be handled duough a du:usted track Kalman filter, where die 
d^ iiiput is eidier the raw angle measure ot pre-processed six-state estimates. In eidier 
case, the tfamsted states must be updated to maintain the track on these objects. In addition, die 
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sensor subsysteois in the post-boost phase will begin to see clusters emerging from their associated 
PB Vs as they aie seedbd along the fly-out path. 

These cluster reports arrive in the form of single object angle measurements or CSO 
measurenibnts. in die later case^ the measurements report ensembles of objects that are too closely 
spaced to be resolved by the sensor signal prcKessing system. No assumption is made here about 
the resolution of the detectors or the resolving capability of such a system^ but a few well-known-^ 
rules should be met First. t|ie system shoulilbe capable of resolving CSOs at the level of 4-5 
obje^pNerncm^ Second, the infortnation 

CSOs shoidd include qualiqr of fit. estimated number of objects, and spixie amplitude estimate or 
eiiisemble to re^sent the intensity of the clurnp. A system using sensor technology limited to 
sin^e and dual object tesoluidbn. with the reriiauider regarded as a clump, will certainly lead to 
system-wide failures (not only in tracking). 

The angle measurements arriving from the sensor must also be specif This 
fcmn is dictated by what the batde management system intends to do with the data in question, and 
not by the simplest method for any one scnspr subsystem. Since the main goal is to merge the 
multi-sensor views, eistablish viable tracks, and use the collection of iiradiance data to discriminate 
the lethal objects, the system must be cast into a conmion frame of reference. This is accomplished 
by specifying, as a requirement, to each sensor subsystem, that it render all of its angle 
nicasurcments and all of its state calculations into a selected Earth-Ceritprcd-Inertial ^CI) 
ccordihate systOTi (no specific definition required; only that there is unique definition). This 
rotation to the global inertial frame is necessary for the batde manager, in order for it to eliminate 
the data's dependence on the sensor's jpointinjg. IMU. and other detailed sensor system operations. 
In rotating the data to this common finune, the sensor must correct for iill of its known variances. 
This includes changes that would naturally result from a non-inertiat coordinate system injdiose 
sensors that are jaSL inertially directed. 

Now. with the data arriving from each sensor in a common fran^ of reference, the cluster map 
tracking sy stem is equipped to deal with this data, whether it be in the form of an^e measurements 
or in the form of six-states. It is more productive and accurate, however, for the system to operate 
off of the transformed raw data, which has been unaffected by any on-board sensor tracking 
algcdthms. Raw data is imiKMrtant to the overall performance of the multi-sensor correlation 
problem, as data modified by tracking tends to in-build non-linearities diat cannot always be 
overcome in stereo-viewing problems. 

As the system moves into the midcoursc-phasc, the angle measurement data continues, but fiom 
other sensor subsystems. It is these sensors tilat tend to supply the required information needed to 
refine the assessments and improve the guidance updates or the interceptors (which are typically in 
the air at this point in the scenario). 

In the cluster msp approach to tracking, the raw angle data is honed at each step of die way by tfie 
cluster's 2D cluster map (one for each sensor view). This data reduction metiiod is based on the 
pattern matching of this 2D image to the subsequent data set, and then using the selected objects 
from that data^set as a replacement image, perpetuatin operation until its conclusion. In 
performing this pattern matching, the System naturally subdivides itself, at the cluster level, into 
parallel operations, gready increasing the possibility of a true parallel architecture within the batde 
manager. In addition, as each cluster map set is seliected. it is passed oh to die objea tracking 
subsystem, where again it can be operated upon, in parallel, to render die desired object track files. 

The overview of the cluster mapping approach is shown in Figure 1, Here the diagram shows die 
inherendy parallel nature of the approach, segmenting the sensor data pipeline into parallel paths 
and the clustered object pipeline into parallel padis. This diagram is designed to demonstrate the 
possible locations for specialized hardware in the cluster mapping approach. 
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The details cxf the cluster mapping method as applied to cluster tracking is shown in Figure 2. 
Here, the diagram shows how the methodology hones the data down from the fiill sensof s view to 
the cluster maps, and then to the new clusters: It also demonstrates how the data is parallel piped, 
by cluster, into the object tracldng system. 

The details of object tracking under the cluster mapping approach are shown in Hgure 3. Here, the 
process is demdnstra^ to be restricted to the (lata asscKnated yath all the sensor views of a single 
cluster. The mechanism for update is simple 2b cdnelarioh, by which the algorithm is capable of 
updating the tracks and simultaneously poform multi-sensor cbrtelatibn at the clustered objea 
level 



2,0 Notable Features 

the cluster mapping approach is capable of 

• Six-state tracking of clusters in a threat (over multiple sensors), 

• Correlating iiewly formed clusters to their parent boosters/PBVs, 

• Creating initial object track states, 

• Updating object track states, 

• Reducing front-end correlation by matching object data to clusters, 

• Reducing front-end correlation by pattern inatching objea groupings to 
existing clusters, and 

« Pdf orming multi-sensor correlation as an automatic consequence of object track 
correlation. 

3.0 Sensor Architecture and Threat Scenario ~ 

Tlie assumptions made in the cluster mapping approach to a battle management system are that the 
sensor has (tone its best to remove clutter, star fields, biases, etc. and that the sensor system is 
capable of resolving taort than just dual CSOs (typically 4-S object CSOs). 

No assumptions have been made about the threat density, configuration, or individual vehicle 
capabilities. The number of objects has been estimated to be above 100,000. 



4,0 Sensor Model/Processing 
The basic infomiation expected is: 

• Hie angle values associated with.the rotation of the azimuth and elevation into 
tiie cominon ECI coordinate system, 

• Hie 3-element symmetric covariance associated with the errors in that 
coordinate system, 

• llie time of measurement, 

• The irradiance data for the different colors of the focal plane array, 

• A status to estimate the nature of the measurement (single, dual, triple, etc.), 
and 

• The sensor^s unique ID number in the battie systeia 
5.0 Track Initiation 

There are two forms of track initiation in the clusta mapping methodology, cluster centroid track 
initiation and object track initiation, bi the case of cluster centroid track initiation, there are also 
two nlethods by which initiaHzation occurs, PBV-seeding and cold start 
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5.1 Cluster State Initiation 

In general, the prefentd method of centioid trade initiation is to seed the traclc from a PBV at the 
time of the clustti's rcleaBe. It is at this point in time that Ae PBV state can be directly transferred 
as the cenmnd state of the associated cluster. In assigning the PBV to the cluster, odier vital (to the 
battle mahkger) infonmdbn can also be transferred for use by otiier algorithms. The associated 
state initiation pitfcess is accon^lished by one of two means. If tfie PBV can be flown directly to 
the point (to within sonoe tolerance) at which the citister is first obsenred, then the PBV state can be 
direcdy used (with only nidnar titne adjustments to the state) as the cluster's centroid state. 

On the other hand, if the cluster detection is sufficiendy late, and the cluster cannot be associated 
witih itoy taown PBV track during its deployment cycl^ 

d^ed to poinit the tfaccbackofthe cluster's an^e position to an ap^ This 
cone is defined by projecting the state at the estimated time of the first deployment forward for a 
fixed len^ of time, and tiieh projecting the state at the estimated time of the last dq;)loyment 
forward a related fixed lengA of time, these two projections form a cohe-like structure, whose 
apex angle represents the deviation of the PBV over its seeding cycle, and whose height is related 
to an estimated rnaximum time-of-flight before tiie duster will be observed. If die clusto- lies 
within diis cone and can be uniquely assigned to a P'BV, then it is assigned to the associated PBV 
and its parent booster. 

If die cluster cannot be assi^ed to a cone ot if a sensor's fint viewing of a cluster is sufficiendy 
late to make assodating it with a PBV mappropriaie, then cold start initiation is used. In tiiis 
approach, all object data rejected frotn the basic pattern matehing correlation is collected as potential 
new clusters. Th6^ grbupings which cannot be correlated to an existing PBV track are tiien 
compared to the existing cluster centroid tracks. This comparison is made because it is pcwsible for 
a newly arriving sensor to view, for the first time, a cluster, that has been seen by anothersensor 
and is abeady in cluster track . 

If a correlation is made, then die associated track is updated as a track maintenance operatitm. The 
remaining object measurements tfiat could not be property corrdated are dien collected over 4-6 
scans of data. This is done on a per-sensor basis as diere is no stereo correlation applied to raw 
data. After it is determinexi diat a suffident number of data sets have been collected, wmdowing is 
used (di^ are only a few objects to consider) to eliminstfe obvious ntisconrelations, a 4-point fit is 
calculated and an initial state is estimated. Once each state is created, diat 6-state is re-checked 
agidhst possible centrdd states to determine if a misdassificarion has occurred. This re-checking 
operation is parfdrmed periodically to purge die system of remnant tracks tiiat really belong to odier 
clusters. 

If diere are several objects widiin grouping distance, dien these objects, and didr new state 
estimates, are used to fonri a new duster and a new cluster centroid state. Aldiough diis later 
mediod will create a valid dustdr track, it does not assodate die data needed from die parent 
booster to die cluster. This lack of detailed infOTmation from a cold started track will affect die 
perfdrmahce of the remainder of the batde manager. 

5.2 Object State Initiation . . , , „ 

The second a^iect of track initiation occurs at die object trackmg Ipvel. Here, die state estimates 
are initiated based on the centroid state estimates of die encompassing cluster. This is 
accomplished by first detamimng which of die cluster maps (or sensor views) associated widi dus 
cluster ttintains die largest nUrhber of objects. This map, by definition, defines die view of die 
cluster where most of die objects are resolved and was chosen because it will render die largest 
number of initial tracks. This supports die philosophy diat it is easier to dispense widi old tracks 
thaii create new ones from scrateh. 
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It is this view of the cluster that is used as the source definition of the objea states within the 
cluster. The object states are fdnned by 

centioid state to the time of ineasurement for the object of intdresL At this point, the centtoid state . 
has been rendered into the same anj^e space as the angle noieasiirement for the object The changcf 
in angl0 recpured to get fibm the cehttoid to the object is tiien calculated and used as a lever-ann to : 
detertnine the deviation in pbsiticmlrpm the centrpid to die obj^ .This deviati<m can only specify 
the two positioh elements that are p^pendictjlar to tfie line-of-sight firom the sensor to the centioid£ 
The thiid ppsidori is determihed as the range along the liiieH>f-siight and is tat^n to be the same 
i^ge as found from the sensor to the c^^ The vdocity estimate is defined to be the 

same as die centtoid's state velocity. - 

^th the six-state of an object in the cluster determined, the error in that estimate must also be : 
calculatcxL This error will represent the covariancc of the object state and can be used in future 
track maintenance filtering opaadohs. The error in the directioil perpendicu^ to the line-of-sight 
is determined by the error in the centrcnd state in that dirtbcdon and the angular measurement error : 
of the diata as defined by the sensor. The error in the range direction is determined by the error in 
the centroid state iii that direction and the estimated spipad rate of the cluster (a parameter calculated 
in the cliisto* mapping process). The error in velociQr is determined based on the covariance of the 
centroid state for velocity and the estimated spread rate of the cluster. 

Once the state for each object within a cluster is estimated, it is not re-calculated by the cluster track 
system. Also, newly formed object tracks, resulting from splitting, are seeded from the original 
state prior to die split, thus elimmating the rieed to perform cold-start initialization on any new 
object tracks. The rernainder of the bpexiations performed on this objefct are determinod by die track 
maintenance cycle. This method dif clustef-inducbd processing guarantees that the requir^ and 
resulting tracks used during the operati 
processing to be irnplemented in paraUeL 



6.0 Track Maintenance 

The maintenance phase of die track occurs at two levels, the cluster track level and the objea track 
level. In the cluster track maintenance, there are two possible problems that must be dealt with: 
the problem of simply updating an existing track file, and the probletn of adding a new sensor's 
view to an existing cluster's track file. This later problem is equivaleiit to mulri -sensor correlation 
at the cluster level and is dealt with explicitiy in the cluster mapping approach. This approach also . 
gurantees that the multi-sensor correlation problem never actu^ly arises at the object track level, 
where it occurs as a natural consequence of the correlation process. 

6.1 Data Association 

The association of the new incoming data to existing tracks can be broken into the association of 
raw angle measurement data to cluster tracks and the association of cluster mapped angle 
measurement data to objea tracks. 

In the case of clusters, the angle measurement data arrives from the sensor, is queued, and is 
eyenmally processed to detehxiine to which cluster(s) the measurmient may belong. This is - 
accomplished by rotating the cluster ceiitrbid state into die fratne of reference of die sensor at the 
time of the measurement In doiiig so, the centroid position is rendered into the same an^e 
coordinate system as the incoming data. A (juick coiriparison is made and a separate Vertical! 
linked list is built foir each cluster that an object measuremeht could be assfign^ to. In doing this, 
the system is oeating the list of possible matches that cah occur as a result of trying to correlate die 
data y^th a given cluster. This Ust is linuted to the extent of die hunibe 
dkta's immedi^ neighborhood (as deterxxiined by input). 
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Once each queued object measurement has been correlated to its list of possible clustos, the queue 
of data is scanmld. As it is sc^ed, each measurement is 'horizontally correlated with the other 
measiotments thkt are also attached to a oomm<» cluster. These horizontal tiucads are woven 
through the nod^ for the cluster unda scrutiny, building a scries of horzontontal lists of virtual. 
cliister inaiKS. On<x these U have been conipletedf^ 

are arbitrated to deteniiine a unique set of assignmm Themethod 
used to artttrate the c6iiq>arisons is called /MOem 

This technique begins with a trial measurement in the horizontal list of the associated cluster, and 
starts by associating the object measurement to the first element in the cluster's clustn map (Note:. 
Both the horizontal list ioid the cluster map are orgs^uzed in time didcr to reduce die mis-match 
possibilities). If a valid correlation is found, bias i^ removed and die next element is considoed. 
A new bias is calculated fiom the nearest neighbor Avithin a predetermined tolerance window. If a 
cortdation can be made, then the pattern of the cluster is adjusted by one-half die calculated bias. 
Once accomplished, die bias for inatching the diird element to its nearest neighbor is calculated, 
except that here the window is now smaller (owing to die concept diat the pattern should be getting 
closer to a correct match). If an object is found widiin diis smaller window, die pattern is shifted 
again, only diis time by one-durd of die required bias. Each succeeding entry in die list is 
considered in diis manner, widi smaller windows md smaller bias corrections (where die 
corrojtions applied are inversely proportional to die number diat have been previously correlated). 
This piocess cohtinueis unt^l either a predetermined number of matdies have been reached, die list 
is exhausted, or there is no matoh within one of the target windows. 

If die process is unable to find a matoh, die algoridim shifts to die next entry in die list as die start 
and repeats die above steps. At die end, eidier a mitch was found diat satisfied die window ^ . 
reduction criterion, no inatehes were found, or a list of partial matthes was found diat sansfied die 
criterion up to some number of matehes. In diis to case, if a nliniinum number of matches is 
satisfied, die algoridim will declare tiiis partial match list as die correct assignment list In dus 
case, or in die case of a list which did satisfy die criterion, die remaining matohes are processed by 
first removing all unique matches, isolating all objiwt measurcnoents or clusta map entries widiout 
a inatoh, and dien performing a chi-square optimization on die rttnainihg mutiple matohes to 
(qitimally resolve &eir multiple correlations. 

If dierc exists a situation where multiple horizontal lists of object measurements are uniquely trying 
to matoh to a sin^e cluster map, dien diat cluster map is scrutinized for die possibility of cluster^ 
splitting. One indication of diis might be die lack of any object measurements near die calculated 
centroid. This can occur if, for example, a bimodal population forms and die centroid is calculated 
to be in die einpty region between die two populations. In dus case, die system assumes a spht has 
occured and dien uses die centroid state, rotated appropriately to each of die new centroid 
positions, to ^nerate die hew states for the ^lit clusters. 

Once a unique assignment is made, die list of assigned object measurements is redirected fiom die 
object measurement queue to die cluster map pointer for diat sensor and cluster. In doing so, die 
space allocated to the existing cluster map is fteed and die existing pointer is set to point to die new 
list, which how bebomes die new cluster m^. Finally, die data associated widi die new cluster 
map is nlerged into a piseudo-measurement iuid diis is used to update tite cluster trade with a 
extended Kalman filter. This process of update arid rq)lac!enaent continues for all new objert 
measurements. 

When updating an existing object ttack, die data must first be correlated to die existing objea 
tracks. This is accomplished very quickly and in parallel (for each cluster) in die cluster mapping 
approach to trackihg. First, because die tracks for the objects widun a cluster are initiated based cm 
die object data in die cluster (see section 5), it mearis diat die system will never start a cluster witii 
track files diat ccmtain data outside die cluster. Thus, each cluster is independent of the odier 
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clusters and can be processed as separate data sets. As a clus^ map is updated (rq)laced) by the 
new object nieasuiements, its data is collectively processed as at set of possible coiielations to 
objea track files. 

Tn armmnvvtote this, the object tracks associated with this clustcr are transformed to the 
associated srasor's view. This renders the object tracks into a 2D image directly comparable to the 
object measureinent daj^ in the cluster map. It is at this point fliat a simple conelation operation is 
performed, widi single ccnreladons b^ing removed and uinmediately {ttOcessed into updates. Next, 
correlations involv&ig equal lists of Object measmements match to equal length lists of tracks are 
optimized via any one of a nun^)er of q)timization algorithms (Muhkres, etc.) to render these short 
lists into a resolvable set these resolved equal length match lists are then processed into updates. 

The remaining possibilities are: more object measurements than tracks^ or more tracks than object 
measurements. In thie case of more object measurements than tracks^ the system must examine the 
nature of the measurements to determine if they represent the splitting of tracks. If this is tiie case, 
then a new track is developed to handle the split object and the old track is updated If more than 
one track is assigned to an object measurement, then the best match should be taken and the 
remaining tfack(s) should be left for later update atte 

Unda no circumstance does the object tracking system ever cold start a track from a sequence of 
stored, or batched, object ineasurenients. This capability may be considered as a last resort,- 
although it introduces siinilar heritage problems as discussed iii the cluster initiation problem. 



6.2 State Estimation ^ 
in the case of cluster tracks, once the data list has been used to replace the existing clustermap, the 
object data in the chx0r map is averaged (unweij^ted) to deterrnine a tiew pseudo-measurement. 
Hie pseudo-measurement covariance for this c^dilated noeasurement is a^ calculated based upon 
the measurement covariances of each object in the cluster map and upon the current spread rate of 
die cluster since last update. This change represents the possible en^ of mis-matched 
rneasurraientis into the data set. With the calculated rneasurernent and measurement covariance in 
hind, the cluster centroid state is updated accrading to a extended Kalinan filter. Here, the 
precision is maintained at 7 digits and a numeric error plant noise (Marquardt parameter) is used to 
weight the ctiagonal covariance elements. This helps inaintain the required convergence of the 
algorithm at tlus precision level and improves the solution steering that occurs in the extended 
filter. 

In the case of object tracks, the object measurement and covariance are direcdy applied according to 
a extended Kalimii filter. Here the precision is maintained at 15 digits to accomiiKxiate the 
requirements of die aiigle changes and coordinate transformations. 



6.3 Track Promotion and Demotion 

As die system maintains track, it handles the birth and death of a track differentiy depending on 
whether the track is defined for clusters or for objects. In the case of cluster tracks, die only 
ndeasure of the track's yalicfity is its existence, for a cluster is valid eyen if it contains no lethals. 
Quster tracks are ni£untained over the entire scenario of the battle to provide for the entry and exit 
of new sensors that may view the cluster. Because of this, a cluster track is riot terminated until.it 
is determined diat die centroki has in^acted. 

However, cluster tracks do undergo different stages of existe^iice. When a cluster is formed in the 
battfe manager, it is initially rnarked 'active' to indidate diat something is happening. However, 
this does not activate the cluster for the remainder of die batde manager. When a cluster is properiy 
associated with the parent booster and PBV, it is noarked ias 'available'. This marking is an 
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indicator to the battle planning systeoi that a cluster with valid heiiGige is availiOiIe for processing. 
Because of tidis, the cluster is also aissessed to detennine its potential inq)act-point and the 
associated possible lethality of die cluster mass. This assessed value is also stored in the track file 
for the planning operation. If the planning fiinctim determines that this cluster contains objects of 
interest to the battle manage, it (Erects the tracking syston to pickup object-level tracking on the 
objects within the clusuar. 

At this ppint, the tracking on this cluster switches from its pure 'data-gathering' mode to the 'object 
tracking' mode and creaites the initial object states fcff the object tracking operation. Upon 
subsequetit updates to this cluster, the cluster tracking system passes the accumulated object 
measurements (in the cluster map) down to the object tracking system, where they are correlated. 
The cluster tracldhg sy stem maintains this mode of (iteration on the clusters that it is tracking until 
the clusters inq>act die Earth. 

b the case of object tracking, the initial status of a track as defined by the cluster tracking system is 
that the track is active. Oiice a track is activated, it must receive continuous update nicasurcments 
frpin tiie cluster tracking systein. If iai specified (by input) numbd: of reports are missed, tiien the 
object track is considned for teiminatibn. Ah object track can only be terminated fiom missing 
data under two ^icumstances, it is not lethal or it has been killed. 

If an object track is determined to be lethal during its discrimination processing, and if engagement 
planning has deemed the object to be under engagement, then an expected kill should occur. The 
track is tiioefore forcibly maintained until such time as die expected kill does occur, or tfie 
expected kill tiiiie has elapsed If an expected kill of die letiial object does not occur, then a miss is 
declared and the old track file is deleted. If a kill did occur, then the track file will show a loss of 
data at tiic expected time of kill If tfiere is a data is loss, the track file is again deleted. The 
priinary reason for not deleting an engaged lethal track no longer deceiving measurement data is tiiat 
this track informatioii tnust be supplied to support in-flight guidance updates and homing views. 
This data and the associated interrogation of the track files would not be possible if die suspected 
object track and its cluster meihbers were deleted. This preventative medicine gready simplifies die 
data insmagement aspects of the weapon contrbller algorithms at the expense of some minor 
additional file rOaintoiahce on die part of objea traclang. 

If an object was not determined to be ledial, dien die loss of data is sufficient grounds for die 
tomination of the track. 

If at any time during die operation, a track which has not been updated for a specified length of 
tiine, recdveis hew angle measurements diatcan be used to update die track, dien die track is 
reactivated and die old discrimination data is purged. This is done to eliminate die possibility of 
old discrimination data polluting future discriminatim. 

Letiiality determination is die primaiy measure by which objects are deleted. In die batde manager, 
an object*s lediality is first determined through a bafteh discrimination process, where 1(^-20 
measurtinents (ftom any sensor) of iiradiahce d^ and state estimates are used to determine die 
lethality value. Subsequent lediality calculations arc performed on a 4-8 measurement cycle basis, 
depending on the diroughput loading, in an iterative manner sb diat new sensor inputs can be used 
to weight die lediality of die object as it progresses in die scenario. This lediality estimate is used 
as a mechanism to determine when tracks can or cannot be dropped. 
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7.0 Track File Maintenance 
The cluster track file contains: 

• The cunent six-state estimate of the centroid, 

• llie associated 21-eIeinent symmetn^ 

« The time at Which the cluster was first seen, 
f The time of the last update, 

• A status to indicate the clusters current mode of operation (normal, a member 
of the object tracking group, being lo^ 

• The as$essed value, its assigned cluster ID, 
r The ^sodated paitnt 

• Sdihe battle plamiing information maintained 

In addition, a sister file is used to maintain pointers to the data lists that make up the associated 
cluster maps for that cluster. This second file is accessed by cluster ED and by sensor ID, The 
actual cluster maps are maintained as part of the queueing mechanism used to manage the input 
sensor data to the battie management system 

The object track file contains: 

• The current six-state estimate of die centroid, 

• Hie associated 21-element symmetric covariance, 

• The time at which the state was initialized, 

• The time of the last update, 

• A status to indicate the object's current mode of operation (normal, under 

engagement, marked for deletion, etc.), 

• Its lethality estimate, ~ 

• Its assessed value, 

• Pointy to the discrimination data (which is kept in another file), 

• Pointers to the last discrimination state (which is kept in another file), 

• The associated cluster's ID, 

• The assigned track ID, and 

• The associated parent booster's ID. 



8.0 Output to BM/C3 and Users 

In general, the battie management system makes use of both the cluster and object tracking data. In 
the case of the clusters, the battie planning function has direct access to the cluster track files of the 
cluster tracking system. It is in these files tfiat the cluster centroid states arc maintained. These 
tracks are then used at the battie planning level to facilitate die global battie planning operations and 
platform-to-cluster assignments. 

The object tracks are used by the engagement planning and weapon controller pnxessing. In 
engagement planning, the letiial tracks are extracted and propagated to the estimated kill point, 
where propd: launch parameters arc generated for the assigned weapon (details not given hoe). 
Once the weapon has been launched and has reached burnout, the associated weapon controller 
extracts fixrai the object track files the guidance corrections tiiat must be specified to intercept the 
lethal object 

In addition, because the object tracking system maintains its tracks in a cluster organized manner 
Avhich is organized by cluster, the associated weapon controller can interrogate the file to determine 
which otiier objects are likely to be traveling witii die (designated object This list of comoving 
objects defines die homing view of the seeker as it moves into its terxninal phase of operation. 
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Thus, although no data is explicitly output to BM/G3, the raw data maintained in the track files of 
die battle manager, itself, ace used by die other parts of the system to detennine which cluster to 
shoot, which object in the cliister to shoot, how to correct for in-flight errors, and what the 
arrangement of the objects will be during die terminal phaise of the kill All of this is made possible 
by the data rechicdon technique in the cluster mapping zpproaclL 



9*0 Computational Requirements 

The shell for this tracking approach is written entirely in Ada and has been in^lemented on a Vax 
88Q0. Timing and paiormznot curves are unavailable and have yet to implemmted and 
deteimined on the tracking subsystenL the tmplicidy parallel nature of this approach must be 
tes^ to determine what gains can be achieved iii a specialized architecture that will allow the full 
painallel operation to be invoked 



10.0 Current Status 

Currendy, a tradi version of the cluster mapping methodology, has been implemented It performs 
the overall cluster birth, maintenance, filtering, and pattcm matching operations. However, track 
and correlation operations are not currendy implemented at the level of an algorithm, but rather are 
modeled at the truth level Future implementations will substitute various track filtering and pattern 
matctting schemes to define die bestpossiblie combination of these elements. 



11.0 Performance Measures and Results 

Individual track performance is based on the specific filtering algorithms used Overall system 
performance has not been evaluated for this system. ^ 



12. Reports. 

The equations, methodology, and the bulk of the discussion of details of the cluster map tracking 
approach are discussed in the paper Cluster Map Tracking in a Battle Management System (contact 
Gary Abercrombie, TRW, (205) 830-3302). 
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The MITRE Experimental Version (EVP) 
BALLISTIC TRACKER 



POINT OF CONTACT; J. A. Krajewski 

(617) 271-4547 
The MTTRE Coiporation 
Burlington Rd. 
Bedford, MA 01730 



ABSTRACT/SUMMARY 

The ballistic tracker accepts as input handover data £rom die boost-phase tracker and 
observajdoiis fironi any number of SSTS sensors. It then atten^ts to maintain PBV tracks 
and initiate RV, decoy, and debris trsicks. Ttite-dimensioal tracks are extrapolated into the 
future, and projected onto the SSTS focal plane. The "expected observations" are then 
associated with die actual SSTS observations using the Hiingaiian algorithm. The six- 
dimensipnal object state Vector is then updated with the two-diidensional observations 
using p extended Kalman filter. Objects, including PBVs, are modelied in die filter as 
ballistic Objects. Due to die relatively smiUI itiagnitu^ of PBV thrusts, a PBV can be 
cohsiderBd to be a ballistic ob^xit with ia sinall perturbation. Adiditional process noise is 
added to the filter to conrpensate for this approximiatioiL RV/decoy observations diat fail to 
be associated widi exiisting RV/decoy tracks are used to spawii new RV/decoy tracks using 
the nearest e»sting PBV track, if such a track is suffldehtly riear^ After each update, die- 
track library is examined to discover iaiiy tracks that should be pruned, basnl on a user- 
input number of failed updates. 

CONTEXT 

This algoridim belongs to the Type IV Sensor-by-Sensor Track Maintenance 
category described by the SDIO Panel on Gitical Isisues in Tra:cking. There is no cold start 
track initiation capability yet. Tracks are initiated by handover from the Boost-Phase 
Tracker and by spawning when a PBV deploys an RV or decoy. This algoridim performs 
track maintainence and spawning of RVs from PBV tracks. Inputs include booster 
handover data from the boost-phase ffackier, SSTS observritions, and a set of algorithm 
parameters. The outputs are object track estimates over time. 

NOTABLE FEATURES 

The ballistic tracker updates a six-state Kalman filter in Keplerian elements widi one 
sensor's azimuth and elevation data at a tiitie. It aicccnhodates data from any number of 
sensors. 

SENSOR ARCHITECTURE AND THREAT SCENARIO 

Tfie EVP can examine tracking performance variations with threat characteristics, 
sensor configurations, and sensor capabilities. The EVP does hot create threat data, but 
radier treats the threat as an input The EVP does model scanning sensors for bodi the 
BSTS and SSTS. 
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SENSOR MODEL/PROCESSING 



The sensor model is a scanning sensor with user-selectable scan rate. The model 
includes quantization error in the azimuth and elevation values output and also the effects of 
platform attitude eiror. The seiisor is assumed to be blinded if viewing any object within a 
plaiietaiy limb defined by the Earth's radius plus a user-iri^ut value representing the height 
of the bright atmosphere. The probability of detection is 1.0, provided that the object is in 
a position whose geometry permits observation (abbve the limb and within range). The 
probability of false alarm is 6. The EVP does hot attempt to remove clutter. The sensor 
model assumes that PB V and RV/decoy observations can be distinguished, but does not 
model or report brighmess to the tracker beyond that distinction. Consideration of a^)ect 
angle of the viewed object relative to the lineK)f-sight between viewer and viewed objea are 
ignored A single observation is generated when twd or more objects fall witiiin the same 
pixel. The SSTS has a user-input range Umit for RVs. 

TRACK INITIATION 

PBV tracks are initiated when the boost-phase tracker produces a handover. The 
boost-phase tracker estimates a bumout time, a burnout position, and a burnout velocity. 
The boost-phase tracko' also ufxiates the coyaiiance matrix in the boost-phase Kalman filter 
to reflect the uncertainty in die estimations. Hie ballistic tracker accepts die bumout time, 
poisition, velocity and covariance data, and initiates a track, the position and velocity of 
the estimated bumout point are transfoimed into Kepler dqiients. The covariance matrix in 
position and velocity is also transfoimed into a coyanance matrix of Kq)ler elements. ~The 
Kakian filter in the ballistic tracker uses Kepler elements as state estimation parameters. 
Thus, an estimated bumout state and covariance in the boost-phase tracker are transfoimed 
into state and covariance estimates in Kepler eleirients to initiidize die filter. 

When an R V observation fails to be associated with an existing RV track, tfie 
distances between tiiat RV observation and all existing PBV tracks are calculated. The 
closest PBV track is then found. If the distance between die unassociated RV observation 
and die nearest PBV track is less tiiari a user-defined direshold, then the PBV track is used 
to spawn a new RV track. A copy of the PBV track parameters is made and entered into 
the track file. This copied PBV state is then updated widi the RV observation, and d»e 
original PBV state is updated with its associated PBV observatiwi. Thus, die PBV track 
has been split into a continued PBV track and a new RV track. 

TRACK MAINTENANCE - DATA ASSOCUTION 

SSTS data are associated widi ballistic tracks in a two-dimensional sense. The 
diree-dimensional tracks are projected onto die SSTS focal plane and, using a scanning 
sensor model, a set of expected obisenrations is calculated. These, together widi die acnial 
observations, are soned by iziinudi. The azimuth values in die list are examined to check 
for any gap in azimudi that may exist di^t is largo: diiaii a dueshold. If such a gap exists, 
dien the two groups of expeaed and actual observation 6ti (ndiier side of die gap are treated 
: completely separately. ITiis is done to take advantage of iuiy structure within the data, a^ 
to reduce on the pipcessing time required for association. In diis way, potential 
associatioiis between an actual observation on one side of the gap and an expecxed 
obsenration on the other side of the gap are not even considered. 
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Each smmuth group of expected and actual observ 
elevadbn angle. Simil^ly, the elevations are exambed to detenxune if there are any gaps in 
elevation angle that are greiatter than a threshold Groups that are separated by sufGciendy 
large gaps in elevation angle are then treated separatel^^ 

A distance matrix for each separate group is calculated Each row corresponds to 
an expected observation, and each coliiinn cone^nds to an actual observation. Each 
entry in the matrix is the angular distance between the corresponding expected observatira 
arid actual observation. The Hungarian algorithm is then performed upon this matrix to 
determine the best set of associations. If, after an input ntimber of iterations, the 
Hpngaiian algorithm fails to detennine the cbmplete solution, die subset of the solution that 
has been unambiguously associated is kept, and the remaining objects are passed to a 
grody-type algorithm to compl^^^ After this process is complete, the 

asspdations are checked against a thre If the angular distance between the expected 
obsemtion and the actual observation are greater than a threshold the association is 
suppressed and the actual observation and tiack are considered to be unassociated in this 
scan. 

An option exists to bypass the Hungarian algorithm, and rely exclusively on the 
greedy algorithm to perform the associations. 

PB Vs and RVs are treated separately. First, RV observations are associated with 
existing RV tracks. Thien, unassociated RV tracks are compared to existing PBV tracks to 
determine which tracks should be spawned. Hnally, PBV observations are associatedwith 
PBV tracks. 

TRACK MAINTENANCE - STATE ESTIMATION 

The estimation process.utilizes an extended Kahnan filter using Keplerian elements 
as the state variables. The six state variables are updated using two-dimensional 
observational data from one SSTS at a time. Thus, azimuth aiui elevation data update the 
six Keplerian elements. 

All objects are modeled using ballistic trajectories. To compensate for the fact that 
PBVs are not iuly ballistic objects, more process noise is added to the filter for PBV tracks 
than for RV tracks. We believe the ballistic formulation for PBVs to be a good one; since- . 
PBV maneuvers are small compared t6 the total velocity of the PBV, die PBVs trajectory 
is nearly ballistic. 

TRACK MAINTENANCE - TRACK PROMOTION/DEMOTION 

Tracks are not promoted and demoted per se, but rather are all considered equally 
valid until pnmed. Each time a track fails tt) be updated, even though it has been calculated 
that the track should have been visible to the SSTS, a counter in the track file is 
incremented. If the track fails to be updated for ah input numbor of successive scans, the 
track is prqned, and no longer is consider^ in any future associations. If a track is 
expected hot to be visible to the SSTS, it is hot included in the asisociation, and its failed 
ufKlate couh^ is not increinented. Unclassified test cases with the current SSTS 
architecture have shown periods of time when some RVs are not seen by any sensor. 
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TRACK HLE MAINTENANCE 



The foUowing informanon is kept in the track fU^ 
Cuirent state estimate: Keplerian elements: 

• Semi-major axis 
Inclination angle 

• E ccentricity 

• MeanAnomaly 
Argmnent of Perigee 

• Argumentof the Ascending Node 

Six-by-six enor-covariance estimate in the state variables 

• Time of the most recent update 

• Object type (PB V or RV/decoy) 

• Number of successive missed updates 
« Track status (active or proned) 

OUTPUT TO BM/C3 AND USERS 

The output of the Ballistic Tracker to other BM/C3 algorithms and user includes: ^ 

• Number of tracks maintained over time 

• Object state estimates for each track over time 

• Lihbvation se()uences for each track over time 

COMPUTATIONAL REQUIREMENTS 

See CURRENT STATUS, below. 

CURRENT STATUS 

All algorithms described above have been coded. Test cases are currently being run 
and results arc being analyzed, the preliminary software for the ballistic tracker was 
written in FORTRAN and runs under the VMS operating system on a VAX 8700. In 
addition, this algorithm is being written in Ada as a component of the EVP Release 4.0 for 
the SDI National test Bed (NTB). EVP Release 4 wUl be submitted to the NTB in April 
1989. 

PERFORMANCE MEASURES AND RESULTS 
Perfoimance Measures: 

• Track purity over time (percentage of active tracks with m out of n consistent 
upd^ 

• Track accuracy over time on a track-by-track basis (distance between estimated 

state and troe object state) 

• SSTS sensor coverage over time on an object-by-object basis (which SSTS 
views a particular object during what time period 



B-68 



• Impact point accuracy over time on a track-by-track basis (distance between 
estimated impaa point and troe impact point) 

Results: 

A very small, unclassified test case (20 RV s) has been run using a notional 
unclassified SStS architecture, architecture, primarily as a 'maiden voyage*. Results 
indicate that SSTS coverage of RVs may be very sparse, even to the point where 0 or 1 
SSTS is observing a particular RV at times. The test showed that the Ballistic Tracker can 
function even with such limited data, but at a degraded level. 

Larger tests cases involving the SDIO's Phase One Test Specification (POTS) threat 
arc currently being performed. 

RELEVANT REPORTS 

In progress. 
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The MITRE Experimental Version Prototype (EVP) 
BOOST-PHASE TRACKER 



POINT OF com ACT: h H. Latimer 

(617)271-4553 
The MITRE (Corporation 
Burlington R(L . 
Bedford, MA 01730 



ABSTRACT/SUMMARY 

The Boost-Phase Tracker uses the inputs from two simulated sensors; associates 
the sensor reports for each scan into three-dimensional position estimates, and supplies 
those three-dimensional estimates as input to the tracking filter for associating with existing 
tracks or initiating new tracks. Three-dimensional position estimates are associated with 
existing tracks by a Greedy-type method based upon the distances between measured and 
expected points. Tracks are initiated from a single three-dixnensional position using a priori- 
assumptions about t^get trajectories. These algorithms have been coded, and have been 
tested using subsets of the SDIO*s Phkse One Threat Scenario (TSCB1*A). 

CONTEXT 2 

This tracking tool belongs to the Type IE (Sensor-to-Sensor then Scan-to-Scaii) 
track Initiation category and the Type in (SensOT-to-Sensor then Scan-to-Scan) Track 
Maintenance category as described by the SDIO Panel on Qitical Issues in Tracking. The- 
tracker performs track initiation arid track maintenance for boost phase. When booster 
targets are observed to bum out, the tiDol hands tracks over to MITRE's midcourse tracker 
(described separately). Inputs to the sensor riiodel axe model oruth target trajectories. 
Output consists of a file of sensor reports including sensor identity, two angles (azimuth 
and elevation), an epoch of observation, and (for analysis only) model truth target identity. 
These, togetiier with sensor orbits, are input to the tracker, which outputs tracking data 
consistirig oif target state estimates, covariances, and auxiliary (status) data plus track 
estimates to hand over to the midcourse tracker. 

NOTABLE FEATURES 

This tracker combines sensor data to derive three-dimensional position data for 
. . targets before associating those points to tracte^ tllis algorithniiriitiates tracks very 
rapidly, using only one triangulated three-dimensional position plus a priori assumptions 
about the target^s velocity and acceleration profile to initialize the nine-veaor state estimate, 
for the target track. 

SENSOR ARCHITECTURE AND THREAT SCENARIO 

The algorithm is designed to exaiiiine trackmg poform 
threats, b) using various senspr configurations, and c) assuming various sensor capabilities, 
(resolution, attitude error, and isampiing rate). Thb algorithm uses the sensor 
measurements from two sehsors of similar capaUlity and in distinctly differing positions.. 
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The two sensors are expected to be located so as to view all the targets to be tracked 
Although there are provisions to detect and remove singly observed measurements from 
further processing, the emphasis is on the estimation process for those targets that actually 
afford the possibility for stereo-optic tracking. The nieasurement em)rs and biases are of 
two categories: a quantization error iii the individual njbasurements in each of the two 
coordinates ~ azimuth and elevation - and an attitude enpr in the orientation of the 
platform which is constant over one frame, but which varies randomly with a given 
variance from one frame (scan) to the next It is assumed that all background clutter is 
removed prior to tracking. 

The boost-phase tracker was intended to use externally supplied threat data. 
Therefore, there are few if any assumptions about target kinematics embedded in the 
tracker. The tool has been tested with roughly 50 targets from one silo complex and also 
with approxiinately lOO targets from a handful of launch fields. 

SENSOR MODEL/PROCESSING 

The model used for the sensor is a scanning sensor that reports the two focal plane 
coordinates (azimuth and altitude) aiid the times of observations. Sensor resolution and 
altitude error are assumed to be variable. Measurements reported by the sensor are 
quantized based on the input resolution. Probability of detection is assumed to be 1.0^ 
provided the object is in a position whose geometry permits observation. The probability 
of false alarm is assumed to be zero. There is no provision to handle clutter. The sensor 
emulator will report only on the first occupant of a pixel, even if there are two or more 
targets present There is no representation of brightness other than the gross distinction 
between die brightness of a burning booster, and the relatively lower brightness of 
anything else. (See the description of the BaUistiic Tracker sensor formore discussion of 
the brightiiess distinctions between post-boost vehicles aiid RVs.) Since there is no 
brightness calculation, no relationship is assumed between the target signature (brighmess) 
and the aspea angle. Frame time is variable; tests have been iuii with fr^e periods firom 
oiie to ten seconds. There is no treatment of wavelength or closely-spaced objects (CSOs). 

TRACK INITIATION 

Track initiation begins when an object is first observed by two sensors. The 
preliminary data association (sensor-to^sensor) is identical to the method described below 
in TRACK MAINTENANCE - DATA ASSOCIATION. Because the sensor observations 
are asynchronous, it is necessary to allow for the target motion in the time interval between 
the two observations. The velocity of the target is estimated by identifying the target as one 
of the existing tracks, and using the track velocity estimate. When ah objea is first 
observed, there is no corresponding ttack, and therefore an assumption is made that the 
object is close to die surface of the Earth. The velocity of the object is taken to be that of 
the Eanh's spin at the point of intersection of the observation line-of-sight with the surface. 
Oiice the velocity is estimated, a reduced three-dimensional position estimate and 
covariance are used to start the track initiation process. 

A reference trajectory is used, together witii the reduced three-dimensional position 
estimate and covariance described above, to generate the initial state estimate for track 
initiation. This reference trajectory is derived fr6m a thfee-degi]ee-of*€reedom rocket 
trajectory model in use at MITRE for about four years. The trajectory model allows for a 
variable pitch aiigle profile, as weU as acceleration magnitude profile. A Newton*s method 
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iterative procedure is used to fit the reference trajectory to the three-dimensional target 
position. Right azimuth is assumed- the valu? of the reference trajectory velocity and 
acceleration at the point of fitting is used to initialize the track state estimate. In this way, a 
three-diniensional trajectory estimate ^i^ This pipcedure is perforaied serially for 
each threesiimensional position tiiat ha$ not been associated with a track, using the data 
association algc^thm describ^ below. This style of rapid track mitiation has the advantage 
that only one scan is needed to generate a track. The technique works well for those targets 
that are still relatively near the ground, and with degraded perfconiance for targets initialized 
at some altitude. 

There is no formal method of track promotion or demotion. Instead, the covariance 
associated witl> the track state estimate, plus the count of associated three-dimensional 
me<asurement points, indicates the level of confidence to be ascribed to each track state 
estimate. There is no specific provision for cluster track initiation, as in boost phase die 
closely-spaced-object problem is thought to be less severe than in midcourse. 

TRACK MAINTENANCE - DATA ASSOCIATION 

Two data association processes are performed as part of track maintenance: 
iassociation of data from two similar, but asynchronously scanning, sensors to geneiate 
three-dimensional object position and covariancie estimates; and association of those three- 
dimensional positions with existing tracks. ^ 

In this tracker, the association of images from two sensors takes advantage of the 
fact that a line-of-sight from a sensor at a known position can be mapped onto the focal 
plane of another sensor and intersect^ with images on the second sensor. This is because 
the positions of the two sensors andi the position of the target at the (assume, for the 
rnpnient, unique) epoch of observation define a plane. Iniagc data from each sensor, plus 
die direction to the other senior define a set of planes. The dihedral angles between tiiese 
planes and sotnt reference plane such as the plane containing the sensors and parallel to the 
Eanh's equator, form scalar keys which are ordered and used to associate the images from 
each sensor (see Figure I). When several images lie on, or nearly on, the same plane, it is 
necessary to invoke further, more involved logic, which we term disambiguation. 

Disambiguation is selectively performed on those groups of images not clearly 
associated by the plane-mapping technique discussed abiove. This technique relies on the 
ability of the track library to associate liiies-of-sight from sensors at known positions and 
epochs with tracks in the track file. When a line-of-sightis associated with a particulau: 
track, the height of the identified mick is intersected with the sight line to give a temporary 
ppsitipn estimate, which can then be used to form ai pscudb-irnage at the companion sensor. 
At the companion sensor, real and pseudo- images are then associated with each other 
using tiie angular positions in the plane defmed by the sensors confidendy and die targets. 
When these in-plane angles are too close to associate real and pseudo- images, then a third 
level of discriihination is invoked for those ambiguous asisociations remaining. This 
method uses both coordinates of the real and pseudo-images and a Greedy-type method of 
associations. 

The two-sensor data correlation algorithm requires lists of target observations. The 
assumption is that each list contains observations of the same targets. If the two lists are 
not of equal length, then the longer list must be examined to determine which observations 
in the longer list are without corresponding observations in die shorter list We denote the 
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algorithm that does this examination the "Discard" algoridim. The algorithm is very 
simple, and is bksed on the principle that the two lists of hinge angles are very similar. It 
looks for a diffdience in hinge angle exceeding a tteeshold, and ascribes that difference to 
an unequal population. 

Because of the asynchronous observations resulting firom the scanning sensor 
model, it is necessary to reduce the triangulated positions obtained ftom the correlated 
sensor lineis-of-^sight for the effect of the target velocity/ An estimate of the target velocity 
is obtained by associating a sensor tiiie-of-sight with one of the estimated track positions at 
the epoch of observation. The relative figure of merit used is the statistical distance (or chi- 
square) between the observation and its associated covariance and the track and its 
covariiance. With the estimated target displacement in the interval between observations, an 
adjustment (reduction) to the triangulated position is performed to yield an unbiased 
estimate of the target position. 

For one scan*s wonh of paired sensor data, a batch of three-dimensional target 
positions is finally obtained, and the process of associating these points witii existing tracks 
is begun. A Greedy-type algorithm is used to assosiate points with tracks, with upper 
bounds for associations observed. Thus there are/no multiple hypothesis techniques 
employed and therefore no problem widi any combinatoric explosion. 

We are not handling stars, false alarms or other fomis of stationary clutter. 

TRACK MAINTENANCE - STATE ESTIMATION 1 

The Kalman filter uses a nine-state model representing the target position, velocity 
and acceleration in Cartesian earth-centered inertial (ECI) coordinates. The model is 
represented as three ganged integrators, and the state is thfe current output of the integrators, 
there is no input to £e integrator prcKiucing the accelCTation, hence the model tries to 
produce constant acceleration. As the acceleration of a booster is certainly not constant, we 
allow the injection of process noise into the acceleration integrator to permit the variable 
acceleration of the target to be tracked. 




STATE VECTOR: 
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STATE TRANSITION MATRIX <D : 



1 0 0^ At 0 0 jAt2 0 0- 
0 1 0 0 At 0 0 jAt2 0 

0 0 1 0 0 At 0 0 jAt2 

0 0 0 1 0 0 At 0 0 

0000 100 AtO 
0000010 OAt 

0000001 00 
00 0 0 000 10 
0000000 01 

Noise measurements are assumed uncorrelated, and the measurement noise level is 
assumai constant. Crossing targets are not a problem unless one sensor sees the two 
targets in one resolution cell, in which case there is apt to be a phantom image for that 
frame. 



TRACK MAINTENANCE - TRACK PROMOTION/DEMOTION 

Tracks are not promoted or demoted. Instead, the covariance estimate and auxiliary 
infomiation (such as the accumulated count of associated observations) indicate the 
confidence to be assigned individual tracks. Tracks are teraiinated when a user-specified 
interval of time has passed without the association of observations to a track. Typically, 
the user might specify three to five scan intervals as the threshold. 

TRACK FILE MAINTENANCE 



Track files consist of the following: 

• Track file serial number (ID) 

• Track state estimate (9-vector) representing the three Cartesian 
components of each of the three quantities: position, velocity and 
acceleration in Earth-centered inertial coordinates 

• Track state estimate error-covariance matrix corresponding to the state 
estimates (9 by 9 matrix) 

• Time of track establishment 

• Time of last update 

• Accumulated count of associated observations (hits) 

• Root sum square (RSS) of last Kalman filter innovation value in 
position 

• RSS of the three current error-covariance matrix diagonal elements in 
position 

COMPUTATIONAL REQUIREMENTS 

See CURRENT STATUS, below. 
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The MITRE MULTI-SENSOR, MULTIPLE TARGET TRACKER 

POINT OF CXDNTACT: R. Varad J. T. McKcman 

(617) 271-4555 (617) 271-4546 

The MTTRE Corporation 
Burlington Road 
Bedford, MA 01730 



^ ABSTRACT/SUMMARY 

MITRE*s multiple-sensor algorithm performs sensor-to-sensor and scan-to-scan 
correlation to initialize and maintain tracks of ballistic missiles in boost phase. Three 
sensors are used in the current software implementation of the algorithin. 

As seen in Figure 1 , the two main sections of the algorithm are the multiple-sensor 
fusion and the scan-tp-scan coirelation and tracking, The first section fuses data from the 
stereo association of sensor pairs AB and AC. Riuige is determined finom the common 
sensor. A, once target lists from each pair are formed and then matched based on (1) hinge 
angles (the angle between each sensor-target-sensor plane and a previously defined 
reference plane), (2) in-plane angles (angles at each sensor in the sensor-target-sensor 
planes), and (3) estimates of the baseline ranges determined fix)m each pair. Rguie 2 and 
Figure 3 show the geometry of the sensor-target-sensor layout 

The second section, correlation and tracking, performs mainly track initialization 
and track maintenance. The key element to die tracing scheme is calculating rate of change 
of hinge angles, in-plane angles, and ranges. When these quantities are available for 
tracks, the rates of change can be propagated ahead to make predictions to be associated 
with the next set of sensor repons. The algorithm assigns a status to each track based on 
age and number of associations (see Section 5), and this status is the basis for 
upkeep/demotion and confirmation/deletion of tracks. Tracks are stored as state vectors. 

Since the algorithm's core processing consists of sorting, rather than mathematical 
optimization, processing requirements are relatively low, on the order of n log(/i), for n 
targets. 

1. CONTEXT 

The processing chain in MITRE's multiple-sensor algorithm most closely fits the 
general description of the **Staric then Dyhainic" chain as defined by the SDIO Panel on 
Critical Issues in Tracking. The two main functions of the algorithm are sensor-to-sensor 
and scan-to-scan correlation. The sensor-to-sensor portion includes the functions of two- 
sensor stereo association, unequal list matdhinjg, and three sensor fusion — fusion of the 
two sets of stereo associated data. The scan-to-scan portion's functions include track 
initiation and track maintenance. 

The algorithm has been implemented to operate on the SDI boost phase, but the 
concepts are directiy applicable to midcourse tracking with some modifications. The inputs 
are non-clustered observations from three scanning sensors. The outputs are 3- 
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dimensional, unsmooihcd tracks (state vectors). Provision is made for pre-ciustering of 
targets, if feasible. 

2. FEATURES 

The multiple-sensor algorithm makes use of a coordinate system that exploits the 
sensor-target-sensor geometry. Figures 2 and 3 show the three components of the 
coordinate system: the hinge angle, the in-plane angle, and tiie baseline range (as measured 
from a sensor to a target). The algorithm associates using two lists of targets sorted by the 
hinge and in-plane angles. The association technique, however, does not ctaploy a priori 
information, matrix optimization algorithms, or Kalman filters. 

3. SENSOR ARCHITECTURE & THREAT SCENARIO 

The software implementation of the multiple-sensor algoritimi has been tested using 
several threats, including a subset of the SDIO's classified TSCB-1 tiireat, and notional 
unclassified threats consisting of lOS and 130 ballistic missiles in both spike and ripple 
launches. (The 105-missile tiireat is launched from six complexes, while the 130-missile 
threat is launched firom two launch sites.) The sensor model does not treat clutter, sun, 
stars, or clouds. 

The sensor models we use are the notional Boost Phase Tracking System (BSTS) senson in a 
geosynchronous orbit. The three sensors are assumed to be centered over the Soviet Union, 
and separated by a variable amount (sensor separation between 40 and 70 degrees). Tests have 
been run using various representative estimates of sensor resolution and jitter. 

The unclassified ballistic missile threat follows a notional powered flight model 
from The MITRE Corporation's sirtiulation, while the classified tiireat follows the TSCB-1 
powered flight models. 

4. SENSOR MODEL/PROCESSING 

The modeling of the sensors is noi a function of this particular algorithm. 

The program can use eitiier scanning or staring sensors witii any user-specified 
frame time. Data are made available throughout powered flight Sensors report azimutii, . 
elevation, and time of observation. 

5. TRACK INITIATION 

Scan-to-scan association consists of two processes, namely, track initialization and 
track maintenance. A track is initialized for a given target when correlation can be obtained 
in hinge angle, in-plane angle and range for data from two conisecutive scans for the chosen 
sensor. Gates are established commensurate vnth sensor resolution and noise 
characteristics. The correlation is performed iridependentiy and sequentially in the three 
coordinates. This method of initialization is suitable for targets with relatively small 
velocities. For targets witii relatively high velocities, hinge rates obtained from non- 
simultaneous observations in the sensor-to-sehisor correlation process are used to 
tentatively initialize tracks. 
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8. OUTPUT TO BM/C3 AND USERS 



Results are in the form of output files with tracking statistics, from which track 
purity, range accuracy compared with model truth, and other quantities can be derived. 
(Track purity is defined as: for a given set of updates to a track, in this case six scan 
upniates, if five of those updates, in any combinadon, have identical model truth IDs, then 
the last scan update to the track is pure.) 

9. COMPUTATIONAL REQUIREMENTS 

The algorithm is implemented in Pascal on a VAXA'MS system. The memory 
requirement for 100 targets is about 300 Kbytes, and should scale linearly for larger 
threats. Throughput resulting from analysis of runs with 100 ballistic missiles is less than 
one second on the VAX 8700, and though not explicidy tested, the processing should be 
on the order of n log(/2). 

10. CURRENT STATUS 

The algorithm is fully designed, developed, and implemented. It has been 
debugged and tested with notional unclassified threats using as many as 130 ballistic 
missiles launched from two complexes, and a subset of the SDIO*s classified TSCB-1 
threat No performaince or hardware optimization has been attempted TTic concept used in 
this boost-phase model is applicable to midcourse as well; no midcourse model has been 
implemented, however. 

XL PERFORMANCE MEASURES AND RESULTS 

We characterized the efficiency and accuracy of this algorithm mainly by using track 
purity tests and by testing measured results against model truth for the calculated ranges, 
positions, and IDs. Tracking results are highly dependent on the threat characteristics. In 
general the ripple-launched ballistic missiles were tracked very well; the spike launch 
performance was poorer. For the case of the 105 boosters spike-launched from six 
complexes, the track purity (5 out of 6), ranges from 70 to 90 percent throughout. 

12. REPORTS 

Rajan Varad, The MITRE Corporation M87-73, Scalar Correlation 
Algorithm r Multi-Target. Multiple Senso r Data Fusion . December 1987. ^ 

Joseph T. McKeman, The MITRE Corporation Memorandum D44-M-337, 
'The Imjplemeritation of the Three Sensor Data Fusion Tracking Program in 
the Scenario Modeling System (SMS)", December 1988. 
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Project ; BMEWS Phased Array Radar Upgrade 
Organization : Raytheon Company 
Sponsor ; U. S. Air Force, Space Command 
Contact: Fred Daum (^-440-8734) 
Status : Operation since June 1987 

Reference ; B-5 Spec, for CPCI-2 (Software Requirements) 

Simulations : Real-time detailed simulation of threat, radar and C^ available 
In JOVIAL on CDC 170-865 machine. 

Purpose ; Process 3-D radar Inputs to estimate position and velocity vectors 
and support discrimination 

Data Processing Architecture ; General purpose digital computer (CDC 170-865). 

Threat ; Dense multiple target environment; objects can be very closely 

spaced; measurements are occasslonally unresolved In range; threat 

density Is approximately the same as the standard CSO SDIO threat for 
LWIR; ballistic or boost phase. 

Sensor : Phased array radar which measures range, elevation, azimuth and 
target amplitude. Data rate varies from 1 pps to 0.25 pps per 
cluster of targets with single pulses or pulse-pairs. Slgnal-to- 
nolse ratio Is adaptive, but nominally 15 dB. Resolution Is 
nominally 15 m. 

A Priori Information : None needed and none used. 

Implicit Assumptions : None 

Characteristics of Algorithms ; 

• 6-state Kalman filter 

• Multiple hypothesis track Initiation 

• Pulse-pair track Initiation and tracking 

• Nearest neighbor chl-square test for return-to-track correlation 

• Up-down chirp pulse-pairs are specifically designed to make track 
Initiation and tracking easy 

• Recognizes unresolved returns using quadrature monopulse test and 
prediction of target ranges 
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Project ; Ground Based Radar (GBR) 
Organization ; Raytheon Company 

Sponsor ; U. S. Army Strategic Defense Command (USASDC) 
Contact : Fred Daum (^iT'-AAO-SySA) 

Status ; Under development and testing In simulation 

Reference ; Software Requirements for Tactical Application Software for GBR 

Simulations : various pieces but no overall simulation currently available 

Purpose : Process 3-D radar inputs to estimate position, velocity and higher 
order rotational dynamics for discrimination 

Data Processing Architecture : General purpose digital computer. Immediate 

target machine Is von Neumann type with vector 
pipeline, but can easily be adapted to 
parallel processing of various types • 

Threat ; Dense multiple target environment; objects can be very closely 
spaced; dense chaff; all ballistic trajectories. 

Sensor: Phased array radar which measures range, elevation, azimuth, target 
amplitude and phase. Data rate is adaptive to local environment 
and Immediate requirements. Slgnal-to-nolse ratio Is adaptive as 
well. Resolution is adaptive to local environment. 

A Priori Information : Can use handover data (position and velocity with 

covariance matrix) or can search autonomously. 

Implicit Assumptions : None 

Characteristics of Algorithm ; 

• Three Kalman filters (6-state, 7-state, and 16-state) 

• Multiple hypothesis track initiation 

• Pulse-pair track Initiation and tracking 

• Nearest neighbor chi-square test for return-to-track correlation 

• Weighted average of monopulse angular measurements over multiple 
returns from a given object 

• Pulse-pairs are specifically designed to make track initiation and 
return-to-track correlation easy 
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SUMMARY DATA ON VELOCITY FILTER ALGORITHM FOR SDI 
DETECTION AND TRACKING WITH PASSIVE ELECTRO-OPTICAL SENSORS 



Space Computer Corporation 

2800 Olympic Boulevard, Suite 104 

Santa Monica, California 90404 

Projects and Sponsors 

1, Advanced Processing for Infrared Sensors (Army Strategic Defense 

« 

Command) 

2. Integrated Track Acquisition and Discrimination Concepts (MIT Lincoln Lab) 

3. Object Acquisition. Tracking and Discrimination via Bulk Processing 

(SDIO/ONR) 

4, Velocity Filter Approach to Boost/Midcourse Tracking (RADC) 

Contact 

Dr. William J. Jacobi or Dr. William B. Kendall 
(213)829-7733 

Status of Algorithms 

Under development and evaluation by Space Computer Corporation. The 
original concept is cue to Prof. Irving S. Reed at the University of Southern 
California [1]. 

Additional Information about Algorithm 

1. Purpose of Algorithm. The velocity filter algorithm performs a combination 
of signal-to-noise raiio enhancement and scan-to-scan correlation functions 
utilizing a "track-bef ore-detect" approach. It utilizes an input sequence of 
image frames (or equivalent data list) which have been processed for detector 
gain/offset correction, TDI, radiation hit suppression (gamma circumvention), 
etc. The filter output is an image (or equivalent data list) of objects 
having velocities within the passband of the filter. A bank of filters tuned 
to different vector velocities provides correlated object positions and 
velocities for track initiation. The vector velocities (and accelerations) to 
which the filters are tuned are derived from cross-correlation of successive 
input image frames. 

2. Data Processing Architecture Assumed." The velocity filter algorithm is 
most efficiently implemented by a special-purpose parallel processor attached 
to a general-purpose host computer. Studies of alternative processor ar- 
chitectures are being carried out by Space Computer Corporation. 
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TRACKING ALGORITHM SUMMARY* 



Submitter's Name: Kenneth Kessler 
Submitter's Company: Systems Control Technology 
Submitter's Phone: 415 494 2233 

Address: 2300 Geng Road. Palo Alto. California 94303 

Author: See Above 



Type of Algorithm: 



Streak Detection (Track before Detect) 



Title of Algorithm: 



Dynamic Programming Algorithm OPA) 



Sponsor: 
Developer: 



DARPA/NOSC 
Yair Baarniv 



ABSTRACT/SUMMARY 

Dynamic Programming can provide an alternative approach to the detection 
and tracking of dim targets using returns from a wide variety of sensors. The 
conventional (classical) approach to this problem is to obtain a detection of 
the target return from a single "frame", then attempt to associate this single 
return with a "track file" previously established on the object- 

Alternatively, the problem of detecting and tracking dim targets in IR 
mosaic-sensor imagery (or other sensor data) can be thought of as one of 
detecting and locating whole target trajectories inside the sensor's field of 
view over some interval of time, ^elected by the user, A batch processing 
algorithm is required in this case. Such an algorithm has the potential of 
obtaining superior performance oyer that of the conventional approach, because 
its performance (measured by its probabilities of detection and false alarm) 
would be based on the integrated energy of a target during its entire stay 
inside the sensor's field of view rather than on the short single-frame time. 

Conceptually, one candidate batch processing algorithm might use 
exhaustive search techniques. The multi-frame data is recognized as being 
three-dimensional (3-D), where the third dimension is associated with the 
frame number, or time of the frame. For the idealized white Gaussian 
uncorrelated noise case, the optimal detection consists of passing the 
"batched" data through a bank of matched-filters (MF). where each filter 
represents a single possible 3-D trajectory. 

The Dynamic Programming Algorithm (DPA) is a practical and feasible 
technique to replace the above exhaustive search technique over all possible 
trajectories and still get similar performance. The main idea is to define 
any possible trajectory as a string of straight-line short segments. Each 
straight-line segment is defined over a smalt number of frames which 
constitute a DP "stage". Exhaustive search can now be performed over one 
stage at a time, and the principles of dynamic programming used to "piece 
together" the results of one stage with that of the next. The number of the 
possible short segments, which is user defined, is relatively small so that 
this single-stage filtering is feasible. 
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performance is obtained when the target trajectories are separated by at least 
1/2 to 1/4 of a pixel and not parallel (in angle) over approximately 6 to 8 
frames of data. 

No specific clutter statistics are assumed, however the matched filters 
are designed for white Gaussian noise (other designs can be used). If the 
standard design is used, the frame data should be differenced, so that the 
background approaches this statistical assumption. 

The target model kinematics is assumed to be nearly straight, constant 
speed model. User inputs do allow, however, for the target to deviate from 
this assumption. Some curvature and acceleration is allowed by input specific 
parameters of the algorithm. Pfbcessing time can be minimized as the target 
characteristics approach those used in the kinematic model, by suitably 
modifying the input parameters. 

4. SENSOR MODEL/PROCESSING 

The algorithm was originally designed for a staring mosaic sensor, 
however in principle, the development can be extended to other mosaics such as 
step starers and scanners as well as radar. The frame time of the sensor is 
variable, but some preprocessing is assumed in order to obtain "apparent" 
target motion of approximately one pixel length per scan. The algorithm 
processes a single waveband of data per batch; additional wavebands can be 
processed independently, then merged to obtain a single track, if desired. 

Zrror parameters or characteristics are not explicitly modeled within the 
tracking algorithm. That is measurement noise, biases etc. are not modeled as 
in a conventional Kalman filter, which explicitly or implicitly models all 
process and measurement errors. The Dynaniiic Programming Algorithm implicitly 
assumes the following: that each pixel return is associated with the center of 
each pixel, that the measuremient noise is Gaussian and uncorrelated In space 
and time (or nearly so), and that the assumed target motion is "nearly" 
straight and of constant speed, the performance of the algorithm is optimal 
with these assumptions but degrades as these assumptions become less valid. 

The precision of the input sensor data has been tested down to 4 or 5 

bits of precision with very small loss of performance. The interface between 

the sensor and the tracker is required to be multiple frames of data, stored 

for subsequent processing. The numbers of frames used for processing are user 
defined and each frame contains a single energy return per pixel. Performance 

is improved if first pr second frames differencing is employed. In an attempt 
to decorrelate the background clutter and noise. 

The performance of the algorithm is largely dependent upon the signal 
strength of the target, background clutter and noise levels, numbers of frames 
processed, etc. The algorithm has been used to process HiCamp data at NOSC. 
and has demonstrated the ability to detect and track targets having a signal 
strength of 2 or so dB. with probabilities of detection of 0.9 or more, and 
false alarms of lOE-5 to -6. 
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The algorithm has been tested against both simulated and actual sensor data 
(HICamp I and II). The algorithm has not been married to a specific hardware 
architecture, rather, it has been used to process post-mission data where the 
specific sensor data has gone through a "pre-processing* stage, typically 
single frames differencing. 

11. PERFORMANCE MEASURES & RESULTS 

According to staff at NOSC. the OPA has been tested to detect and track 
dim targets having a SNR of approximately 1-3 in magnitude. Test cases have 
also demonstrated simultaneous track of 72 very dim, crossing targets, all 
within the field of view (simulated data) and actual data from the HiCamp I 
and II sensors. This latter set is with a variety of geometries and SNR's and 
targets. 

There have been theoretical performance bound computed during the past 
several years and contained in several reports and papers. The algorithm has 
not been sufficiently tested to determine under what conditions the algorithm 
performs poorly 

12. REPORT 

The most recent project report is: 

Dynamic Programming Algorithm (OPA) Analysis 
Systems Control Technology, Inc. 
Final Report under contract TLJ-5734-1 866-87 for 
Titan Systems, Inc. 

Authors: Kenneth Kessler, Michael Sutphen, Thomas Holmes, Barbara 

McQuiston 

November 1988 
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SYSTEMS. INC. 



3 April 1989 



Dr. Gabriel Frenkel, IDA 
IDA Tracking Panel 
Washington, DC 

.Subject: Panel Survey 

Dr. Frenkel: 

This letter is in response to the IDA Tracking Panel Survey. Before getting into 
the actual survey items, I want to provide you with some background on the 
Knowledge-Based Sensor Fusion (KBSF) program. It Is my desire to establish the 
proper context withih which the survey was completed. I have tried to fit a description 
of our system into your format, which is really designed for algorithm <jescriptions in a 
strict sense, so please bear these thoughts ih mind as you review the survey. 

As I indicated to you earlier, TITAN Systems, Huntsville is not in the process of 
designing tracking algorithms perse. The KBSF program was initiated by the SDI to 
determine whether knowledge / rule-based techniques could be applied to the 
strategic sensor fusion problem arena. Specifically, the intent is to demonstrate a 
rule-based approach to multisensor discrimination. As you know, it is not entirely 
possible to decouple tracking issues frorti discrimination issues, hence the KBSF 
program's interest is in hosting and evaluating "competing" algorithms (primarily those 
designed for platform-to-platform association and track maintenance). While some of 
these algorithms may have been designed for essentially the same track-related 
purpose, each may contribute differently td the discrimination process. Our goal is to 
build a machine that automatically selects the "best" algorithm for fusing multiple / 
single platform, multiple sensor track data, given the peculiarities of the scenario at 
hand. 

Of particular interest to our program are the methods by which uncertainties in 
the tracking functions (and hence uncertiainties in the radiometric / RF quantities 
important to discrimination) are characterized, propagated and dealt with at the Battle / 
Sensor Management level. We are experimenting with fuzzy techniques for 
determining "degrees of membership" in the set of valid tracks that is maintained by 
the system, and for uncertainty management in general. 
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5. Track Initiation 



Currently, a fuzzy variation on the Sequential Probability Ratio Test is employed 
to determine the degree of membership of each sensor report in the 
system-maintained set of "valid tracks". This degree of membership effectively 
weights each report. The weight is propagated through the system to the fusion and 
then to the discrimination process, where it impacts lethality decisions. The same 
procedure is used for track maintenance, in conjunction with an impact point prediction 
routine for both single and multiple platform association. RSS scoring is used. All 
target and sensor states are maintained in Earth-Centered Inertial coordinates. 

7. Tracic F iie Maintenance 

See attachment 1 . for a description of the Fused Data Frame which is utilized by 
the KBSF machine. 

8. Output to BfA I ca and Users 

See attachment 2. for an example display screen which shows 
Platform-to-Platform association results and discrimination results. 

9. Computationai Requirements 

The current system is implemented in a multi-SUN wori<station architecture. All 
KBSF components are written in C, However, algorithms written in other languages 
supported by SUN Unix may be called by the rule basses (the trajectory model used 
for impact point prediction is an example written in FORTRAN). This feature provides 
a powerful testbed capability. 

10. Current Status 

A multiple computer system has been demonstrated and effort is underway to 
enhance baseline algorithmic, rule - base and uncertainty management capabilities. 
Simulated sensor data has been used to exercise the system to date. Real sensor 
data is desired for system testing. 
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NAME: J. T. Lawson DATE: 14 February 1989 

COMPANY: TRW System Development Division, Huntsville Operations 
PHONE: (205) 830-3213 

ADDRESS: 213 Wynn Drive, Huntsville AL 35805 M/SHSV8-1414 
ALGORITHM TITLE: ADO? Scan-to-Scan and Track Algorithm Set 
SPONSOR: U. S. Army Strategic Defense Command (USASDQ 

DEVELOPER: TRW Huntsville Operations 
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message is sent Unassociated observations (those not assigned to a track) are then 
submitted to Track Association, which attempts to form strings of observations to generate 
new tracks over a period of several scan cycles. Successful strings of associations are 
reported to the Track process via New Track messages* Only unique assignments are 
reported 

The Track process is responsible for initializing track filters, updating track filters and 
making predictions for the next sighting of the tracked object Two-types of track filters are 
supported: angle-<mly and six-state. The angle*only filter is a five or six point cubic fit 
least squares filter in azimuth and elevation* The six-state filter is a Kalman filter in Earth 
Centered Bdresight Inertial (ECBI) coordinates. The ECBI system axes are parallel to the 
fixed boresight coordinate system of the sensor and are a fixed transformation from Earth 
Centered Inertial (ECI) coordinates. 

1. CONTEXT 

The ADOP algorithm set addressed the problem of single sensor onboard acquisition 
and track for a ballistic probe (e.g. GSTS) during die midcourse engagement phase. The 
sensor is assumed to hold a fixed orientation along an inertial boresight for the duration of a 
mission and executes a racetrack scan pattern with a slight overlap zone. Hence, each 
sensor *'£rame** consists of a top scan and a bottom scan within the sensor frame time. 
Detections fix>m sensor signal processing in three LWIR bands are processed by ADOP 
Measurement Pnx:essing algorithms (not discussed here) to perform color correlation, 
velocity estimation, irradiance calibration, bulk filtering of stars and scan smoothing, 
among other functions. The input to the Scan-to-Scan process from Measurement 
Processing consists of a set of object observations from the latest scan. Each observation 
record includes sighting time, azimuth, elevation, estimated azimuth and elevation rates, 
standard deviations for the angles and rates, and irradiances measured in the three color 
bands. Each record is marked as an isolated single or a CSO. For CSOs, azimuth and 
elevation extent data are also included. Outputs from the Track process to a ground-based 
batde manager conisist of threat leports that include object six-state and covariance, time of 
applicability, and estimated lethality. A threat report is provided every two to three frames 
for each object deteniiined to be lethal by discrimination. 

2. NOTABLE FEATURES 

The ADOP algorithm set was the first, to our knowledge, to incorporate a pattern 
matching approach to the midcourse track association problem. Moreover, both pattern 
matching algorithms and classical gating algorithms v/ert used in a complementary fashion 
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Figure 2: Scan-to-Scan and Track Tasks and Flows 
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of group members are written as a list to the PM_ReLMembers file. Degenerate groups of 
one member (isolated singles) are discarded and their scan file record numbers are written 
to the Unmatched Initiators file. (Altemadvely they can be retained for group matching and 
continue through pattern matching as aii option.) The Grouper also operates on the second 
scan (Frame 2) of the scan-to-scan cycle when it arrives. For this scan the outputs are 
written to the Master.PMjGroups and Master_PM_Members files. The Group Matcher 
task is then activated 

Group Matchen This task attempts to match Frame 1 groups to Frame 2 groups. The 
primary criteria are that the Frame 2 centroid lies in the predicted centroid gate and the 
extents of the groups are reasonably close. The comparative numbers of members in both 
groups could also be used. Where a unique match is found, information for both groups 
from the . . .(jroups and . . .Members files is written to the PM_Matched .Groups file for use 
in detection matching. Where there are either no matches or more than one match to a 
Frame 1 group, the scan file record numbers of the members of the Frame 1 group are 
written to the Unmatched_Iniriators file. When all Frame 1 / Frame 2 groups have been 
processed, the Frame 2 Detection Matcher is enabled. 

Frame 2 Detection Matchen This task operates in Frame 2 after the Group Matcher, hence 
its name. Its purpose is to match individual observations within the Frame 2 groups with 
those in the matched Irame 1 groups. An empirical multipass algorithm is currendy used 
to do the matching. This algorithm seeks to establish an initial match between one Frame 1 
group member and one Frame 2 group member, and use that match to narrow acceptance 
gates for matching the remaining members. Other techniques may give better performance 
with less effort. Exacdy which detection matching technique is best under what conditions 
is still an open research issue. The ratio of successful matches to group size is used as a 
criteria to accept or reject the match data. Engineering evaluations during the ADOP project 
showed tiiat when the ratio slipped below 0.34, die matches were usually incorrect When 
the matches are rejected, the scan record numbers for all members of the Frame 1 group are 
written to the Unmatched Initiators file. When the matches are accepted for the group, the 
individual matches art written to the Track Candidates file for use by the Frame 3 Detection 
Matcher. 

Frame 3 Detection Matcher This task operates during Frame 3 and attempts to match 
Frame 3 observations to the two-point candidate tracks established by the Frame 2 
Detection Matcher, Candidate tracks are processed by Frame 2 group. The group centroid 
and velocity are calculated and a Frame 3 extent gate is established. All Frame 3 
observations within this gate are processed in a manner sinoilar to that used by the Frame 2 
Detection Matcher. For candidate tracks which pass the match criteria, a predicted position 
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Track Arbiter : This task is the last stop in the track association chain. Its job is to ensure 
that only one track claims each of the constituent observations in the chain of observations 
that result in a new track. The arbitration scheme is simple, but effective - first come, first 
served. The Track Arbiter operates first when the scan observations from Frame 5 arc 
available, again immediately after the Frame 6 observations are available, and once again 
when all proposed tracks from the scan-to-scan cycle have been extracted and the system is 
ready to begin a new cycle* The Track Arbiter takes each instance in sequence from the 
Track Requests file, and p2X)pag9ites the candidate track to the most recent fr3me. If any 
ambiguities arise (more than one obseryadon per acceptance gate) the nearest to the 
predicted sighting is chosen. At this point the Track Arbiter checks all observations fiiom 
the previous scan files that belong to the candidate track being processed If any of the 
observations have been previously been taken the candidate is discarded If none of the 
observations have been taken, the Track Arbiter removes (or marics as taken) the 
observations from the previous scan files and sends a New Track message to the 
appropriate Track process. This message contains all needed data from the string of 
sightings that is needed to initiate angle-only track. 

At the end of the scan-to-scan cycle (six frames) the Scan Manago* task decides which 
ft^e will be Frame 1 of a new cycle. If pattern matching is to be employed, the new 
cycle must begin with the first frame following the cycle just completed However, if use 
of the gating algorithm alone is acceptable, the Scan Manager can designate either Frame 3 
or Frame 4 of the old cycle as Frame 1 of the new cycle. 

6. TRACK MAINTENANCE 

Most of the *'track maintenance" frmctions are performed in the Track process shown in 
Figure 2. The exception is the Track Continue task in the Scan-to-Scan process. 

6.1 TRACK MAINTENANCE DATA ASSOCIATION 

The Track Continue task in the Scan-to-Scan process associates scan observations with 
predicted object positions forwarded by the Track process(es), A simple gating technique 
was used witii gates computed by the Track process. If no observations lay in the gate, a 
Missed Track message is returned to the Track process. When one or more observations 
lay in the gate a Track Update message is returned to track with tiie angular measurement 
data for the chosen observation and the track id with which it has been associated When 
more than one observation is in the gate, the one closest to the predicted position is chosen. 
All predictions are processed sequentially and independentiy. Optimization across multiple 
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?• TRACK FILE MAINTENANCE 



The data maintained in the track file is dependent on the track mcxie. During angle-only 
tracking an observation history of up to fifteen sightings is maintained. For each sighting, 
the time, azimuth, elevation, az-sigma, el-sigma, class (single, CSO), and irradiances are 
niaintained In the precision track mode, the time, six-state, covaiiance, lethality, next 
predicted sighting, last measured sighting, and a radiant intensity history are maintained 
In both modes additional flags and counters are used to indicate track status. All track 
records are considered independent and there was no attempt to organize them by group, 
* cluster, etc. Several improvements could be made to support more sophisticated track 
continuation techniques, but none haive been implemented 

8. OUTPUT TO BM/C^ AND USERS 

Outputs from the Track process to a ground-based battie manager consist of threat 
reports that include object six-state and covariahce, time of applicability, and estimated 
lethality. A threat report is provided every two to three frames (i.e every twenty to thirty 
seconds) for each object determined to be lethal by discrimination. Each threat report is 
formed by extracting information frora the object's track file when the object id is passed to 
the Threat Reporting task by the Track Update task. 

9. COMPUTATIONAL REQUIREMENTS 

The ADOP software was implemented in Pascal to run on a distributed network of 
Honeywell ADOP processors. Each node consisted of tiiree 1-MIP MIL-STD-1750A 
CPUs sharing up to one megabyte of memory. A development test version, also in Pascal 
executes on DEC VAX hardware. The design goal was to support 300 to 400 objects per 
MEP witii one of tiie ADOP nodes whetiier it be Scan-to-Scan or Track. A 1000 object test 
case was executed in real time on a three-node ADOP network, but detailed internal data 
were not logged as it was primarily a hardware test A 65 object test case has executed 
successfully in real time on a single VAX 1 1/780. 

10. CURRENT STATUS 

Versions of the ADOP software are available at the USASDC Advanced Research 
Center (ARC) for botii tiie Honeywell and DEC VAX configurations. Ko significant 
improvements to the implementation have been made since 1985. A number of 
performance and capability improvements were identified under the APIRS program, but 
could not be implemented due to funding constraints. Some of these upgrades may be 
made under the GSTS program if an option is elected to convert the software to Ada for use 
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